Topology Testing and Demographic Modeling Illuminate a Novel Speciation Pathway in the Greater Caribbean Sea Following the Formation of the Isthmus of Panama.

Recent genomic analyses have highlighted the prevalence of speciation with gene flow in many taxa and have underscored the importance of accounting for these reticulate evolutionary processes when constructing species trees and generating parameter estimates. This is especially important for deepening our understanding of speciation in the sea where fast moving ocean currents, expanses of deep water, and periodic episodes of sea level rise and fall act as soft and temporary allopatric barriers that facilitate both divergence and secondary contact. Under these conditions, gene flow is not expected to cease completely while contemporary distributions are expected to differ from historical ones. Here we conduct range-wide sampling for Pederson's cleaner shrimp (Ancylomenes pedersoni), a species complex from the Greater Caribbean that contains three clearly delimited mitochondrial lineages with both allopatric and sympatric distributions. Using mtDNA barcodes and a genomic ddRADseq approach, we combine classic phylogenetic analyses with extensive topology testing and demographic modeling (10 site frequency replicates x 45 evolutionary models x 50 model simulations/replicate = 22,500 simulations) to test species boundaries and reconstruct the evolutionary history of what was expected to be a simple case study. Instead, our results indicate a history of allopatric divergence, secondary contact, introgression, and endemic hybrid speciation that we hypothesize was driven by the final closure of the Isthmus of Panama and the strengthening of the Gulf Stream Current ~3.5 million years ago. The history of this species complex recovered by model-based methods that allow reticulation differs from that recovered by standard phylogenetic analyses and is unexpected given contemporary distributions. The geologically and biologically meaningful insights gained by our model selection analyses illuminate what is likely a novel pathway of species formation not previously documented that resulted from one of the most biogeographically significant events in Earth's history.

Individual Variability in Bothropsatrox Snakes Collected from Different Habitats in the Brazilian Amazon: New Findings on Venom Composition and Functionality.

Differences in snake venom composition occur across all taxonomic levels and it has been argued that this variation represents an adaptation that has evolved to facilitate the capture and digestion of prey and evasion of predators. Bothrops atrox is a terrestrial pitviper that is distributed across the Amazon region, where it occupies different habitats. Using statistical analyses and functional assays that incorporate individual variation, we analyzed the individual venom variability in B. atrox snakes from four different habitats (forest, pasture, degraded area, and floodplain) in and around the Amazon River in Brazil. We observed venom differentiation between spatially distinct B. atrox individuals from the different habitats, with venom variation due to both common (high abundance) and rare (low abundance) proteins. Moreover, differences in the composition of the venoms resulted in individual variability in functionality and heterogeneity in the lethality to mammals and birds, particularly among the floodplain snakes. Taken together, the data obtained from individual venoms of B. atrox snakes, captured in different habitats from the Brazilian Amazon, support the hypothesis that the differential distribution of protein isoforms results in functional distinctiveness and the ability of snakes with different venoms to have variable toxic effects on different prey.

Tracking the recruitment and evolution of snake toxins using the evolutionary context provided by the Bothrops jararaca genome.

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option.

Phylogenetically diverse diets favor more complex venoms in North American pitvipers

The role of natural selection in the evolution of trait complex-ity can be characterized by testing hypothesized links betweencomplex forms and their functions across species. Predatory ven-oms are composed of multiple proteins that collectively function toincapacitate prey. Venom complexity fluctuates over evolutionarytimescales, with apparent increases and decreases in complexity,and yet the causes of this variation are unclear. We tested alterna-tive hypotheses linking venom complexity and ecological sourcesof selection from diet in the largest clade of front-fanged ven-omous snakes in North America: the rattlesnakes, copperheads,cantils, and cottonmouths. We generated independent transcrip-tomic and proteomic measures of venom complexity and collatedseveral natural history studies to quantify dietary variation. Wethen constructed genome-scale phylogenies for these snakes forcomparative analyses. Strikingly, prey phylogenetic diversity wasmore strongly correlated to venom complexity than was overallprey species diversity, specifically implicating prey species’ diver-gence, rather than the number of lineages alone, in the evolutionof complexity. Prey phylogenetic diversity further predicted tran-scriptomic complexity of three of the four largest gene familiesin viper venom, showing that complexity evolution is a concertedresponse among many independent gene families. We suggest thatthe phylogenetic diversity of prey measures functionally relevantdivergence in the targets of venom, a claim supported by sequencediversity in the coagulation cascade targets of venom. Our resultssupport the general concept that the diversity of species in an eco-logical community is more important than their overall number indetermining evolutionary patterns in predator trait complexity.

Individual variability in bothrops atrox snakes collected from different habitats in the Brazilian Amazon: new findings on venom composition and functionality

Differences in snake venom composition occur across all taxonomic levels and it has been argued that this variation represents an adaptation that has evolved to facilitate the capture and digestion of prey and evasion of predators. Bothrops atrox is a terrestrial pitviper that is distributed across the Amazon region, where it occupies different habitats. Using statistical analyses and functional assays that incorporate individual variation, we analyzed the individual venom variability in B. atrox snakes from four different habitats (forest, pasture, degraded area, and floodplain) in and around the Amazon River in Brazil. We observed venom differentiation between spatially distinct B. atrox individuals from the different habitats, with venom variation due to both common (high abundance) and rare (low abundance) proteins. Moreover, differences in the composition of the venoms resulted in individual variability in functionality and heterogeneity in the lethality to mammals and birds, particularly among the floodplain snakes. Taken together, the data obtained from individual venoms of B. atrox snakes, captured in different habitats from the Brazilian Amazon, support the hypothesis that the differential distribution of protein isoforms results in functional distinctiveness and the ability of snakes with different venoms to have variable toxic effects on different prey.

Tracking the recruitment and evolution of snake toxins using the evolutionary context provided by the Bothrops jararaca genome

Venom is a key adaptive innovation in snakes, and how nonvenom genes were co-opted to become part of the toxin arsenal is a significant evolutionary question. While this process has been investigated through the phylogenetic reconstruction of toxin sequences, evidence provided by the genomic context of toxin genes remains less explored. To investigate the process of toxin recruitment, we sequenced the genome of Bothrops jararaca, a clinically relevant pitviper. In addition to producing a road map with canonical structures of genes encoding 12 toxin families, we inferred most of the ancestral genes for their loci. We found evidence that 1) snake venom metalloproteinases (SVMPs) and phospholipases A2 (PLA2) have expanded in genomic proximity to their nonvenomous ancestors; 2) serine proteinases arose by co-opting a local gene that also gave rise to lizard gilatoxins and then expanded; 3) the bradykinin-potentiating peptides originated from a C-type natriuretic peptide gene backbone; and 4) VEGF-F was co-opted from a PGF-like gene and not from VEGF-A. We evaluated two scenarios for the original recruitment of nontoxin genes for snake venom: 1) in locus ancestral gene duplication and 2) in locus ancestral gene direct co-option. The first explains the origins of two important toxins (SVMP and PLA2), while the second explains the emergence of a greater number of venom components. Overall, our results support the idea of a locally assembled venom arsenal in which the most clinically relevant toxin families expanded through posterior gene duplications, regardless of whether they originated by duplication or gene co-option

Venom variation in Bothrops asper lineages from North-Western South America.

Bothrops asper is a venomous pitviper that is widely distributed and of clinical importance in Mesoamerica and northern South America, where it is responsible for 50-80% of all envenomations by Viperidae species. Previous work suggests that B. asper has a complex phylogeographic structure, with the existence of multiple evolutionarily distinct lineages, particularly in the inter-Andean valleys of north South America. To explore the impact of the evolutionary history of B. asper on venom composition, we have investigated geographic variation in the venom proteome of this species from the populations from the Pacific side of Ecuador and south-western Colombia. Among the 21 classes of venom components identified, proteins from mainly four major toxin families, snake venom metalloproteases (PI- and PII-SVMP), phospholipases A2 (K49- and D49-PLA2s), serine proteinases (SVSP), and C-type lectins-like (CTL) proteins are major contributors to the geographic variability in venom. Principal component analyses demonstrate significant differences in venom composition between B. asper lineages previously identified through combination of molecular, morphological and geographical data, and provide additional insights into the selection pressures modulating venom phenotypes on a geographic scale. In particular, altitudinal zonation within the Andean mountain range stands out as a key ecological factor promoting diversification in venom. In addition, the pattern of distribution of PLA2 molecules among B. asper venoms complements phylogenetic analysis in the reconstruction of the dispersal events that account for the current biogeographic distribution of the present-day species' phylogroups. Ontogenic variation was also evident among venoms from some Ecuadorian lineages, although this age-related variation was less extreme than reported in B. asper venoms from Costa Rica. The results of our study demonstrate a significant impact of phylogenetic history on venom composition in a pitviper and show how analyses of this variation can illuminate the timing of the cladogenesis and ecological events that shaped the current distribution of B. asper lineages. BIOLOGICAL SIGNIFICANCE: Bothrops asper, called "the ultimate pitviper" due to its defensive behavior, large body size, and medical importance, represents a species complex that is widely distributed from southern México southwards across north-western South America to north-western Perú. This work reports the characterization of the venom proteomes of B. asper lineages from the Pacific sides of Ecuador and south-western Colombia. Multivariate analyses indicate that variability in venom composition among the B. asper lineages is driven by proteins from four major toxin families, presumably in response to selection pressures created by recent and historical ecological conditions created by geological and climatic events from the Pliocene-Pleistocene to the present along the Central and South American Continental Divide. The emerging biogeographic pattern of venom variation, interpreted in the context of the current phylogenetic hypotheses, support and complement previously proposed evolutionary Plio-Pleistocene dispersal events that shaped the present-day distribution range of B. asper lineages. In addition, our venomics data indicate the occurrence of genetic exchange between Colombian and Pacific Costa Rican populations, which may have occurred during the second wave of B. asper migration into Mesoamerica. Our work represents a foundation for a future broader sampling and more complete "-omics" analyses to deepen our understanding of the patterns and causes of venom variation in this medically important pitviper.

Molecular mechanisms underlying intraspecific variation in snake venom.

Elucidating the molecular mechanisms underlying snake venom variability provides important clues for understanding how the biological functions of this powerful toxic arsenal evolve. We analyzed in detail individual transcripts and venom protein isoforms produced by five specimens of a venomous snake (Bothrops atrox) from two nearby but genetically distinct populations from the Brazilian Amazon rainforest which show functional similarities in venom properties. Individual variation was observed among the venoms of these specimens, but the overall abundance of each general toxin family was conserved both in transcript and in venom protein levels. However, when expression of independent paralogues was analyzed, remarkable differences were observed within and among each toxin group, both between individuals and between populations. Transcripts for functionally essential venom proteins ("core function" proteins) were highly expressed in all specimens and showed similar transcription/translation rates. In contrast, other paralogues ("adaptive" proteins) showed lower expression levels and the toxins they coded for varied among different individuals. These results provide support for the inferences that (a) expression and translational differences play a greater role in defining adaptive variation in venom phenotypes than does sequence variation in protein coding genes and (b) convergent adaptive venom phenotypes can be generated through different molecular mechanisms. SIGNIFICANCE: Analysis of individual transcripts and venom protein isoforms produced by specimens of a venomous snake (Bothrops atrox), from the Brazilian Amazon rainforest, revealed that transcriptional and translational mechanisms contribute to venom phenotypic variation. Our finding of evidence for high expression of toxin proteins with conserved function supports the hypothesis that the venom phenotype consists of two kinds of proteins: conserved "core function" proteins that provide essential functional activities with broader relevance and less conserved "adaptive" proteins that vary in expression and may permit customization of protein function. These observations allowed us to suggest that genetic mechanisms controlling venom variability are not restricted to selection of gene copies or mutations in structural genes but also to selection of the mechanisms controlling gene expression, contributing to the plasticity of this important phenotype for venomous snakes.

Molecular mechanisms underlying intraspecific variation in snake venom

Elucidating the molecular mechanisms underlying snake venom variability provides important clues for understanding how the biological functions of this powerful toxic arsenal evolve. We analyzed in detail individual transcripts and venom protein isoforms produced by five specimens of a venomous snake (Bothrops atrox) from two nearby but genetically distinct populations from the Brazilian Amazon rainforest which show functional similarities in venom properties. Individual variation was observed among the venoms of these specimens, but the overall abundance of each general toxin family was conserved both in transcript and in venom protein levels. However, when expression of independent paralogues was analyzed, remarkable differences were observed within and among each toxin group, both between individuals and between populations. Transcripts for functionally essential venom proteins ("core function" proteins) were highly expressed in all specimens and showed similar transcription/translation rates. In contrast, other paralogues ("adaptive" proteins) showed lower expression levels and the toxins they coded for varied among different individuals. These results provide support for the inferences that (a) expression and translational differences play a greater role in defining adaptive variation in venom phenotypes than does sequence variation in protein coding genes and (b) convergent adaptive venom phenotypes can be generated through different molecular mechanisms. Significance: Analysis of individual transcripts and venom protein isoforms produced by specimens of a venomous snake (Bothrops atrox), from the Brazilian Amazon rainforest, revealed that transcriptional and translational mechanisms contribute to venom phenotypic variation. Our finding of evidence for high expression of toxin proteins with conserved function supports the hypothesis that the venom phenotype consists of two kinds of proteins: conserved "core function" proteins that provide essential functional activities with broader relevance and less conserved "adaptive" proteins that vary in expression and may permit customization of protein function. These observations allowed us to suggest that genetic mechanisms controlling venom variability are not restricted to selection of gene copies or mutations in structural genes but also to selection of the mechanisms controlling gene expression, contributing to the plasticity of this important phenotype for venomous snakes.

Recent lineage diversification in a venomous snake through dispersal across the Amazon River

Identifying the evolutionary and ecological mechanisms that drive lineage diversification in the species-rich tropics is of broad interest to evolutionary biologists. Here, we use phylogeographical and demographic analyses of genome-scale RADseq data to assess the impact of a large geographical feature, the Amazon River, on lineage formation in a venomous pitviper, Bothrops atrox. We compared genetic differentiation in samples from four sites near Santarem, Brazil, that spanned the Amazon and represented major habitat types. A species delimitation analysis identified each population as a distinct evolutionary lineage while a species tree analysis with populations as taxa revealed a phylogenetic tree consistent with dispersal across the Amazon from north to south. Phylogenetic analyses of mitochondrial DNA variation confirmed this pattern and suggest that all lineages originated during the mid- to late Pleistocene. Historical demographic analyses support a population model of lineage formation through isolation between lineages with low ongoing migration between large populations and reject a model of differentiation through isolation by distance alone. The results provide a rare example of a phylogeographical pattern demonstrating dispersal over evolutionary timescales across a large tropical river and suggest a role for the Amazon River as a driver of in situ divergence both by impeding (but not preventing) gene flow and through parapatric differentiation along an ecological gradient.

Marshes as "Mountain Tops": Genetic Analyses of the Critically Endangered São Paulo Marsh Antwren (Aves: Thamnophilidae).

Small populations of endangered species can be impacted by genetic processes such as drift and inbreeding that reduce population viability. As such, conservation genetic analyses that assess population levels of genetic variation and levels of gene flow can provide important information for managing threatened species. The São Paulo Marsh Antwren (Formicivora paludicola) is a recently-described and critically endangered bird from São Paulo State (Brazil) whose total estimated population is around 250-300 individuals, distributed in only 15 isolated marshes around São Paulo metropolitan region. We used microsatellite DNA markers to estimate the population genetic characteristics of the three largest remaining populations of this species all within 60 km of each other. We detected a high and significant genetic structure between all populations (overall FST = 0.103) which is comparable to the highest levels of differentiation ever documented for birds, (e.g., endangered birds found in isolated populations on the tops of African mountains), but also evidence for first-generation immigrants, likely from small local unsampled populations. Effective population sizes were small (between 28.8-99.9 individuals) yet there are high levels of genetic variability within populations and no evidence for inbreeding. Conservation implications of this work are that the high levels of genetic structure suggests that translocations between populations need to be carefully considered in light of possible local adaptation and that remaining populations of these birds should be managed as conservation units that contain both main populations studied here but also small outlying populations which may be a source of immigrants.

Proteomic and toxicological profiling of the venom of Bothrocophias campbelli, a pitviper species from Ecuador and Colombia.

Detailed snake venom proteomes for nearly a hundred species in different pitviper genera have accumulated using 'venomics' methodologies. However, venom composition for some lineages remains poorly known. Bothrocophias (toad-headed pitvipers) is a genus restricted to the northwestern portion of South America for which information on venom composition is lacking. Here, we describe the protein composition, toxicological profiling, and antivenom neutralization of the venom of Bothrocophias campbelli, a species distributed in Colombia and Ecuador. Our analyses show that its venom mainly consists of phospholipases A2 (43.1%), serine proteinases (21.3%), and metalloproteinases (15.8%). The low proportion of metalloproteinases and high amount of a Lys49 phospholipase A2 homologue correlate well with the low hemorrhagic and high myotoxic effects found. Overall, B. campbelli venom showed a simpler composition compared to other crotalines in the region. A polyvalent antivenom prepared with a mixture of Bothrops asper, Crotalus simus, and Lachesis stenophrys venoms cross-recognized B. campbelli venom and neutralized its lethal effect in mice, albeit with a lower potency than for B. asper venom. Additional work comparing B. campbelli venom properties with those of related species could help understand the evolution of different venom protein families during the South American radiation of New World pitvipers.

Genetic structure in a tropical lek-breeding bird, the blue manakin (Chiroxiphia caudata) in the Brazilian Atlantic Forest.

Determining the genetic structure of tropical bird populations is important for assessing potential genetic effects of future habitat fragmentation and for testing hypotheses about evolutionary mechanisms promoting diversification. Here we used 10 microsatellite DNA loci to describe levels of genetic differentiation for five populations of the lek-mating blue manakin (Chiroxiphia caudata), sampled along a 414-km transect within the largest remaining continuous tract of the highly endangered Atlantic Forest habitat in southeast Brazil. We found small but significant levels of differentiation between most populations. F(ST) values varied from 0.0 to 0.023 (overall F(ST)=0.012) that conformed to a strong isolation by distance relationship, suggesting that observed levels of differentiation are a result of migration-drift equilibrium. N(e)m values estimated using a coalescent-based method were small (

Lineage origin and expansion of a Neotropical migrant songbird after recent glaciation events.

Birds of the Northern Hemisphere often harbour the genetic signature of postglaciation expansion but analyses identifying the location of refugia and the directionality of expansions are rare. Here we explore the evolutionary history of yellow warbler lineages, focusing on how these lineages recolonized their current range. We genotyped samples from 696 yellow warblers via direct sequencing of a 333-bp control region I mitochondrial DNA fragment or lineage-specific genotyping. Phylogenetic analysis revealed two monophyletic clades: a highly migratory group including previously identified eastern and western lineages and a less migratory group including a lineage consisting of tropical residents and a new 'southern' lineage localized in southwest United States. We then modelled the expansion of the eastern and western lineages, identified the location of potential refugia and assessed the importance of migration as a historical factor promoting gene flow. The expansion of the eastern lineage proceeded from a main refugia in the eastern United States, with possible contribution of an additional local refugia. In the western lineage, the expansion proceeded from a single refugia possibly located in western United States. Because two lineages overlapped to varying degrees in central North America, we suggest that the Canadian Prairies offered a bridge of riparian habitats where the lineages met after glacier retreat, while the US Central Great Plains acted as a barrier that limited secondary contact. Finally, gene flow was more important along the north-south axis of migration than away from it, suggesting spring migration played a role in the dispersal of lineages.