De Novo Genome Assembly Highlights the Role of Lineage-Specific Gene Duplications in the Evolution of Venom in Fea's Viper (Azemiops feae).

Despite the medical significance to humans and important ecological roles filled by vipers, few high-quality genomic resources exist for these snakes outside of a few genera of pitvipers. Here we sequence, assemble, and annotate the genome of Fea's Viper (Azemiops feae). This taxon is distributed in East Asia and belongs to a monotypic subfamily, sister to the pitvipers. The newly sequenced genome resulted in a 1.56 Gb assembly, a contig N50 of 1.59 Mb, with 97.6% of the genome assembly in contigs >50 Kb, and a BUSCO completeness of 92.4%. We found that A. feae venom is primarily composed of phospholipase A2 (PLA2) proteins expressed by genes that likely arose from lineage-specific PLA2 gene duplications. Additionally, we show that renin, an enzyme associated with blood pressure regulation in mammals and known from the venoms of two viper species including A. feae, is expressed in the venom gland at comparative levels to known toxins and is present in the venom proteome. The cooption of this gene as a toxin may be more widespread in viperids than currently known. To investigate the historical population demographics of A. feae, we performed coalescent-based analyses and determined that the effective population size has remained stable over the last 100 kyr. This suggests Quaternary glacial cycles likely had minimal influence on the demographic history of A. feae. This newly assembled genome will be an important resource for studying the genomic basis of phenotypic evolution and understanding the diversification of venom toxin gene families.

An integrative view of the toxic potential of Conophis lineatus (Dipsadidae: Xenodontinae), a medically relevant rear-fanged snake.

Most traditional research on snake venoms has focused on front-fanged snake families (Viperidae, Elapidae, and Atractaspididae). However, venom is now generally accepted as being a much more broadly possessed trait within snakes, including species traditionally considered harmless. Unfortunately, due to historical inertia and methodological challenges, the toxin repertoires of non-front-fanged snake families (e.g., Colubridae, Dipsadidae, and Natricidae) have been heavily neglected despite the knowledge of numerous species capable of inflicting medically relevant envenomations. Integrating proteomic data for validation, we perform a de novo assembly and analysis of the Duvernoy's venom gland transcriptome of the Central American Road Guarder (Dipsadidae: Xenodontinae: Conophis lineatus), a species known for its potent bite. We identified 28 putative toxin transcripts from 13 toxin families in the Duvernoy's venom gland transcriptome, comprising 63.7% of total transcriptome expression. In addition to ubiquitous snake toxin families, we proteomically confirmed several atypical venom components. The most highly expressed toxins (55.6% of total toxin expression) were recently described snake venom matrix metalloproteases (svMMPs), with 48.0% of svMMP expression contributable to a novel svMMP isoform. We investigate the evolution of the new svMMP isoform in the context of rear-fanged snakes using phylogenetics. Finally, we examine the morphology of the venom apparatus using µCT and explore how the venom relates to autecology and the highly hemorrhagic effects seen in human envenomations. Importantly, we provide the most complete venom characterization of this medically relevant snake species to date, producing insights into the effects and evolution of its venom, and point to future research directions to better understand the venoms of 'harmless' non-front-fanged snakes.

An integrative view of the toxic potential of Conophis lineatus (Dipsadidae: Xenodontinae), a medically relevant rear-fanged snake

Most traditional research on snake venoms has focused on front-fanged snake families (Viperidae, Elapidae, and Atractaspididae). However, venom is now generally accepted as being a much more broadly possessed trait within snakes, including species traditionally considered harmless. Unfortunately, due to historical inertia and methodological challenges, the toxin repertoires of non-front-fanged snake families (e.g., Colubridae, Dipsadidae, and Natricidae) have been heavily neglected despite the knowledge of numerous species capable of inflicting medically relevant envenomations. Integrating proteomic data for validation, we perform a de novo assembly and analysis of the Duvernoy's venom gland transcriptome of the Central American Road Guarder (Dipsadidae: Xenodontinae: Conophis lineatus), a species known for its potent bite. We identified 28 putative toxin transcripts from 13 toxin families in the Duvernoy's venom gland transcriptome, comprising 63.7% of total transcriptome expression. In addition to ubiquitous snake toxin families, we proteomically confirmed several atypical venom components. The most highly expressed toxins (55.6% of total toxin expression) were recently described snake venom matrix metalloproteases (svMMPs), with 48.0% of svMMP expression contributable to a novel svMMP isoform. We investigate the evolution of the new svMMP isoform in the context of rear-fanged snakes using phylogenetics. Finally, we examine the morphology of the venom apparatus using μCT and explore how the venom relates to autecology and the highly hemorrhagic effects seen in human envenomations. Importantly, we provide the most complete venom characterization of this medically relevant snake species to date, producing insights into the effects and evolution of its venom, and point to future research directions to better understand the venoms of ‘harmless’ non-front-fanged snakes.

The Tiger Rattlesnake genome reveals a complex genotype underlying a simple venom phenotype.

Variation in gene regulation is ubiquitous, yet identifying the mechanisms producing such variation, especially for complex traits, is challenging. Snake venoms provide a model system for studying the phenotypic impacts of regulatory variation in complex traits because of their genetic tractability. Here, we sequence the genome of the Tiger Rattlesnake, which possesses the simplest and most toxic venom of any rattlesnake species, to determine whether the simple venom phenotype is the result of a simple genotype through gene loss or a complex genotype mediated through regulatory mechanisms. We generate the most contiguous snake-genome assembly to date and use this genome to show that gene loss, chromatin accessibility, and methylation levels all contribute to the production of the simplest, most toxic rattlesnake venom. We provide the most complete characterization of the venom gene-regulatory network to date and identify key mechanisms mediating phenotypic variation across a polygenic regulatory network.

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.

Replacement and Parallel Simplification of Nonhomologous Proteinases Maintain Venom Phenotypes in Rear-Fanged Snakes.

Novel phenotypes are commonly associated with gene duplications and neofunctionalization, less documented are the cases of phenotypic maintenance through the recruitment of novel genes. Proteolysis is the primary toxic character of many snake venoms, and ADAM metalloproteinases, named snake venom metalloproteinases (SVMPs), are largely recognized as the major effectors of this phenotype. However, by investigating original transcriptomes from 58 species of advanced snakes (Caenophidia) across their phylogeny, we discovered that a different enzyme, matrix metalloproteinase (MMP), is actually the dominant venom component in three tribes (Tachymenini, Xenodontini, and Conophiini) of rear-fanged snakes (Dipsadidae). Proteomic and functional analyses of these venoms further indicate that MMPs are likely playing an "SVMP-like" function in the proteolytic phenotype. A detailed look into the venom-specific sequences revealed a new highly expressed MMP subtype, named snake venom MMP (svMMP), which originated independently on at least three occasions from an endogenous MMP-9. We further show that by losing ancillary noncatalytic domains present in its ancestors, svMMPs followed an evolutionary path toward a simplified structure during their expansion in the genomes, thus paralleling what has been proposed for the evolution of their Viperidae counterparts, the SVMPs. Moreover, we inferred an inverse relationship between the expression of svMMPs and SVMPs along the evolutionary history of Xenodontinae, pointing out that one type of enzyme may be substituting for the other, whereas the general (metallo)proteolytic phenotype is maintained. These results provide rare evidence on how relevant phenotypic traits can be optimized via natural selection on nonhomologous genes, yielding alternate biochemical components.

Replacement and parallel simplification of non-homologous proteinases maintain venom phenotypes in rear-fanged snakes

Novel phenotypes are commonly associated with gene duplications and neofunctionalization, less documented are the cases of phenotypic maintenance through the recruitment of novel genes. Proteolysis is the primary toxic character of many snake venoms, and ADAM metalloproteinases, named Snake Venom Metalloproteinases (SVMPs), are largely recognized as the major effectors of this phenotype. However, by investigating original transcriptomes from 58 species of advanced snakes (Caenophidia) across their phylogeny, we discovered that a different enzyme, matrix metalloproteinase (MMP), is actually the dominant venom component in three tribes (Tachymenini, Xenodontini, and Conophiini) of rear-fanged snakes (Dipsadidae). Proteomic and functional analyses of these venoms further indicate that MMPs are likely playing an ‘SVMP-like’ function in the proteolytic phenotype. A detailed look into the venom-specific sequences revealed a new highly expressed MMP subtype, named snake venom MMP (svMMP), which originated independently on at least three occasions from an endogenous MMP-9. We further show that by losing ancillary non-catalytic domains present in its ancestors, svMMPs followed an evolutionary path toward a simplified structure during their expansion in the genomes, thus paralleling what has been proposed for the evolution of their Viperidae counterparts, the SVMPs. Moreover, we inferred an inverse relationship between the expression of svMMPs and SVMPs along the evolutionary history of Xenodontinae, pointing out that one type of enzyme may be substituting for the other, while the general (metallo)proteolytic phenotype is maintained. These results provide rare evidence on how relevant phenotypic traits can be optimized via natural selection on non-homologous genes, yielding alternate biochemical components.