The benefits of contributing to the citizen science platform iNaturalist as an identifier.

As the number of observations submitted to the citizen science platform iNaturalist continues to grow, it is increasingly important that these observations can be identified to the finest taxonomic level, maximizing their value for biodiversity research. Here, we explore the benefits of acting as an identifier on iNaturalist.

Genomics reveals broad hybridization in deeply divergent Palearctic grass and water snakes (Natrix spp.).

Understanding speciation is one of the cornerstones of biological diversity research. Currently, speciation is often understood as a continuous process of divergence that continues until genetic or other incompatibilities minimize or prevent interbreeding. The Palearctic snake genus Natrix is an ideal group to study speciation, as it comprises taxa representing distinct stages of the speciation process, ranging from widely interbreeding parapatric taxa through parapatric species with very limited gene flow in narrow hybrid zones to widely sympatric species. To understand the evolution of reproductive isolation through time, we have sequenced the genomes of all five species within this genus and two additional subspecies. We used both long-read and short-read methods to sequence and de-novo-assemble two high-quality genomes (Natrix h. helvetica, Natrix n. natrix) to their 1.7 Gb length with a contig N50 of 4.6 Mbp and 1.5 Mbp, respectively, and used these as references to assemble the remaining short-read-based genomes. Our phylogenomic analyses yielded a well-supported dated phylogeny and evidence for a surprisingly complex history of interspecific gene flow, including between widely sympatric species. Furthermore, evidence for gene flow was also found for currently allopatric species pairs. Genetic exchange among these well-defined, distinct, and several million-year-old reptile species emphasizes that speciation and maintenance of species distinctness can occur despite continued genetic exchange.

Mohave Rattlesnake (Crotalus scutulatus) Identification Revisited.

Crotalus scutulatus (Mohave rattlesnake) is a clinically significant pit viper broadly distributed across much of the arid southwestern United States and mainland Mexico. Identification of C scutulatus is a concern among emergency medical service and emergency department personnel owing to its reputation for severe envenomations and difficulty in visually differentiating between C scutulatus and other species, primarily Crotalus atrox (western diamond-backed rattlesnake). We contrast distinctive characteristics of C scutulatus, C atrox, and 3 other sympatric species: Crotalus molossus, Crotalus ornatus, and Crotalus viridis (western and eastern black-tailed rattlesnakes and prairie rattlesnake, respectively). Greenish coloration eliminates C atrox but does not confirm C scutulatus. Obvious coarse and fine speckling of the dorsal pattern and a pale postocular stripe intersecting the mouth characterize C atrox. Dorsal speckling is insignificant or absent in the other species, whereas the pale postocular stripe passes above the mouth in C scutulatus and C viridis and is absent in C molossus and C ornatus. Tails boldly ringed with alternating black and white or contrasting shades of gray are shared by C atrox and C scutulatus, respectively, but a lack of boldly ringed tails characterizes the other species. The proximal rattle segment is yellow and black, or entirely yellow, in C scutulatus but black in the others. The most reliable visual identifications are based on evaluations of multiple traits, all of which are variable to some extent. Traits such as tail ring width and the size and number of crown scales have frequently been overemphasized in the past.

King or royal family? Testing for species boundaries in the King Cobra, Ophiophagus hannah (Cantor, 1836), using morphology and multilocus DNA analyses.

In widespread species, the diverse ecological conditions in which the populations occur, and the presence of many potential geographical barriers through their range are expected to have created ample opportunities for the evolution of distinct, often cryptic lineages. In this work, we tested for species boundaries in one such widespread species, the king cobra, Ophiophagus hannah (Cantor, 1836), a largely tropical elapid snake distributed across the Oriental realm. Based on extensive geographical sampling across most of the range of the species, we initially tested for candidate species (CS) using Maximum-Likelihood analysis of mitochondrial genes. We then tested the resulting CS using both morphological data and sequences of three single-copy nuclear genes. We used snapclust to determine the optimal number of clusters in the nuclear dataset, and Bayesian Phylogenetics and Phylogeography (BPP) to test for likely species status. We used non-metric multidimensional scaling (nMDS) analysis for discerning morphological separation. We recovered four independently evolving, geographically separated lineages that we consider Confirmed Candidate Species: (1) Western Ghats lineage; (2) Indo-Chinese lineage (3) Indo-Malayan lineage; (4) Luzon Island lineage, in the Philippine Archipelago. We discuss patterns of lineage divergence, particularly in the context of low morphological divergence, and the conservation implications of recognizing several endemic king cobra lineages.

Unexpected lack of specialisation in the flow properties of spitting cobra venom.

Venom spitting is a defence mechanism based on airborne venom delivery used by a number of different African and Asian elapid snake species ('spitting cobras'; Naja spp. and Hemachatus spp.). Adaptations underpinning venom spitting have been studied extensively at both behavioural and morphological level in cobras, but the role of the physical properties of venom itself in its effective projection remains largely unstudied. We hereby provide the first comparative study of the physical properties of venom in spitting and non-spitting cobras. We measured the viscosity, protein concentration and pH of the venom of 13 cobra species of the genus Naja from Africa and Asia, alongside the spitting elapid Hemachatus haemachatus and the non-spitting viper Bitis arietans By using published microCT scans, we calculated the pressure required to eject venom through the fangs of a spitting and a non-spitting cobra. Despite the differences in the modes of venom delivery, we found no significant differences between spitters and non-spitters in the rheological and physical properties of the studied venoms. Furthermore, all analysed venoms showed a Newtonian flow behaviour, in contrast to previous reports. Although our results imply that the evolution of venom spitting did not significantly affect venom viscosity, our models of fang pressure suggests that the pressure requirements to eject venom are lower in spitting cobras than in non-spitting cobras.

Unexpected lack of specialisation in the flow properties of spitting cobra venom.

Venom spitting is a defence mechanism based on airborne venom delivery used by a number of different African and Asian elapid snake species ('spitting cobras'; Naja spp. and Hemachatus spp.). Adaptations underpinning venom spitting have been studied extensively at both behavioural and morphological level in cobras, but the role of the physical properties of venom itself in its effective projection remains largely unstudied. We hereby provide the first comparative study of the physical properties of venom in spitting and non-spitting cobras. We measured the viscosity, protein concentration and pH of the venom of 13 cobra species of the genus Naja from Africa and Asia, alongside the spitting elapid Hemachatus haemachatus and the non-spitting viper Bitis arietans. By using published microCT scans, we calculated the pressure required to eject venom through the fangs of a spitting and a non-spitting cobra. Despite the differences in the modes of venom delivery, we found no significant differences between spitters and non-spitters in the rheological and physical properties of the studied venoms. Furthermore, all analysed venoms showed a Newtonian flow behaviour, in contrast to previous reports. Although our results imply that the evolution of venom spitting did not significantly affect venom viscosity, our models of fang pressure suggests that the pressure requirements to eject venom are lower in spitting cobras than in non-spitting cobras.

When one phenotype is not enough: divergent evolutionary trajectories govern venom variation in a widespread rattlesnake species.

Understanding the origin and maintenance of phenotypic variation, particularly across a continuous spatial distribution, represents a key challenge in evolutionary biology. For this, animal venoms represent ideal study systems: they are complex, variable, yet easily quantifiable molecular phenotypes with a clear function. Rattlesnakes display tremendous variation in their venom composition, mostly through strongly dichotomous venom strategies, which may even coexist within a single species. Here, through dense, widespread population-level sampling of the Mojave rattlesnake, Crotalus scutulatus, we show that genomic structural variation at multiple loci underlies extreme geographical variation in venom composition, which is maintained despite extensive gene flow. Unexpectedly, neither diet composition nor neutral population structure explain venom variation. Instead, venom divergence is strongly correlated with environmental conditions. Individual toxin genes correlate with distinct environmental factors, suggesting that different selective pressures can act on individual loci independently of their co-expression patterns or genomic proximity. Our results challenge common assumptions about diet composition as the key selective driver of snake venom evolution and emphasize how the interplay between genomic architecture and local-scale spatial heterogeneity in selective pressures may facilitate the retention of adaptive functional polymorphisms across a continuous space.

Deconstructing compassionate conservation.

Compassionate conservation focuses on 4 tenets: first, do no harm; individuals matter; inclusivity of individual animals; and peaceful coexistence between humans and animals. Recently, compassionate conservation has been promoted as an alternative to conventional conservation philosophy. We believe examples presented by compassionate conservationists are deliberately or arbitrarily chosen to focus on mammals; inherently not compassionate; and offer ineffective conservation solutions. Compassionate conservation arbitrarily focuses on charismatic species, notably large predators and megaherbivores. The philosophy is not compassionate when it leaves invasive predators in the environment to cause harm to vastly more individuals of native species or uses the fear of harm by apex predators to terrorize mesopredators. Hindering the control of exotic species (megafauna, predators) in situ will not improve the conservation condition of the majority of biodiversity. The positions taken by so-called compassionate conservationists on particular species and on conservation actions could be extended to hinder other forms of conservation, including translocations, conservation fencing, and fertility control. Animal welfare is incredibly important to conservation, but ironically compassionate conservation does not offer the best welfare outcomes to animals and is often ineffective in achieving conservation goals. Consequently, compassionate conservation may threaten public and governmental support for conservation because of the limited understanding of conservation problems by the general public.

Deconstrucción de la Conservación Compasiva Resumen La conservación compasiva se enfoca en cuatro principios: no causar daño; los individuos importan; la integración de los animales individualmente; y la coexistencia pacífica entre los humanos u los animales. Recientemente, la conservación compasiva ha sido promovida como una alternativa a la filosofía convencional de la conservación. Creemos que los ejemplos presentados por los conservacionistas compasivos han sido elegidos arbitraria o deliberadamente por estar enfocados en los mamíferos; por ser inherentes y no compasivos; y por ofrecer soluciones de conservación poco efectivas. La conservación compasiva se enfoca arbitrariamente en las especies carismáticas, principalmente los grandes depredadores y los megaherbívoros. La filosofía no es compasiva cuando deja que los depredadores invasores dentro del ambiente causen daño a un vasto número de individuos nativos o usa el miedo al daño por superdepredadores para aterrorizar a los mesodepredadores. El entorpecimiento del control de especies exóticas (megafauna, depredadores) in situ no mejorará las condiciones de conservación de la mayoría de la biodiversidad, incluso si los conservacionistas compasivos no dañan a los individuos exóticos. Las posiciones que toman los llamados conservacionistas compasivos sobre especies particulares y sobre las acciones de conservación podrían extenderse para entorpecer otros tipos de conservación, incluyendo las reubicaciones, el encercado para la conservación y el control de la fertilidad. El bienestar animal es increíblemente importante para la conservación e irónicamente, la conservación compasiva no ofrece los mejores resultados de bienestar para los animales y comúnmente es poco efectiva en el logro de los objetivos de conservación. Como consecuencia, la conservación compasiva puede poner en peligro el apoyo público y del gobierno que tiene la conservación debido al entendimiento poco limitado que tiene el público general sobre los problemas de conservación.

Widespread convergence in toxin resistance by predictable molecular evolution.

The question about whether evolution is unpredictable and stochastic or intermittently constrained along predictable pathways is the subject of a fundamental debate in biology, in which understanding convergent evolution plays a central role. At the molecular level, documented examples of convergence are rare and limited to occurring within specific taxonomic groups. Here we provide evidence of constrained convergent molecular evolution across the metazoan tree of life. We show that resistance to toxic cardiac glycosides produced by plants and bufonid toads is mediated by similar molecular changes to the sodium-potassium-pump (Na(+)/K(+)-ATPase) in insects, amphibians, reptiles, and mammals. In toad-feeding reptiles, resistance is conferred by two point mutations that have evolved convergently on four occasions, whereas evidence of a molecular reversal back to the susceptible state in varanid lizards migrating to toad-free areas suggests that toxin resistance is maladaptive in the absence of selection. Importantly, resistance in all taxa is mediated by replacements of 2 of the 12 amino acids comprising the Na(+)/K(+)-ATPase H1-H2 extracellular domain that constitutes a core part of the cardiac glycoside binding site. We provide mechanistic insight into the basis of resistance by showing that these alterations perturb the interaction between the cardiac glycoside bufalin and the Na(+)/K(+)-ATPase. Thus, similar selection pressures have resulted in convergent evolution of the same molecular solution across the breadth of the animal kingdom, demonstrating how a scarcity of possible solutions to a selective challenge can lead to highly predictable evolutionary responses.

Multi-locus phylogeny and species delimitation of Australo-Papuan blacksnakes (Pseudechis Wagler, 1830: Elapidae: Serpentes).

Genetic analyses of Australasian organisms have resulted in the identification of extensive cryptic diversity across the continent. The venomous elapid snakes are among the best-studied organismal groups in this region, but many knowledge gaps persist: for instance, despite their iconic status, the species-level diversity among Australo-Papuan blacksnakes (Pseudechis) has remained poorly understood due to the existence of a group of cryptic species within the P. australis species complex, collectively termed "pygmy mulga snakes". Using two mitochondrial and three nuclear loci we assess species boundaries within the genus using Bayesian species delimitation methods and reconstruct their phylogenetic history using multispecies coalescent approaches. Our analyses support the recognition of 10 species, including all of the currently described pygmy mulga snakes and one undescribed species from the Northern Territory of Australia. Phylogenetic relationships within the genus are broadly consistent with previous work, with the recognition of three major groups, the viviparous red-bellied black snake P. porphyriacus forming the sister species to two clades consisting of ovoviviparous species.

Medically important differences in snake venom composition are dictated by distinct postgenomic mechanisms.

Variation in venom composition is a ubiquitous phenomenon in snakes and occurs both interspecifically and intraspecifically. Venom variation can have severe outcomes for snakebite victims by rendering the specific antibodies found in antivenoms ineffective against heterologous toxins found in different venoms. The rapid evolutionary expansion of different toxin-encoding gene families in different snake lineages is widely perceived as the main cause of venom variation. However, this view is simplistic and disregards the understudied influence that processes acting on gene transcription and translation may have on the production of the venom proteome. Here, we assess the venom composition of six related viperid snakes and compare interspecific changes in the number of toxin genes, their transcription in the venom gland, and their translation into proteins secreted in venom. Our results reveal that multiple levels of regulation are responsible for generating variation in venom composition between related snake species. We demonstrate that differential levels of toxin transcription, translation, and their posttranslational modification have a substantial impact upon the resulting venom protein mixture. Notably, these processes act to varying extents on different toxin paralogs found in different snakes and are therefore likely to be as important as ancestral gene duplication events for generating compositionally distinct venom proteomes. Our results suggest that these processes may also contribute to altering the toxicity of snake venoms, and we demonstrate how this variability can undermine the treatment of a neglected tropical disease, snakebite.

Phylogeny and diversification of mountain vipers (Montivipera, Nilson et al., 2001) triggered by multiple Plio-Pleistocene refugia and high-mountain topography in the Near and Middle East.

The Near and Middle East is a hotspot of biodiversity, but the region remains underexplored at the level of genetic biodiversity. Here, we present an extensive molecular phylogeny of the viperid snake genus Montivipera, including all known taxa. Based on nuclear and mitochondrial data, we present novel insights into the phylogeny of the genus and review the status of its constituent species. Maximum likelihood methods revealed a montane origin of Montivipera at 12.3Mya. We then analyzed factors of mountain viper diversity. Our data support substantial changes in effective population size through Plio-Pleistocene periods. We conclude that climatic oscillations were drivers of allopatric speciation, and that mountain systems of the Near and Middle East have strongly influenced the evolution and survival of taxa, because climatic and topographical heterogeneities induced by mountains have played a crucial role as filters for dispersal and as multiple refugia. The wide diversity of montane microhabitats enabled mountain vipers to retain their ecological niche during climatic pessima. In consequence the varied geological and topographical conditions between refugia favoured genetic isolation and created patterns of species richness resulting in the formation of neoendemic taxa. Our data support high concordance between geographic distributions of Montivipera haplotypes with putative plant refugia.

The king cobra genome reveals dynamic gene evolution and adaptation in the snake venom system.

Snakes are limbless predators, and many species use venom to help overpower relatively large, agile prey. Snake venoms are complex protein mixtures encoded by several multilocus gene families that function synergistically to cause incapacitation. To examine venom evolution, we sequenced and interrogated the genome of a venomous snake, the king cobra (Ophiophagus hannah), and compared it, together with our unique transcriptome, microRNA, and proteome datasets from this species, with data from other vertebrates. In contrast to the platypus, the only other venomous vertebrate with a sequenced genome, we find that snake toxin genes evolve through several distinct co-option mechanisms and exhibit surprisingly variable levels of gene duplication and directional selection that correlate with their functional importance in prey capture. The enigmatic accessory venom gland shows a very different pattern of toxin gene expression from the main venom gland and seems to have recruited toxin-like lectin genes repeatedly for new nontoxic functions. In addition, tissue-specific microRNA analyses suggested the co-option of core genetic regulatory components of the venom secretory system from a pancreatic origin. Although the king cobra is limbless, we recovered coding sequences for all Hox genes involved in amniote limb development, with the exception of Hoxd12. Our results provide a unique view of the origin and evolution of snake venom and reveal multiple genome-level adaptive responses to natural selection in this complex biological weapon system. More generally, they provide insight into mechanisms of protein evolution under strong selection.

Phylogeography of the widespread African puff adder (Bitis arietans) reveals multiple Pleistocene refugia in southern Africa.

Evidence from numerous Pan-African savannah mammals indicates that open-habitat refugia existed in Africa during the Pleistocene, isolated by expanding tropical forests during warm and humid interglacial periods. However, comparative data from other taxonomic groups are currently lacking. We present a phylogeographic investigation of the African puff adder (Bitis arietans), a snake that occurs in open-habitat formations throughout sub-Saharan Africa. Multiple parapatric mitochondrial clades occur across the current distribution of B. arietans, including a widespread southern African clade that is subdivided into four separate clades. We investigated the historical processes responsible for generating these phylogeographic patterns in southern Africa using species distribution modelling and genetic approaches. Our results show that interior regions of South Africa became largely inhospitable for B. arietans during glacial maxima, whereas coastal and more northerly areas remained habitable. This corresponds well with the locations of refugia inferred from mitochondrial data using a continuous phylogeographic diffusion model. Analysis of data from five anonymous nuclear loci revealed broadly similar patterns to mtDNA. Secondary admixture was detected between previously isolated refugial populations. In some cases, this is limited to individuals occurring near mitochondrial clade contact zones, but in other cases, more extensive admixture is evident. Overall, our study reveals a complex history of refugial isolation and secondary expansion for puff adders and a mosaic of isolated refugia in southern Africa. We also identify key differences between the processes that drove isolation in B. arietans and those hypothesized for sympatric savannah mammals.

Medically important snakes and snakebite envenoming in Iran.

Snakebite is a relatively common health condition in Iran with a diverse snake fauna, especially in tropical southern and mountainous western areas of the country with a plethora of snake species. The list of medically important snakes, circumstances and effects of their bite, and necessary medical care require critical appraisal and should be updated regularly. This study aims to review and map the distributions of medically important snake species of Iran, re-evaluate their taxonomy, review their venomics, describe the clinical effects of envenoming, and discuss medical management and treatment, including the use of antivenom. Nearly 350 published articles and 26 textbooks with information on venomous and mildly venomous snake species and snakebites of Iran, were reviewed, many in Persian (Farsi) language, making them relatively inaccessible to an international readership. This has resulted in a revised updated list of Iran's medically important snake species, with taxonomic revisions of some, compilation of their morphological features, remapping of their geographical distributions, and description of species-specific clinical effects of envenoming. Moreover, the antivenom manufactured in Iran is discussed, together with treatment protocols that have been developed for the hospital management of envenomed patients.

One Size Fits All-Venomics of the Iberian Adder (Vipera seoanei, Lataste 1878) Reveals Low Levels of Venom Variation across Its Distributional Range.

European vipers (genus Vipera) are medically important snakes displaying considerable venom variation, occurring at different levels in this group. The presence of intraspecific venom variation, however, remains understudied in several Vipera species. Vipera seoanei is a venomous snake endemic to the northern Iberian Peninsula and south-western France, presenting notable phenotypic variation and inhabiting several diverse habitats across its range. We analysed the venoms of 49 adult specimens of V. seoanei from 20 localities across the species' Iberian distribution. We used a pool of all individual venoms to generate a V. seoanei venom reference proteome, produced SDS-PAGE profiles of all venom samples, and visualised patterns of variation using NMDS. By applying linear regression, we then assessed presence and nature of venom variation between localities, and investigated the effect of 14 predictors (biological, eco-geographic, genetic) on its occurrence. The venom comprised at least 12 different toxin families, of which five (i.e., PLA2, svSP, DI, snaclec, svMP) accounted for about 75% of the whole proteome. The comparative analyses of the SDS-PAGE venom profiles showed them to be remarkably similar across the sampled localities, suggesting low geographic variability. The regression analyses suggested significant effects of biological and habitat predictors on the little variation we detected across the analysed V. seoanei venoms. Other factors were also significantly associated with the presence/absence of individual bands in the SDS-PAGE profiles. The low levels of venom variability we detected within V. seoanei might be the result of a recent population expansion, or of processes other than directional positive selection.

Museum DNA reveals a new, potentially extinct species of rinkhals (Serpentes: Elapidae: Hemachatus) from the Eastern Highlands of Zimbabwe.

Genetic information plays a pivotal role in species recognition and delimitation, but rare or extinct animals can be difficult to obtain genetic samples from. While natural history wet collections have proven invaluable in the description of novel species, the use of these historical samples in genetic studies has been greatly impeded by DNA degradation, especially because of formalin-fixation prior to preservation. Here, we use recently developed museum genomics approaches to determine the status of an isolated population of the elapid snake genus Hemachatus from Zimbabwe. We used multiple digestion phases followed by single strand sequencing library construction and hybridisation capture to obtain 12S and 16S rDNA sequences from a poorly preserved tissue sample of this population. Phylogenetic and morphological analyses in an integrated taxonomic framework demonstrate that the Zimbabwean rinkhals population represents an old and highly distinct lineage, which we describe as a new species, Hemachatus nyangensis sp. nov. Our phylogenetic dating analysis is compatible with venom spitting having evolved in response to the threat posed by early hominins, although more data are required for a robust test of this hypothesis. This description demonstrates the power of museum genomics in revealing rare or even extinct species: Hemachatus from Zimbabwe are only known from a small area of the Eastern Highlands known for high endemism. No living specimens have been seen since the 1980s, most likely due to dramatic land-use changes in the Eastern Highlands, suggesting that the species could be extinct. In view of its recognition as a highly distinct lineage, urgent action is required to determine whether any populations survive, and to safeguard remaining habitat.

Multilevel Comparison of Indian Naja Venoms and Their Cross-Reactivity with Indian Polyvalent Antivenoms.

Snake envenoming is caused by many biological species, rather than a single infectious agent, each with a multiplicity of toxins in their venom. Hence, developing effective treatments is challenging, especially in biodiverse and biogeographically complex countries such as India. The present study represents the first genus-wide proteomics analysis of venom composition across Naja species (N. naja, N. oxiana, and N. kaouthia) found in mainland India. Venom proteomes were consistent between individuals from the same localities in terms of the toxin families present, but not in the relative abundance of those in the venom. There appears to be more compositional variation among N. naja from different locations than among N. kaouthia. Immunoblotting and in vitro neutralization assays indicated cross-reactivity with Indian polyvalent antivenom, in which antibodies raised against N. naja are present. However, we observed ineffective neutralization of PLA2 activities of N. naja venoms from locations distant from the source of immunizing venoms. Antivenom immunoprofiling by antivenomics revealed differential antigenicity of venoms from N. kaouthia and N. oxiana, and poor reactivity towards 3FTxs and PLA2s. Moreover, there was considerable variation between antivenoms from different manufacturers. These data indicate that improvements to antivenom manufacturing in India are highly desirable.

Gene tree parsimony of multilocus snake venom protein families reveals species tree conflict as a result of multiple parallel gene loss.

The proliferation of gene data from multiple loci of large multigene families has been greatly facilitated by considerable recent advances in sequence generation. The evolution of such gene families, which often undergo complex histories and different rates of change, combined with increases in sequence data, pose complex problems for traditional phylogenetic analyses, and in particular, those that aim to successfully recover species relationships from gene trees. Here, we implement gene tree parsimony analyses on multicopy gene family data sets of snake venom proteins for two separate groups of taxa, incorporating Bayesian posterior distributions as a rigorous strategy to account for the uncertainty present in gene trees. Gene tree parsimony largely failed to infer species trees congruent with each other or with species phylogenies derived from mitochondrial and single-copy nuclear sequences. Analysis of four toxin gene families from a large expressed sequence tag data set from the viper genus Echis failed to produce a consistent topology, and reanalysis of a previously published gene tree parsimony data set, from the family Elapidae, suggested that species tree topologies were predominantly unsupported. We suggest that gene tree parsimony failure in the family Elapidae is likely the result of unequal and/or incomplete sampling of paralogous genes and demonstrate that multiple parallel gene losses are likely responsible for the significant species tree conflict observed in the genus Echis. These results highlight the potential for gene tree parsimony analyses to be undermined by rapidly evolving multilocus gene families under strong natural selection.

Domain loss facilitates accelerated evolution and neofunctionalization of duplicate snake venom metalloproteinase toxin genes.

Gene duplication is a key mechanism for the adaptive evolution and neofunctionalization of gene families. Large multigene families often exhibit complex evolutionary histories as a result of frequent gene duplication acting in concordance with positive selection pressures. Alterations in the domain structure of genes, causing changes in the molecular scaffold of proteins, can also result in a complex evolutionary history and has been observed in functionally diverse multigene toxin families. Here, we investigate the role alterations in domain structure have on the tempo of evolution and neofunctionalization of multigene families using the snake venom metalloproteinases (SVMPs) as a model system. Our results reveal that the evolutionary history of viperid (Serpentes: Viperidae) SVMPs is repeatedly punctuated by domain loss, with the single loss of the cysteine-rich domain, facilitating the formation of P-II class SVMPs, occurring prior to the convergent loss of the disintegrin domain to form multiple P-I SVMP structures. Notably, the majority of phylogenetic branches where domain loss was inferred to have occurred exhibited highly significant evidence of positive selection in surface-exposed amino acid residues, resulting in the neofunctionalization of P-II and P-I SVMP classes. These results provide a valuable insight into the mechanisms by which complex gene families evolve and detail how the loss of domain structures can catalyze the accelerated evolution of novel gene paralogues. The ensuing generation of differing molecular scaffolds encoded by the same multigene family facilitates gene neofunctionalization while presenting an evolutionary advantage through the retention of multiple genes capable of encoding functionally distinct proteins.