Monkeying around with venom: an increased resistance to α-neurotoxins supports an evolutionary arms race between Afro-Asian primates and sympatric cobras

Background: Snakes and primates have a multi-layered coevolutionary history as predators, prey, and competitors with each other. Previous work has explored the Snake Detection Theory (SDT), which focuses on the role of snakes as predators of primates and argues that snakes have exerted a selection pressure for the origin of primates’ visual systems, a trait that sets primates apart from other mammals. However, primates also attack and kill snakes and so snakes must simultaneously avoid primates. This factor has been recently highlighted in regard to the movement of hominins into new geographic ranges potentially exerting a selection pressure leading to the evolution of spitting in cobras on three independent occasions. Results: Here, we provide further evidence of coevolution between primates and snakes, whereby through frequent encounters and reciprocal antagonism with large, diurnally active neurotoxic elapid snakes, Afro-Asian primates have evolved an increased resistance to α-neurotoxins, which are toxins that target the nicotinic acetylcholine receptors. In contrast, such resistance is not found in Lemuriformes in Madagascar, where venomous snakes are absent, or in Platyrrhini in the Americas, where encounters with neurotoxic elapids are unlikely since they are relatively small, fossorial, and nocturnal. Within the Afro-Asian primates, the increased resistance toward the neurotoxins was significantly amplified in the last common ancestor of chimpanzees, gorillas, and humans (clade Homininae). Comparative testing of venoms from Afro-Asian and American elapid snakes revealed an increase in α-neurotoxin resistance across Afro-Asian primates, which was likely selected against cobra venoms. Through structure-activity studies using native and mutant mimotopes of the α-1 nAChR receptor orthosteric site (loop C), we identified the specific amino acids responsible for conferring this increased level of resistance in hominine primates to the α-neurotoxins in cobra venom. Conclusion: We have discovered a pattern of primate susceptibility toward α-neurotoxins that supports the theory of a reciprocal coevolutionary arms-race between venomous snakes and primates.


Background
Snakes and primates have a sophisticated coevolutionary history as predators, prey, and competitors with each other. Prior work has focused on two distinct areas: (1) the Snake Detection Theory (SDT), in which encounters with snakes have shaped primate evolution, and (2) how primate diversification has shaped the evolution of venomous snakes' defensive behaviors and morphological innovations.
The SDT of primate evolution encapsulates coevolutionary interactions between primates and snakes that have shaped primate neurobiology, psychology, and physiology [1,2]. This theory proposes that constricting snakes selected for changes in visual systems that led to the differentiation of primates as a separate order of mammals. It also suggests that further visual system expansion later evolved in anthropoids, including Catarrhini (African/Asian monkeys, apes, and humans) and Platyrrhini (American monkeys), as a response to the evolution of venomous snakes in the families Elapidae and Viperidae [1,2]. Predation by venomous snakes would become rarer events as primate body size increased, but defensive envenomation by these snakes would still exert a strong selection pressure, according to this theory. In support of the SDT, a visual bias leading to faster or more reliable detection of snakes than to other animals/objects has been found in those primates tested, i.e., catarrhines, including Homo sapiens [3][4][5][6][7]. However, the SDT suggests that rapid and reliable detection might be variable among primates because visual adaptations differ among major lineages and coincide with their evolutionary co-existence with venomous snakes. For example, primates of the suborder Strepsirrhini (Lorisiformes, Lemuriformes) have poorer central vision and acuity systems than the suborder Haplorhini (Tarsiiformes, Platyrrhini, and Catarrhini) because they do not have a true fovea [1,2]. Lemuriformes colonized Madagascar between 50 and 65 MA [8], likely before the evolution of vipers~50 MA [9], but certainly before their pan Afro-Asian radiation, and long before the evolution of elapid snakes and their characteristic neurotoxic venom~37 MA [10]. Madagascar has remained devoid of venomous snakes; thus, the SDT proposes that there has never been a selection pressure on this group to evolve survival adaptations against venomous snakes. Although platyrrhines and catarrhines are anthropoid primates, which are thought to have evolved in Asia or Africa, platyrrhines have more variable visual systems than catarrhines [1,2]. The platyrrhines diverged from other anthropoids after they dispersed to South Americã 35 MA [11], long before the arrival of venomous snakes in South America between 3 and 25 MA, based on both molecular data and the formation of the Panamanian land bridge [9,[12][13][14]. Thus, unlike catarrhines, platyrrhines have had an interrupted evolutionary time scale with venomous snakes [1,2]. Most platyrrhines are highly arboreal, while the recently sympatric venomous snakes are mostly terrestrial, with neurotoxic elapid species in particular being fossorial and crepuscular/nocturnal. Thus, selection pressures from venomous snakes are also expected to be weaker on platyrrhines than on catarrhines.
Conversely, there is evidence that venomous elapids have evolved some adaptations in response to primate diversification. Spitting in cobras has evolved on three convergent occasions (once in the genus Hemachatus and on two independent occasions within Naja) as a defensive trait that causes intense ocular pain and inflammation, which might have evolved in response to unique behavioral adaptations by anthropoids, particularly early hominines [15]. Predation events by natural predators (e.g., birds, mammals, and other reptiles) of spitting cobras often occur. Thus, it seems that spitting is not an effective method against most of their natural predators [15]. Intriguingly, the emergence of spitting in Hemachatus (< 17 MA), the African spitting Naja (~6.7 MA), and the Asian spitting Naja (~2.5 MA) [10], all coincided with the diverse emergence and radiations of terrestrial dwelling apes, including early hominines [16,17], all of which occurred in habitats also frequented by largely terrestrial spitting cobras [18]. These hominines showed a myriad of features indicating a tendency toward upright bipedal movement and increased precision and opposability of the hands and thumbs. It has been inferred that the last common ancestors of chimpanzees and humans already used wood and stone tools and had the ability to hunt vertebrates [17]. Furthermore, terrestrial dwelling apes, through their locomotor style, foraging habits, and group movements, were more likely than other terrestrial taxa to come into frequent contact with venomous snakes [19]. The emergence of Asian spitting cobras at~2.5 MA [10] coincided with the appearance of Homo erectus in Asia (~1.8 MA), the latter of which is characterized by advancements in stone tool technology [20]. Primates in general show overt reactions to snakes, including snake-specific alarm calls, and may mob snakes with improvised tools and projectiles [21][22][23]. Thus, the use of projectile weaponry in primates, particularly in anthropoids that stand upright when using such weapons, would provide a selection pressure on snakes to evolve a ranged defense such as spitting venom at the eyes of upright primates. Additionally, the evolution of injury-inducing defensive cytotoxins is also correlated with that of defensive hooding displays and aposematic marking in African and Asian cobras [24]. The evolution of defensive cytotoxins also preceding spitting suggests that a high selection pressure might have been exerted upon cobras for defense against primates before the emergence of spitting when hominins evolved bipedalism. Hence, there is some fundamental evidence to consider that reciprocal coevolution occurred between primates and snakes, particularly venomous snakes, driving novel traits in both taxa.
Non-primate mammals that co-exist with venomous snakes have some resistance toward snake venom [25][26][27][28][29][30], while populations absent of venomous snakes do not [31]. A prominently evolved resistance mechanism is seen at the orthosteric site (acetylcholine (ACh) binding region loop C [32][33][34]) of the α-1 nicotinic acetylcholine receptor (nAChR) by which specific amino acid residues reduce the binding of α-neurotoxins (such as threefinger toxins; 3FTxs) found within some snakes venoms, particularly that of Elapidae [35,36]. The α-1 nAChR is located at the neuromuscular junction, ultimately controlling muscle contraction and function, and thus is a prime target for many toxins that cause paralysis, such as predatory venoms. The α-1 nAChR orthosteric site of humans has been shown to be much less susceptible to binding of α-neurotoxins [37][38][39][40] in contrast to the same site binding of other taxa [41][42][43][44][45]. Thus, while neurotoxic effects of elapid venoms may be a lethal outcome in some human envenomations (but not always), it is possible that in these instances, binding allosterically (to a receptor region other than the orthosteric site) might be a more derived method [46], rather than binding to the orthosteric site.
Despite this increased resistance in H. sapiens and given the mounting evidence in support of their coevolution with venomous snakes, no studies have attempted to investigate how α-neurotoxin binding might affect a wide range of primate species. Afro-Asian primates coexist with large, abundant, diurnally active cobras that readily defend themselves, leading to death. Thus, we set out to test the hypothesis that if resistance to these paralytic toxins evolved, it would be present in Afro-Asian primates and lacking in primates from other regions that either lack venomous snakes (Madagascar) or where neurotoxic elapid snakes are small, fossorial, and crepuscular/nocturnal and encounters are rare (the Americas).

Results
Investigations into the publicly available α-1 nAChR sequences of a range of primate species revealed that the orthosteric site loop C region of H. sapiens is conserved across the Homininae (chimpanzees, gorillas, humans), while other primates such as Cercopithecidae (African/ Asian monkeys), Galagidae (bushbabies), Lemuriformes (lemurs), Platyrrhini (American monkeys) and Ponginae (orangutans) all have clade-specific orthosteric sequences (Table S1). Given these differences among primate clades and the theory that primates and snakes might have coevolved certain traits through reciprocal antagonism, we tested if there were any differences in αneurotoxin binding between primate lineages, particularly given evidence that H. sapiens nAChRs are less susceptible to α-neurotoxins [37][38][39][40]. To assess these differences, we tested seven venoms that are rich in αneurotoxins (3FTxs): five from species of cobra (Naja spp.) that inhabit Africa (N. haje, N. mossambica, and N. nubiae) and Asia (N. kaouthia and N. siamensis) and two species of Micrurus from South America (M. browni and M. corallinus). We further categorized the cobras into spitting (N. mossmabica, N. nubiae, and N. siamensis) and non-spitting (N. haje and N. kaouthia) species. We tested the binding of these venoms upon mimotopes that represent the orthosteric site loop C region of α-1 nAChRs from different primate clades (Table S1). We utilized a biolayer interferometry (BLI) assay, which has been previously validated to assess the binding of αneurotoxins to taxa-specific orthosteric mimotopes [37,[46][47][48], to determine if primates might have evolved some resistance elements toward α-neurotoxins.
The results indicate a consistent binding pattern across all cobra venoms tested toward the primate lineage mimotopes (Figs. 1 and 2; Additional file 1: Figs S1-S3). The non-spitting species seem to have a much lower binding overall than the spitting, which would suggest lower proportions of orthosteric-site binding αneurotoxic 3FTxs within their venom [47]. The difference observed in binding of the crude venoms is a consequence of the toxin interactions with the different biochemical properties of each amino acid present within the sequences [39,[49][50][51]. The binding patterns across all venoms tested showed a significantly increased resistance toward α-neurotoxins across the Homininae (Figs. 1 and 2; Additional file 1: Figs S1-S3). These data additionally support previous observations of weakly binding α-neurotoxins toward human nAChRs [37][38][39][40]. Further, the venoms also bound relatively weakly to the Cercopithecidae and Ponginae mimotopes (Figs. 1 and 2; Additional file 1: Figs S1-S3). The three clades, Cercopithecidae, Homininae, and Ponginae, occur in Africa and Asia, coinciding with the geographical distribution of cobras. The representatives of Lorisiformes and Tarsiiformes also have relatively low susceptibility, and both of these clades also occupy Africa and/or Asia. Conversely, the Lemuriformes (with the exception of Eulemur flavifrons) and Platyrrhini were the most susceptible to binding by the venoms (Figs. 1 and 2; Additional file 1: Figs S1-S3). Both of these clades are geographically separated from the African and Asian primate clades and also from diurnal elapids.
To further support the hypothesis that primates have evolved adaptations to reduce mortality from cobra encounters in Africa and Asia, we tested the binding of venoms from two species of American coral snakes (Micrurus browni and M. corallinus). The data show that there is no significance in binding between all primate groups of M. browni and a very small difference in binding between the clades of M. corallinus (Figs. 3 and 4) in comparison to that of cobras (Figs. 1 and 2; Additional file 1: Figs S1-S3). These data suggest that since there is a greater disparity between binding of both Platyrrhini and Lemuriformes for cobras but very little or no significant binding between these clades for coral snakes, then the strong selection pressure for the resistance across the African and Asian primates is likely from cobras.
In order to ascertain the structure-activity relationships leading to hominines having the most reduced level of sensitivity to α-neurotoxins, we constructed B Bar graphs represent the mean AUC values of the adjacent curve graphs. Statistical significance is indicated by matching letters with the colors of letter indicating the level of significance; brown p< 0.001, green p< 0.01, pink p< 0.05. C Curve graphs show the mean wavelength shift (nm) in light with binding of venoms over a 120-s association phase. The venom was tested in triplicate (n=3). Error bars on all graphs represent the SEM. AUC values were statistically analyzed using a one-way ANOVA with a Tukey's comparisons multiple comparisons test comparing to the native mimotope. All raw data and statistical analyses outputs can be found in Additional file 2 mimotope mutants to elucidate key amino acid positions in the native sequence that are responsible for resistance. The bindings of representative African and Asian cobra venoms were tested against a series of orthosteric loop C site mutants with amino acid substitutions of the native hominin sequences. Each amino acid substituted was for that of the ancestral amino acid (Fig. 5). The results indicate that the amino acid positions 187 and 189 together seem to conform the bulk of the increased resistance since there was a significant increase in binding when these positions were substituted for the ancestral versions. These results were consistent for the representative African and Asian cobra venoms (Fig. 5). The amino acid substitutions within the Homininae sequence Statistical significance is indicated by matching letters with the colors of letter indicating the level of significance; brown p< 0.001, green p< 0.01, pink p< 0.05. C Curve graphs show the mean wavelength shift (nm) in light with binding of venoms over a 120-s association phase. The venom was tested in triplicate (n=3). Error bars on all graphs represent the SEM. AUC values were statistically analyzed using a one-way ANOVA with a Tukey's comparisons multiple comparisons test comparing to the native mimotope. All raw data and statistical analyses outputs can be found in Additional file 2 are biochemically polar, while the ancestral substitutions are non-polar. These polar amino acids likely provide enough of an electrostatic static repulsion in some instances to reduce the effective binding of the charged surface of the 3FTxs within the venoms. These results are also consistent with the previously untested hypothesis that lower binding of cobra venoms to humans is likely due to combinational substitutions at the aromatic subsite (positions 187 and 189) [39]. Further, there was an increase in susceptibility with all three amino acids at positions 187, 189, and 191 removed. All three of these positions are key binding sites for α-neurotoxins [28,39,50,51], further confirming the importance of these positions in the binding of these toxin types to nAChRs.

Discussion
All these data combined suggest that there was a strong selection pressure that has led to African and Asian primates becoming less susceptible toward cobra αneurotoxin that bind to the orthosteric site, while B Bar graphs represent the mean AUC values of the adjacent curve graphs. C Curve graphs show the mean wavelength shift (nm) in light with binding of venoms over a 120-s association phase. The venom was tested in triplicate (n=3). Error bars on all graphs represent the SEM. No statistical significance was detected using a one-way ANOVA. All raw data and statistical analyses outputs can be found in Additional file 2 primates living in other geographical areas were not subjected to these same selection pressures. The emerging patterns are consistent with the SDT [1] in that the primate clades that evolved in Africa and/or remained sympatric to venomous snakes show mechanisms to cope with certain selective pressures. Catarrhini, Lorisiformes, and Tarsiiformes have evolved an increased resistance toward α-neurotoxins whereas the Malagasy Lemuriformes have never co-existed with venomous snakes and the Platyrrhini occupy Central and South America where the elapids are small, crepuscular/nocturnal and fossorial, and thus pose little threat. The data further support the hypothesis that cobras and early hominins have been entwined in a coevolutionary armsrace resulting in the evolution of specific traits for both clades [15]. From an evolutionary perspective, it would seem that increased resistance convergently evolved on three occasions: once within Tarsiiformes (Carlito syrichta), once within Catarrhini (with resistance further amplified within Homininae) and once within Lorisiformes (Otolemur garnettii) (Figs. 1 and 2; Additional file 1: Figs S1-S3). Tarsiiformes and Lorisiformes likely evolved an increased resistance since they are small and prospective prey and in some cases predators to venomous snakes, whereas the larger Afro-Asian primates likely evolved the increased resistance due to an increase in terrestriality thus potentiating the number of possible encounters leading to defensive bites. Lemuriformes and Platyrrhini were never in contact with large venomous elapids as previously mentioned, and thus have never had a selection pressure to evolve any form of resistance. In the case of E. flavifrons, it is uncertain why there seems to be an increased resistance for this species while the Fig. 5. The effects of venom from a representative African cobra species (Naja mossambica) and a representative Asian cobra species (N. siamensis) against native and mutant Homininae mimotope sequences. A Amino acid sequences of native and mutant mimotopes. B Bar graphs represent the mean area under the curve (AUC) values of the adjacent curve graphs. Curve graphs show the mean wavelength (nm) shift in light with increased binding of venoms over a 120-s association phase. Each venom was tested in triplicate (n=3). Error bars on all graphs represent the SEM. AUC values were statistically analyzed using a one-way ANOVA with a Dunnett's multiple comparisons post hoc test comparing to the native mimotope. A statistical significance is annotated above bars by ** (p< 0.01), *** (p< 0.001), or **** (p< 0.0001). All raw data and statistical analyses outputs can be found in Additional file 2 other Lemuriformes did not evolve any since there is a lack of selection pressure. Conspicuously some dietary items of Eulemur, such as mushrooms and millipedes, contain high levels of neurotoxins [52]. Future work should investigate if the toxins from these sources target nAChRs (e.g., high concentrations of nicotine) and provide enough of a selection to maintain the lower susceptibility toward nAChR targeting toxins. There is likely an evolutionary trade-off disadvantage to fitness with evolving specific residue changes that lowers α-neurotoxin binding such as a reduced binding of the endogenous ACh [53]; thus, there must be an evolutionary selective advantage for Eulemur to have evolved this when the rest of the Lemuriformes did not.
Primates are largely an arboreal order, and the period during which spitting in cobras evolved occurred during a major climatic shift linked to a reduction in continuous forests [17]. During this time, diurnal primates would have needed to descend to the ground, and a number of species adapted to be fully terrestrial, increasing their chances of encountering and being defensively bitten by venomous elapid snakes. It should be noted that Ponginae are slightly less resistant than the Homininae clade, and this is possibly due to them being more arboreal than Homininae, reinforcing the importance of terrestriality in being exposed to elapids. Considering that snakes are can be prey items for some primates, and venomous snakes are a danger if startled and defensively bite, it is argued that snakes placed a tremendous selection pressure on anthropoid primates, to visually detect and avoid venomous snakes [1,2].
One aspect of the SDT proposes that nocturnal strepsirrhines in Africa and Asia should have venom resistance since they were unable to expand their visual systems to the same extent as diurnal primates [2], which our data supports. However, this hypothesis went one step further suggesting that they would also have a greater venom resistance than other primates, which according to our data is not quite the case, but not entirely far from our observations either. Indeed, Lorisiformes and Tarsiiformes have relatively low susceptibility, with both of these clades also occupying Africa and/or Asia. They both also have particularly high orbital convergence, which has been associated with their coevolution with snakes [1]. Both of these clades consist of relatively small primates that are potential prey items for large venomous snakes. However, within the Lorisidae, the Sunda slow loris (Nycticebus coucang) contains a different type of resistance, that of the N-glycosylation motif (Table S1), which has been associated with extremely effective α-neurotoxin resistance [54][55][56][57]. This trait might have evolved due to a different selection pressure than that of Otolemur garnettii, since slow lorises have been documented eating venomous snakes. Slow lorises are also venomous themselves [58,59], but given the extreme cytotoxicity of their venom and lack of neurotoxic symptoms, the N-glycosylation motif is unlikely to have evolved for autoresistance.
The greater increase in resistance within hominins coincides not only with obligate bipedalism and terrestriality, but also with the reduction of continuous forests requiring their movements into more open habitats. Foraging and traveling hominins were likely to come into contact much more often with terrestrial, diurnally active neurotoxic elapid snakes such as cobras, which would defend themselves against perceived threats using their venom [22]. The likely increased resistance toward α-neurotoxins from cobras might have driven cobras to develop a different kind of defense against hominins in the form of cytotoxins and spitting [15,24]. An evolved form of partial resistance meant that cobra venom was not an effective method to reduce hominin interactions thus cytotoxins and spitting might have evolved in response as a deterrent. It has been proposed that spitting cobra clades evolved in response to hominines, and this might explain the difference we observed in binding between spitting and non-spitting cobras. A lower orthosteric site binding of α-neurotoxins in the venom of nonspitting than spitting cobras (likely due to the proportion of α-neurotoxins specific for this site within the venom) might be that the increased resistance across African and Asian primates has meant that spitters have had to evolve a high proportion of α-neurotoxins to overcome this increased resistance. Yet this proposed hypothesis needs careful additional consideration and investigation.
Further, since we have highlighted the selective coevolutionary pressures that might have come through hominin-snake interactions, it is entirely plausible that α-neurotoxin susceptibility differs between populations within humans. For example, populations that are highly prone to venomous snake conflicts and high rates of envenomations, such as Sub-Saharan Africa or rural India, may have a greater level of partial resistance than populations that never encounter wild elapids, such as within Europe. Though hypothetical, further research into variations of the orthosteric site sequences between different human populations could allow us to understand the effects of envenomations between different populations and aid our understanding of the snakebite crisis.
Thus, given that commonly encountered neurotoxic snakes exert a strong selection pressure in their defensive envenomations, why does there only seem to be an increased resistance rather than full resistance? Firstly, our aforementioned proposition that there might be an evolutionary trade-off disadvantage to evolving full resistance at the nAChR orthosteric site [53,60] might mean that evolving partial resistance could allow for a balance between reduced α-neurotoxin susceptibility and a somewhat efficient ACh binding. This is supported by previous work revealing that increased resistance toward the nAChR agonist epibatidine decreased ACh sensitivity to the site [53]. Further, some vipers that are prey to αneurotoxic snakes often display resistance, but this resistance is secondarily lost in populations which have dispersed into geographic areas devoid of neurotoxic snakes [57]. This differential resistance between sympatric and allopatric populations of venomous snakes and prey is also paralleled in ground squirrels and their rattlesnake predators, albeit a different mechanism of resistance [31]. Secondly, it might also be that partial resistance allows for a greater chance of survival against low levels of α-neurotoxins, since venomous snakes can control their injection volumes, often producing lower venom volumes during defensive bites to conserve the energetically costly venom for prey items rather than wasting it on a non-prey item [61][62][63]. Further to this point, human envenomations do not always result in death (even without antivenom intervention), thus partial resistance could also provide enough of an advantage to aid in survival in enough envenomation instances that the reduced susceptibility trait is effective enough to be maintained within a population. Thirdly, since primates are largely social animals [64,65], it is possible that increased resistance might have been successful enough for kin selection [66][67][68], to maintain the partial resistance trait reciprocally. Having an increased resistance would allow the envenomed individual to fight against the snake for longer than more susceptible primates and also warn other members of the group that a venomous snake is present. Warning calls toward snakes occur within both nocturnal and diurnal primate groups [69,70], thus being able to marginally extend an individual's survival time to warn the group, allowing the safety of the relatives and, ultimately, their successful reproduction. This kin selection hypothesis certainly has its merits, while also raising more questions in terms of exactly how kin selection would be favored in light of partial resistance.
Although we have provided potential hypotheses for the coevolution of increased resistance within hominins that is well supported by data in this study and previous work, a competing hypothesis is that the resistance is due to simple genetic drift rather than being the result of evolutionary selection pressures. Under the evolutionary drift hypothesis, early hominins and primates in general tended to occur in small groups where genetic drift is most effective, thus the trait may have been maintained through this scenario. However, given the consilience of evidence suggesting a deeper coevolutionary history between snakes and primates, we think that our aforementioned hypotheses might hold a greater standing than something as innocuous as the genetic drift hypothesis. This thus provides a foundation upon which to conduct additional research to further strengthen and explore this fascinating evolutionary hypothesis.

Conclusion
We have revealed an increased resistance across Afro-Asian primates toward snake venom α-neurotoxins that target nicotinic acetylcholine receptor orthosteric sites that seems to have convergently evolved on three occasions. There was a further amplification of resistance within the Homininae (chimpanzee, gorilla, human clade) when tested across African and Asian spitting and non-spitting cobras. This finding is consistent with clinical symptoms of neurotoxicity being reduced for αneurotoxins that bind to the orthosteric site, relative to the natural prey types. Conversely, the Lemuriformes and Platyrrhini have the highest susceptibility of primates tested toward α-neurotoxic venoms. This coincides with the evolutionary biogeography of these lineages in Madagascar and South/Central America respectively, evolving in the absence of large neurotoxic snake species. Our data also further support one of the hypotheses proposed by the Snake Detection Theory, whereby venomous snakes and primates share a history of coevolution through their continuous co-existence across Africa and Asia [1], as evidenced by the venom of two species of coral snakes tested here show little or no disparity between primate group mimotope binding. These data further support the hypothesis that cobras and hominins have evolved certain traits through a coevolutionary arms race. The evidence in support of the coevolution of primates and snakes seems to be increasing, yet despite this, fundamental gaps in knowledge remain both regarding the ecology of snakes and primates and how their shared coevolutionary history might have brought about some of their most distinguishing traits. Future work should investigate why partial rather than full resistance has evolved, including our hypothesis of kin selection.

Venom collection and preparation
All venom work was undertaken under the auspices of UQ IBSC approval #IBC134BSBS2015. Venoms were pools of adult snakes (N=3) to minimize the effect of individual variation from the long-term cryogenic collection of the Venom Evolution Lab. All venom samples were lyophilized and reconstituted in double deionized water (ddH 2 O), and centrifuged (4°C, 10 min at 14,000 relative centrifugal force (RCF)). The supernatant was made into a working stock (1 mg/mL) in 50% glycerol at − 20°C. The concentrations of experimental stocks were determined using a NanoDrop 2000 UV-Vis Spectrophotometer (Thermo Fisher, Sydney, Australia) at an absorbance wavelength of 280 nm.
The C-C of the native mimotope is replaced during peptide synthesis with S-S to avoid uncontrolled postsynthetic thiol oxidation. The C-C bond in the nAChR binding region does not participate directly in analyteligand binding [51,71,72], thus replacement to S-S is not expected to have any effect on the analyte-ligand complex formation. However, the presence of the C-C bridge is key in the conformation of the interaction site of whole receptors [73]. As such, we suggest direct comparisons of kinetics data, such as Ka or KD, between nAChR mimotopes and whole receptor testing should be avoided, or at least approached with caution. Mimotopes were further synthesized to a biotin linker bound to two aminohexanoic acid (Ahx) spacers, forming a 30-Å linker. Mimotope dried stocks were solubilized in 100% dimethyl sulfoxide (DMSO) and diluted in ddH 2 O at 1:10 dilution to obtain a stock concentration of 50 μg/mL. Stocks were stored at − 80°C until required.

Biolayer Interferometry (BLI)
Full details of the developed assay, including all methodology and data analysis, can be found in the validated protocol [48] and further data using this protocol [37,47]. In brief, the BLI assay was performed on the Octet HTX system (ForteBio, Fremont, CA, USA). Venom samples were diluted 1:20, making a final concentration of 50 μg/mL per well. Mimotope aliquots were diluted 1: 50, with a final concentration of 1 μg/mL per well. The assay running buffer was 1X DPBS with 0.1% BSA and 0.05% Tween-20. Preceding experimentation, Streptavidin biosensors were hydrated in the running buffer for 30-60 min, while on a shaker at 2.0 revolutions per minute (RPM). The dissociation of analytes occurred using a standard acidic glycine buffer solution (10 mM glycine (pH 1.5-1.7) in ddH 2 O). Raw data are provided in Additional file 1.

Data processing and analysis
All data obtained from BLI on Octet HTX system (For-teBio) were processed in exact accordance to the validation of this assay [48]. The association step data were obtained and imported into Prism8.0 software (GraphPad Software Inc., La Jolla, CA, USA) where area under the curve (AUC) and one-way ANOVA with Tukey's multiple comparisons analyses were conducted and graphs produced.
Phylogenetic trees were obtained from timetree.org. The obtained phylogenetic trees were then further analyzed in RStudio (R Core Team, 2015) for all comparative analyses using the Ape package [74]. Heat-mapping of AUC values over the phylogenetic trees was achieved using the con-tMap function of the R package phytools [75].
Additional file 1: Fig. S1. The effects of venom from an additional African cobra species (Naja nubiae) against the nAChR orthosteric site mimotopes from seven clades of primates. Fig. S2. The effects of venom from an additional African cobra species (Naja haje) against the nAChR orthosteric site mimotopes from seven clades of primates. Fig. S3. The effects of venom from an additional Asian cobra species (Naja kaouthia) against the nAChR orthosteric site mimotopes from seven clades of primates. Table S1. Orthosteric site sequences for each species and public database accessions codes.
Additional file 2:. All the raw experimental data output and statistical analyses. Availability of data and materials All data generated or analyzed during this study are included in this published article and Additional files 1 and 2.

Declarations
Ethics approval and consent to participate As no live animals were used in this project, this is not applicable.

Consent for publication
Not applicable