Phylogenetically diverse diets favor more complex venoms in North American pitvipers.

The role of natural selection in the evolution of trait complexity can be characterized by testing hypothesized links between complex forms and their functions across species. Predatory venoms are composed of multiple proteins that collectively function to incapacitate prey. Venom complexity fluctuates over evolutionary timescales, with apparent increases and decreases in complexity, and yet the causes of this variation are unclear. We tested alternative hypotheses linking venom complexity and ecological sources of selection from diet in the largest clade of front-fanged venomous snakes in North America: the rattlesnakes, copperheads, cantils, and cottonmouths. We generated independent transcriptomic and proteomic measures of venom complexity and collated several natural history studies to quantify dietary variation. We then constructed genome-scale phylogenies for these snakes for comparative analyses. Strikingly, prey phylogenetic diversity was more strongly correlated to venom complexity than was overall prey species diversity, specifically implicating prey species' divergence, rather than the number of lineages alone, in the evolution of complexity. Prey phylogenetic diversity further predicted transcriptomic complexity of three of the four largest gene families in viper venom, showing that complexity evolution is a concerted response among many independent gene families. We suggest that the phylogenetic diversity of prey measures functionally relevant divergence in the targets of venom, a claim supported by sequence diversity in the coagulation cascade targets of venom. Our results support the general concept that the diversity of species in an ecological community is more important than their overall number in determining evolutionary patterns in predator trait complexity.

Sequence Divergence in Venom Genes Within and Between Montane Pitviper (Viperidae: Crotalinae: Cerrophidion) Species is Driven by Mutation-Drift Equilibrium.

Snake venom can vary both among and within species. While some groups of New World pitvipers-such as rattlesnakes-have been well studied, very little is known about the venom of montane pitvipers (Cerrophidion) found across the Mesoamerican highlands. Compared to most well-studied rattlesnakes, which are widely distributed, the isolated montane populations of Cerrophidion may facilitate unique evolutionary trajectories and venom differentiation. Here, we describe the venom gland transcriptomes for populations of C. petlalcalensis, C. tzotzilorum, and C. godmani from Mexico, and a single individual of C. sasai from Costa Rica. We explore gene expression variation in Cerrophidion and sequence evolution of toxins within C. godmani specifically. Cerrophidion venom gland transcriptomes are composed primarily of snake venom metalloproteinases, phospholipase A[Formula: see text]s (PLA[Formula: see text]s), and snake venom serine proteases. Cerrophidion petlalcalensis shows little intraspecific variation; however, C. godmani and C. tzotzilorum differ significantly between geographically isolated populations. Interestingly, intraspecific variation was mostly attributed to expression variation as we did not detect signals of selection within C. godmani toxins. Additionally, we found PLA[Formula: see text]-like myotoxins in all species except C. petlalcalensis, and crotoxin-like PLA[Formula: see text]s in the southern population of C. godmani. Our results demonstrate significant intraspecific venom variation within C. godmani and C. tzotzilorum. The toxins of C. godmani show little evidence of directional selection where variation in toxin sequence is consistent with evolution under a model of mutation-drift equilibrium. Cerrophidion godmani individuals from the southern population may exhibit neurotoxic venom activity given the presence of crotoxin-like PLA[Formula: see text]s; however, further research is required to confirm this hypothesis.

RESUMEN: El veneno de las serpientes puede variar entre y dentro de las especies. Mientras algunos grupos de viperidos del Nuevo Mundo­como las cascabeles­han sido bien estudiadas, muy poco se sabe acerca del veneno de las nauyacas de frío (Cerrophidion) que se encuentran en las zonas altas de Mesoamérica. Comparadas con las extensamente estudiadas cascabeles, que estan ampliamente distribuidas, las poblaciones de Cerrophidion, aisladas en montañas, pueden poseer trayectorias evolutivas y diferenciación en su veneno unicos. En el presente trabajo, describimos el transcriptoma de las glándulas de veneno de poblaciones de C. petlalcalensis, C. tzotzilorum, y C. godmani de México, y un individuo de C. sasai de Costa Rica. Exploramos la variación en la expresión de toxinas en Cerrophidion y la evolución en las secuencias geneticas en C. godmani específicamente. El transcriptoma de la glándula de veneno de Cerrophidion esta compuesto principalmente de Metaloproteinasas de Veneno de Serpiente, Fosfolipasas A[Formula: see text] (PLA[Formula: see text]s), y Serin Proteasas de Veneno de Serpiente. Cerrophidion petlalcalensis presenta poca variación intraespecífica; sin embargo, los transcriptomas de la glandula de veneno de C. godmani y C. tzotzilorum difieren significativamente entre poblaciones geográficamente aisladas. Curiosamente, la variación intraespecífica estuvo atribuida principalmente a la expresión de las toxinas ya que no encontramos señales de selección en las toxinas de C. godmani. Adicionalmente, encontramos miotoxinas similares a PLA[Formula: see text] en todas las especies excepto C. petlalcalensis, y PLA[Formula: see text]s similares a crotoxina en la población sureña de C. godmani. Nuestros resultados demuestran la presencia de variacion intraespecífica presente en el veneno de C. godmani y C. tzotzilorum. Las toxinas de Cerrophidion godmani muestran poca evidencia de selección direccional, y la variación en la secuencias de las toxinas es consistente con evolucion bajo un modelo de equilibrio de mutación-deriva. Algunos individuos de C. godmani de la población del sur potencialmente tienen un veneno neurotóxico dada la presencia de PLA[Formula: see text]s similares a la crotoxina, sin embargo, se necesita más evidencia para corroborar esta hipótesis.

Cryptic genetic diversity, population structure, and gene flow in the Mojave rattlesnake (Crotalus scutulatus).

The Mojave rattlesnake (Crotalus scutulatus) inhabits deserts and arid grasslands of the western United States and Mexico. Despite considerable interest in its highly toxic venom and the recognition of two subspecies, no molecular studies have characterized range-wide genetic diversity and population structure or tested species limits within C. scutulatus. We used mitochondrial DNA and thousands of nuclear loci from double-digest restriction site associated DNA sequencing to infer population genetic structure throughout the range of C. scutulatus, and to evaluate divergence times and gene flow between populations. We find strong support for several divergent mitochondrial and nuclear clades of C. scutulatus, including splits coincident with two major phylogeographic barriers: the Continental Divide and the elevational increase associated with the Central Mexican Plateau. We apply Bayesian clustering, phylogenetic inference, and coalescent-based species delimitation to our nuclear genetic data to test hypotheses of population structure. We also performed demographic analyses to test hypotheses relating to population divergence and gene flow. Collectively, our results support the existence of four distinct lineages within C. scutulatus, and genetically defined populations do not correspond with currently recognized subspecies ranges. Finally, we use approximate Bayesian computation to test hypotheses of divergence among multiple rattlesnake species groups distributed across the Continental Divide, and find evidence for co-divergence at this boundary during the mid-Pleistocene.

[Temporal diversity dynamics of the arbuscular mycorrhizal fungi of Larrea tridentata (Sesse & Mocino ex DC) Coville in a semi-arid ecosystem]. / Dinámica de la diversidad temporal de los hongos micorrícicos arbusculares de Larrea tridentata (Sesse & Mocino ex DC) Coville en un ecosistema semiárido.

Arbuscular mycorrhizal fungi (AMF) of arid and semiarid ecosystems are important for the development of plants that grow under biotic stress in wild or in agro-ecosystems. There is little information on the temporal diversity of these organisms in perennial plants from arid ecosystems in northern Mexico. On this study, the mycorrhizal colonization and the temporal diversity of AMF in the rhizosphere of Larrea tridentata, perennial plant abundant in the Chihuahuan Desert region were explored. Samples of the rhizosphere and roots of fifteen plants in each of the three sampling dates during the 2015 year were obtained. A total of 17 species of HMA belonging to 12 genera and 7 families within the phylum Glomeromycota in all three sampling dates were found. Funneliformis geosporum was the dominant species belonging to the family Glomeraceae which possess the highest genera number on L. tridentata. The highest mycorrhization percentage was in February with 83.22, followed by September and May with 75.27 and 65.27%, respectively. A maximum of 16 AM fungal species were isolated and identified from L. tridentata rhizosphere in February, 15 species in May and 12 species in September. Statistical analysis showed significant differences between sampling dates in the spores number.

Species delimitation of the blue-spotted spiny lizard within a multilocus, multispecies coalescent framework, results in the recognition of a new Sceloporus species.

Species delimitation is a major topic in systematics. Species delimitation methods based on molecular data have become more common since this approach provides insights about species identification via levels of gene flow, the degree of hybridization and phylogenetic relationships. Also, combining multilocus mitochondrial and nuclear DNA leads to more reliable conclusions about species limits. Coalescent-based species delimitation methods explicitly reveal separately evolving lineages using probabilistic approaches and testing the delimitation hypotheses for several species. Within a multispecies, multilocus, coalescent framework, we were able to clarify taxonomic uncertainties within S. cyanostictus, an endangered lizard that inhabits a narrow strip of the Chihuahuan Desert in Mexico. We included, for the first time in a phylogenetic analysis, lizards from the three populations of S. cyanostictus recognized so far (East Coahuila, West Coahuila and Nuevo León). Phylogenetic analysis corroborates the hypothesis of two separately evolving lineages, i.e. the East and West Coahuila populations, as proposed in a previous study. We also found a distant phylogenetic relationship between the lizards from Nuevo León and those of East and West Coahuila. Finally, based on the species delimitation results, we propose and describe a new species of Sceloporus: S. gadsdeni sp. nov.

Where the "ruber" Meets the Road: Using the Genome of the Red Diamond Rattlesnake to Unravel the Evolutionary Processes Driving Venom Evolution.

Understanding the proximate and ultimate causes of phenotypic variation is fundamental in evolutionary research, as such variation provides the substrate for selection to act upon. Although trait variation can arise due to selection, the importance of neutral processes is sometimes understudied. We presented the first reference-quality genome of the Red Diamond Rattlesnake (Crotalus ruber) and used range-wide 'omic data to estimate the degree to which neutral and adaptive evolutionary processes shaped venom evolution. We characterized population structure and found substantial genetic differentiation across two populations, each with distinct demographic histories. We identified significant differentiation in venom expression across age classes with substantially reduced but discernible differentiation across populations. We then used conditional redundancy analysis to test whether venom expression variation was best predicted by neutral divergence patterns or geographically variable (a)biotic factors. Snake size was the most significant predictor of venom variation, with environment, prey availability, and neutral sequence variation also identified as significant factors, though to a lesser degree. By directly including neutrality in the model, our results confidently highlight the predominant, yet not singular, role of life history in shaping venom evolution.

Ontogenetic change in the venom composition of one Mexican black-tailed rattlesnake (Crotalus molossus nigrescens) from Durango, Mexico.

To corroborate the ontogenetic shift in the venom composition of the Mexican Black-tailed Rattlesnake (Crotalus molossus nigrescens) previously reported through the census approach, we evaluated the shift in the protein profile, lethality, and proteolytic and phospholipase activities of four venom samples obtained in 2015, 2018, 2019, and 2021 from one C. m. nigrescens individual (CMN06) collected in Durango, Mexico. We demonstrated that the venom of C. m. nigrescens changed from a myotoxin-rich venom to a phospholipase A2 and snake venom metalloproteinase-rich venom. Additionally, the proteolytic and phospholipase activities increased with age, but the lethality decreased approximately three times.

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 A 2 s (PLA 2 s), 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 LD 50 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 present.

Ontogenetic Change in the Venom of Mexican Black-Tailed Rattlesnakes (Crotalus molossus nigrescens).

Ontogenetic changes in venom composition have important ecological implications due the relevance of venom in prey acquisition and defense. Additionally, intraspecific venom variation has direct medical consequences for the treatment of snakebite. However, ontogenetic changes are not well documented in most species. The Mexican Black-tailed Rattlesnake (Crotalus molossus nigrescens) is large-bodied and broadly distributed in Mexico. To document venom variation and test for ontogenetic changes in venom composition, we obtained venom samples from twenty-seven C. m. nigrescens with different total body lengths (TBL) from eight states in Mexico. The primary components in the venom were detected by reverse-phase HPLC, western blot, and mass spectrometry. In addition, we evaluated the biochemical (proteolytic, coagulant and fibrinogenolytic activities) and biological (LD50 and hemorrhagic activity) activities of the venoms. Finally, we tested for recognition and neutralization of Mexican antivenoms against venoms of juvenile and adult snakes. We detected clear ontogenetic venom variation in C. m. nigrescens. Venoms from younger snakes contained more crotamine-like myotoxins and snake venom serine proteinases than venoms from older snakes; however, an increase of snake venom metalloproteinases was detected in venoms of larger snakes. Venoms from juvenile snakes were, in general, more toxic and procoagulant than venoms from adults; however, adult venoms were more proteolytic. Most of the venoms analyzed were hemorrhagic. Importantly, Mexican antivenoms had difficulties recognizing low molecular mass proteins (<12 kDa) of venoms from both juvenile and adult snakes. The antivenoms did not neutralize the crotamine effect caused by the venom of juveniles. Thus, we suggest that Mexican antivenoms would have difficulty neutralizing some human envenomations and, therefore, it may be necessary improve the immunization mixture in Mexican antivenoms to account for low molecular mass proteins, like myotoxins.

Evidence for divergent patterns of local selection driving venom variation in Mojave Rattlesnakes (Crotalus scutulatus).

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.

Dinámica de la diversidad temporal de los hongos micorrícicos arbusculares de Larrea tridentata (Sesse & Mocino ex DC) Coville en un ecosistema semiárido / Temporal diversity dynamics of the arbuscular mycorrhizal fungi of Larrea tridentata (Sesse & Mocino ex DC) Coville in a semi-arid ecosystem

Los hongos micorrícicos arbusculares (HMA) de los ambientes áridos y semiáridos son importantes para el desarrollo de las plantas que crecen bajo estrés biótico y abiótico en áreas naturales o en agroecosistemas. Existe poca información sobre la diversidad temporal de los HMA en plantas perennes de ecosistemas áridos en el norte de México. En este estudiose evaluaron la colonización micorrícica y la diversidad temporal de los HMA en la rizosfera de Larrea tridentata, planta perenne de amplia distribución en el Desierto Chihuahuense. Se obtuvieron muestras de la rizosfera y de raíces de 15 plantas, en 3 fechas de muestreo del año 2015. Se encontró un total de 17 especies de HMA, distribuidas en 12 géneros y 7 familias dentro del phylum Glomeromycota. La especie más abundante fue Funneliformis geosporum. Esta especie pertenece a la familia Glomeraceae, la única que estuvo representada con varios géneros en L. tridentata. El mayor porcentaje de micorrización se presentó en febrero, con un 83,22%, en septiembre fue del 75,27% y en mayo del 65,27%. El muestreo realizado en febrero presentó el mayor número de especies (16), seguido por el de mayo (15) y el de septiembre (12). El análisis estadístico mostró diferencias significativas en el número de esporas entre los muestreos de febrero y los de mayo y septiembre.

Arbuscular mycorrhizal fungi (AMF) of arid and semiarid ecosystems are important for the development of plants that grow under biotic stress in wild or in agro-ecosystems. There is little information on the temporal diversity of these organisms in perennial plants from arid ecosystems in northern Mexico. On this study, the mycorrhizal colonization and the temporal diversity of AMF in the rhizosphere of Larrea tridentata, perennial plant abundant in the Chihuahuan Desert region were explored. Samples of the rhizosphere and roots of fifteen plants in each of the three sampling dates during the 2015 year were obtained. A total of 17 species of HMA belonging to 12 genera and 7 families within the phylum Glomeromycota in all three sampling dates were found. Funneliformis geosporum was the dominant species belonging to the family Glomeraceae which possess the highest genera number on L. tridentata.The highest mycorrhization percentage was in February with 83.22, followed by September and May with 75.27 and 65.27%, respectively. A maximum of 16 AM fungal species were isolated and identified from L. tridentata rhizosphere in February, 15 species in May and 12 species in September. Statistical analysis showed significant differences between sampling dates in the spores number.