Evolution of snake venom

Venom in snakes and some lizards is a form of saliva that has been modified into venom over its evolutionary history.[1] In snakes, venom has evolved to kill or subdue prey, as well as to perform other diet-related functions.[2] While snakes occasionally use their venom in self defense, this is not believed to have had a strong effect on venom evolution.[3] The evolution of venom is thought to be responsible for the enormous expansion of snakes across the globe.[4][5][6]

The rattlesnake Crotalus oreganus eating its prey, which it uses its venom to subdue

The evolutionary history of snake venom is a matter of debate. Historically, snake venom was believed to have evolved once, at the base of the Caenophidia, or derived snakes. Molecular studies published beginning in 2006 suggested that venom originated just once among a putative clade of reptiles, called Toxicofera, approximately 170 million years ago.[7] Under this hypothesis, the original toxicoferan venom was a very simple set of proteins that were assembled in a pair of glands. Subsequently, this set of proteins diversified in the various lineages of toxicoferans, including Serpentes, Anguimorpha, and Iguania: several snake lineages also lost the ability to produce venom.[8][9] The Toxicoferan hypothesis was challenged by studies in the mid-2010s, including a 2015 study which found that venom proteins had homologs in many other tissues in the Burmese python.[10][11] The study therefore suggested that venom had evolved independently in different reptile lineages, including once in the Caenophid snakes.[10] Venom containing most extant toxin families is believed to have been present in the last common ancestor of the Caenophidia: these toxins subsequently underwent tremendous diversification, accompanied by changes in the morphology of venom glands and delivery systems.[12]

Snake venom evolution is thought to be driven by an evolutionary arms race between venom proteins and prey physiology.[13] The common mechanism of evolution is thought to be gene duplication followed by natural selection for adaptive traits.[14] The adaptations produced by this process include venom more toxic to specific prey in several lineages,[15][16][17] proteins that pre-digest prey,[18] and a method to track down prey after a bite.[19] These various adaptations of venom have also led to considerable debate about the definition of venom and venomous snakes.[20] Changes in the diet of a lineage have been linked to atrophication of the venom.[8][9]

Evolutionary history

edit

The origin of venom is thought to have provided the catalyst for the rapid diversification of snakes in the Cenozoic period,[21] particularly to the Colubridae and their colonization of the Americas.[5] Scholars suggest that the reason for this huge expansion was the shift from a mechanical to a biochemical method of subduing prey.[22][23] Snake venoms attack biological pathways and processes that are also targeted by venoms of other taxa; for instance, calcium channel blockers have been found in snakes, spiders, and cone snails, thus suggesting that venom exhibits convergent evolution.[24] Venom is common among derived snake families.[13] Venom containing most extant toxin families is believed to have been present in the last common ancestor of the Caenophidia, also called Colubroidea. These toxins subsequently underwent tremendous diversification, accompanied by changes in the morphology of venom glands and delivery systems. This diversification is linked to the rapid global radiation of the advanced snakes.[12] The tubular or grooved fangs snakes use to deliver their venom to their target have evolved multiple times, and are an example of convergent evolution. The tubular fangs common to front-fanged snakes are believed to have evolved independently in Viperidae, Elapidae, and Atractaspidinae.[25]

Serpentes

A cladogram adapted from Fry et al. (2012) showing a subset of suggested protein recruitment events. [A]: 13 toxin families, including 3FTx and metalloprotease. [B]: 2 toxin families, including PLA2 Type IIA and the P-II class of snake venom metalloproteases [C]: 2 toxin families, including PLA2 Type IB.[26][27]

Until the use of gene sequencing to create phylogenetic trees became practical, phylogenies were created on the basis of morphology. Such traditional phylogenies suggested that venom originated along multiple branches among Squamata approximately 100 million years ago: in the Caenophidia, or derived snakes, and in the lizard genus Heloderma.[7][10][28] Studies using nuclear gene sequences in the mid-2000s and early 2010s found the presence of venom proteins in the lizard clades Anguimorpha and Iguania similar to those of snakes, and suggested that together with Serpentes, these formed a clade, which they named "Toxicofera".[7][28] This led to the theory that venom originated only once within the entire lineage approximately 170 million years ago. This ancestral venom was described as consisting of a very simple set of proteins, assembled in a pair of glands.[7] The venoms of the different lineages then diversified and evolved independently, along with their means of injecting venom into prey.[7] This diversification included the independent evolution of front-fanged venom delivery from the ancestral rear-fanged venom delivery system.[29] The single origin hypothesis also suggests that venom systems subsequently atrophied, or were completely lost, independently in a number of lineages.[20] The phylogenetic position of Iguania within Toxicofera is supported by most molecular studies, but not by morphological ones.[30]

The "Toxicoferan hypothesis" was subsequently challenged. A study performed in 2014 found that homologs of 16 venom proteins, which had been used to support the single origin hypothesis, were all expressed at high levels in a number of body tissues.[11] The authors therefore suggested that previous research, which had found venom proteins to be conserved across the supposed Toxicoferan lineage, might have misinterpreted the presence of more generic "housekeeping" genes across this lineage, as a result of only sampling certain tissues within the reptiles' bodies. Therefore, the authors suggested that instead of evolving just once in an ancestral reptile, venom evolved independently in multiple lineages, including once prior to the radiation of the "advanced" snakes.[11] A 2015 study found that homologs of the so-called "toxic" genes were present in numerous tissues of a non-venomous snake, the Burmese python. One of the authors stated that they had found homologs to the venom genes in many tissues outside the oral glands, where venom genes might be expected. This demonstrated the weaknesses of only analyzing transcriptomes (the total messenger RNA in a cell).[10] The team suggested that pythons represented a period in snake evolution before major venom development. The researchers also found that the expansion of venom gene families occurred mostly in highly venomous caenophidian snakes (also referred to as "colubrid snakes"), thus suggesting that most venom evolution took place after this lineage diverged from other snakes.[10] The debate over the Toxicoferan hypothesis is driven in part by disagreement over the definition of a venom.[10][31] As of 2022, the Toxicoferan hypothesis remains a prevalent view.[32]

Mechanisms of evolution

edit
 
Phospholipase A2, an enzyme found in normal tissue that has been adapted in certain snake venoms. The example pictured here is found in bee stings.

The primary mechanism for the diversification of venom is thought to be the duplication of gene coding for other tissues, followed by their expression in the venom glands. The proteins then evolved into various venom proteins through natural selection. This process, known as the birth-and-death model, is responsible for several of the protein recruitment events in snake venom.[33][13] These duplications occurred in a variety of tissue types with a number of ancestral functions. Notable examples include 3FTx, ancestrally a neurotransmitter found in the brain, which has adapted into a neurotoxin that binds and blocks acetylcholine receptors. Another example is phospholipase A2 (PLA2) type IIA, ancestrally involved with inflammatory processes in normal tissue, which has evolved into venom capable of triggering lipase activity and tissue destruction.[21] The change in function of PLA2, in particular, has been well documented; there is evidence of several separate gene duplication events, often associated with the origin of new snake species.[34] Non-allelic homologous recombination induced by transposon invasion (or recombination between DNA sequences that are similar, but not alleles) has been proposed as the mechanism of duplication of PLA2 genes in rattlesnakes, as an explanation for its rapid evolution.[35] These venom proteins have also occasionally been recruited back into tissue genes.[36]

Gene duplication is not the only way that venom has become more diverse. There have been instances of new venom proteins generated by alternative splicing. The Elapid snake Bungarus fasciatus, for example, possesses a gene that is alternatively spliced to yield both a venom component and a physiological protein.[36] Further diversification may have occurred by gene loss of specific venom components. For instance, the rattlesnake ancestor is believed to have had the PLA2 genes for a heterodimeric neurotoxin now found in Crotalus scutulatus, but those genes are absent in modern non-neurotoxic Crotalus species; the PLA2 genes for the Lys49-myotoxin supposedly existing in the common ancestor of rattlesnakes were also lost several times on recent lineages to extant species [35] Domain loss has also been implicated in venom neofunctionalization. Investigation of the evolutionary history of viperid SVMP venom genes revealed repeated occasions of domain loss, coupled with significant positive selection in most of the phylogenetic branches where domain loss was thought to have occurred.[37] Venom toxins have also evolved via the gene "hijacking" or "co-opting", or the change in function of unrelated genes. A 2021 study suggested that co-opting explained the evolution of most types of toxins, but not that of the toxins that are most abundant in snake venom.[27]

Protein recruitment events have occurred at different points in the evolutionary history of snakes. For example, the 3FTX protein family is absent in the viperid lineage, suggesting that it was recruited into snake venom after the viperid snakes branched off from the remaining colubroidae.[38] PLA2 is thought to have been recruited at least two separate times into snake venom, once in elapids and once in viperids, displaying convergent evolution of this protein into a toxin.[39][40] A 2019 study suggested that gene duplication could have allowed different toxins to evolve independently, allowing snakes to experiment with their venom profiles and explore new and effective venom formulations.[40] This was proposed as one of the ways snakes have diversified their venom formulations through millions of years of evolution.[40] The various recruitment events had resulted in snake venom evolving into a very complex mixture of proteins. The venom of rattlesnakes, for example, includes nearly 40 different proteins from different protein families,[41] and other snake venoms have been found to contain more than 100 distinct proteins.[22] The composition of this mixture has been shown to vary geographically, and to be related to the prey species available in a certain region.[17] Snake venom has generally evolved very quickly, with changes occurring faster in the venom than in the rest of the organism.[42]

Selection pressure

edit

Long-standing hypotheses of snake venom evolution have assumed that most snakes inject far more venom into their prey than is required to kill them; thus, venom composition would not be subject to natural selection. This is known as the "overkill" hypothesis.[43] However, recent studies of the molecular history of snake venom have contradicted this, instead finding evidence of rapid adaptive evolution in many different clades, including the carpet vipers, Echis,[44] the ground rattlesnakes, Sistrurus,[41] and the Malayan pit viper,[45] as well as generally in the diversification of PLA2 proteins.[34] There is phylogenetic evidence of positive selection and rapid rates of gene gain and loss in venom genes of Sistrurus taxa feeding on different prey.[46] As of 2019, evidence existed both of "overkill" occurring in some lineages, and rapid adaptive evolution, and an evolutionary arms race with prey physiology, in many others.[47]

The genes that code for venom proteins in some snake genera have a proportion of synonymous mutations that is lower than would be expected if venom were evolving through neutral evolutionary processes; the non-synonymous mutation rate, however, was found higher in many cases, indicating directional selection. [34] In addition, snake venom is metabolically costly for a snake to produce, which scientists have suggested as further evidence that a selection pressure exists on snake venom (in this case, to minimize the volume of venom required).[43] The use of model organisms, rather than snakes' natural prey, to study prey toxicity, has been suggested as a reason why the "overkill" hypothesis may have been overemphasized.[48] However, the pitviper genus Agkistrodon has been found to be an exception to this; the composition of venom in Agkistrodon has been found to be related to the position of the species within the phylogeny, suggesting that those venoms have evolved mostly through neutral processes (mutation and genetic drift), and that there may be significant variation in the selection pressure upon various snake venoms.[49]

Several studies have found evidence that venom and resistance to venom in prey species have evolved in a coevolutionary arms race.[43] For example, wood rats of the genus Neotoma have a high degree of resistance to the venom of rattlesnakes, suggesting that the rats have evolved in response to the snake venom, thus renewing selection pressure upon the snakes.[41] Resistance to venoms of sympatric predatory snake species has been found in eels,[50] ground squirrels,[51] rock squirrels,[52] and Eastern gray squirrels.[53] All these studies suggested a co-evolutionary arms race between prey and predator, indicating another potential selection pressure on snake venom to increase or innovate toxicity. The selection pressure on snake venom is thought to be selecting for functional diversity within the proteins in venom, both within a given species, and across species.[54] In addition to prey physiology, evidence exists that snake venom has evolved in response to the physiology of predators.[13]

Besides diet, there are other possible pressures on snake venom composition. A 2019 study found that larger body mass and smaller ecological habitats were correlated with increased venom yield.[47] Another study found that weather and temperature had stronger correlations with snake venom than diets or types of prey.[55] While venomous snakes use their venoms in defence (hence the problem of snakebite in humans), it is not well known to what extent natural selection for defence has driven venom evolution. The venoms of the Texas coral snake, Micrurus tener, and other species of Micrurus have been found to contain toxins with specific pain-inducing activity, suggesting a defensive function.[56] However, a questionnaire survey of snakebite patients bitten by a wide variety of venomous species showed that pain after most snakebites is of slow onset, arguing against widespread selection for defence.[3] The spitting of venom displayed by some species of spitting cobra is solely a defensive adaptation. A 2021 study showed that the venoms of all three lineages of spitting cobra convergently evolved higher levels of sensory neuron activation (i.e., cause more pain) than the venoms of non-spitting cobras, through the synergistic action of cytotoxins and Phospholipase A2 toxins, indicating selection for a defensive function.[57]

In contrast to venom composition and toxicity to specific lineages, venom yield, or the quantity of venom produced by an individual snake, has not been found to vary with the body-mass of prey animals, and instead to vary with the body-mass of snakes producing it. Yield increases with snake body-mass in a consistent with the hypothesis that snakes invest a constant proportion of metabolic output into producing venom, which is metabolically costly.[47]

Functional adaptations

edit
 
The Mangrove snake Boiga dendrophila, whose venom is toxic to birds.

Snakes use their venom to kill or subdue prey, as well as for other diet-related functions, such as digestion. Current scientific theory suggests that snake venom is not used for defense or for competition between members of the same species, unlike in other taxa. Thus adaptive evolution in snake venom has resulted in several adaptations with respect to these diet-related functions that increase the fitness of the snakes that carry them.[58][43][17] This is also reflected in variation in venom composition within a species; venom is known to vary geographically, and by age and size, likely reflecting variation in the prey consumed by different groups within a species.[13] Geographic variation is also present at the species level; island snakes tend to have less complex venoms, while those living in highly productive habitats have more complex venoms, suggesting a biogeographic pattern.[59]

Prey-specific venom toxicity

edit
 
Echis carinatus, one of the many species of the widespread genus Echis. The toxicity of Echis venom to scorpions has been found to vary with the proportion of arthropods in the snake's diet.

Venom that is toxic only to a certain taxon, or strongly toxic only to a certain taxon, has been found in a number of snakes, suggesting that these venoms have evolved via natural selection to subdue preferred prey species. Examples of this phenomenon have been found in the Mangrove snake Boiga dendrophila, which has a venom specifically toxic to birds,[48] as well as in the genera Echis and Sistrurus, and in sea snakes. The venom of Spilotes sulphureus which has two components, one of which is toxic to lizards but non-toxic in mammals, while the other is toxic in mammals and non-toxic in lizards.[60] However, while several snakes possess venom that is highly toxic to their preferred prey species, the reverse correlation is not necessarily true: the venoms of several snakes are toxic to taxa that they do not consume in high proportions. Most snake venom, for instance, is highly toxic to lizards, although the proportion of lizard prey varies among snake species. This has led researchers to suggest that toxicity to a certain taxon is nearly independent of toxicity to another distantly related taxon.[41]

 
Hydrophis cyanocinctus, a member of a clade that has a greatly simplified venom that evolved in response to a diet of fish

The natural diets of snakes in the widespread viper genus Echis are highly varied, and include arthropods, such as scorpions, as well as vertebrates. Various Echis species consume different quantities of arthropods in their diet.[61] A 2009 study injected scorpions with the venom of various Echis species, and found a high correlation between the proportion of arthropods that the snakes consumed in their natural habitat, and the toxicity of their venom to scorpions.[44] The researchers also found evidence that the evolution of venom more toxic to arthropods was related to an increase in the proportion of arthropods in the snakes' diet, and that diet and venom may have evolved by a process of coevolution. A phylogeny of the genus constructed using mitochondrial DNA showed that one instance of a change in venom composition in the species ancestral to all Echis snakes was correlated with a shift to an arthropod based diet, whereas another shift in a more recent lineage was correlated with a shift to a diet of vertebrates.[62] Despite the higher toxicity of the venom of arthropod-consuming species, it was not found to incapacitate or kill prey any faster than that of species with fewer arthropods in their diet. Thus, the venom is thought to have evolved to minimize the volume required, as the production of venom carries a significant metabolic cost, thus providing a fitness benefit.[62] This pattern is also found in other lineages.[63] Similar results were obtained by a 2012 study which found that the venom of arthropod-consuming Echis species was more toxic to locusts than that of vertebrate-consuming species.[64]

A 2009 study of the venom of four Sistrurus pit viper species found significant variation in the toxicity to mice.[41] This variation was related to the proportion of small mammals in the diet of those species.[41] The idea that Sistrurus venom had evolved to accommodate a mammal-based diet was supported by phylogenetic analysis. The researchers suggested that the basis for the difference in toxicity was the difference in muscle physiology in the various prey animals.[41] Two lineages of elapid snakes, common sea snakes and Laticauda sea kraits, have independently colonized marine environments, and shifted to a very simple diet based on teleosts, or ray-finned fish.[17] A 2005 study found that both these lineages have a much simpler set of venom proteins than their terrestrial relatives on the Australian continent, which have a more varied and complex diet.[17] These findings were confirmed by a 2012 study, which compared the venoms of Toxicocalamus longissimus, a terrestrial species, and Hydrophis cyanocinctus, a marine species, both within the subfamily Hydrophiinae (which is also within the Elapid family). Despite being closely related to one another, the marine species had a significantly simpler set of venom proteins.[16] The venoms of the sea snakes are nonetheless among the most toxic venoms known. It has been argued that since sea snakes are typically unable to prevent the escape of bitten prey, their venoms have evolved to act very rapidly.[65]

Pre-digestion of prey

edit
The venom of the prairie rattlesnake, Crotalus viridis (left) includes metalloproteinases (example on the right) which help digest the prey before the snake eats it.

The various subspecies of the rattlesnake genus Crotalus, produce venoms that carry out two conflicting functions. The venom immobilizes prey after a bite, and also helps digestion by breaking down tissues before the snake eats its prey.[18] As with other members of the family Viperidae, the venoms of Crotalus disrupt the homeostatic processes of prey animals. However, there is a wide variety of venom compositions among the species of Crotalus. A 2010 study found a 100-fold difference in the amount of metalloproteinase activity among the various snakes, with Crotalus cerberus having the highest activity and Crotalus oreganus concolor having the lowest. There was also a 15-fold variation in the amount of protease activity, with C. o. concolor and C. cerberus having the highest and lowest activities, respectively.[18]

Metalloproteinase activity causes hemorrhage and necrosis following a snake bite, a process which aids digestion. The activity of proteases, on the other hand, disrupts platelet and muscle function and damages cell membranes, and thus contributes to a quick death for the prey animal.[18] The study found that the venoms of Crotalus fell into two categories; those that favored metalloproteinases (Type I) and those that favored proteases (Type II). The study stated that these functions were essentially mutually exclusive; venoms had been selected for based on either their toxicity or their tenderizing potential. The researchers also hypothesized that the reason for this dichotomy was that a venom high in neurotoxicity, such as a type II venom, kills an animal quickly, preventing the relatively slower acting metalloproteinase from digesting tissue.[18]

Tracking bitten prey

edit
The western diamondback rattlesnake, Crotalus atrox (left), whose venom contains disintegrins (right) which allow it to track bitten prey

Another example of an adaptive function other than prey immobilization is the role of viperid venom in allowing the snake to track a prey animal it has bitten, a process known as "prey relocalization." This important adaptation allowed rattlesnakes to evolve the strike-and-release bite mechanism, which provided a huge benefit to snakes by minimizing contact with potentially dangerous prey animals.[19] However, this adaptation then requires the snake to track down the bitten animal in order to eat it, in an environment full of other animals of the same species. A 2013 study found that western diamondback rattlesnakes (Crotalus atrox) responded more actively to mouse carcasses that had been injected with crude rattlesnake venom. When the various components of the venom were separated out, the snakes responded to mice injected with two kinds of disintegrins. The study concluded that these disintegrin proteins were responsible for allowing the snakes to track their prey, by changing the odor of the bitten animal.[19]

Diet-based atrophication

edit
 
A shift to a diet of eggs has resulted in an atrophied venom system in the common egg-eater Dasypeltis scabra

Venom in a number of lineages of snakes is thought to have atrophied in response to dietary shifts.[7] A 2005 study in the marbled sea snake, Aipysurus eydouxii found that the gene for a three-fingered protein found in the venom had undergone a deletion of two nucleotide bases which made the venom 50–100 times less toxic than it had been previously. This change was correlated with a change in diet from fish to a diet consisting almost entirely of fish eggs, suggesting that the adaptation to an egg diet had removed the selection pressure needed to maintain a highly toxic venom, allowing the venom genes to accumulate deleterious mutations.[17][66] A similar venom degradation following a shift to an egg-based diet has been found in the Common egg-eater Dasypeltis scabra, whose diet consists entirely of birds' eggs, meaning that the snake had no use for its venom. This has led biologists to suggest that if venom is not used by a species, it is rapidly lost.[67]

References

edit

Citations

edit
  1. ^ Hargreaves et al. (a) 2014.
  2. ^ Casewell et al. 2013, pp. 218–220.
  3. ^ a b Ward-Smith et al. 2020.
  4. ^ Fry et al. 2012a, pp. 441–442.
  5. ^ a b Wuster et al. 2008.
  6. ^ Lomonte et al. (a) 2014, p. 326.
  7. ^ a b c d e f Fry et al. 2012a, pp. 434–436.
  8. ^ a b Fry et al. 2012a, pp. 424–436.
  9. ^ a b Casewell et al. 2013, pp. 224–227.
  10. ^ a b c d e f Reyes-Velasco et al. 2015.
  11. ^ a b c Hargreaves et al. (b) 2014, pp. 153–155.
  12. ^ a b Xie et al. 2022.
  13. ^ a b c d e Casewell et al. 2020, pp. 570–581.
  14. ^ Casewell et al. 2013, pp. 222–223.
  15. ^ Barlow et al. 2009, pp. 2447–2448.
  16. ^ a b Calvete et al. 2012, pp. 4094–4098.
  17. ^ a b c d e f Li et al. 2005.
  18. ^ a b c d e Mackessy 2010.
  19. ^ a b c Saviola et al. 2013.
  20. ^ a b Fry et al. 2012a, p. 443.
  21. ^ a b Fry et al. 2012a.
  22. ^ a b Lomonte et al. (a) 2014, pp. 326–327.
  23. ^ Mackessy 2010, p. 1464.
  24. ^ Casewell et al. 2013, pp. 225–227.
  25. ^ Palci et al. 2021.
  26. ^ Fry et al. 2012a, p. 435.
  27. ^ a b Almeida et al. 2021.
  28. ^ a b Sunagar & Abraham 2021.
  29. ^ Fry & Wuster 2004, p. 870.
  30. ^ Mount & Brown 2022, pp. 973–985.
  31. ^ Hargreaves et al. (b) 2014.
  32. ^ Rao et al. 2022.
  33. ^ Casewell et al. 2013, p. 223.
  34. ^ a b c Lynch 2007.
  35. ^ a b Dowell et al. 2016.
  36. ^ a b Casewell et al. 2013, p. 223–224.
  37. ^ Casewell et al. 2011.
  38. ^ Fry & Wuster 2004, p. 871.
  39. ^ Fry et al. 2012b.
  40. ^ a b c Mikheyev & Barua 2019.
  41. ^ a b c d e f g Gibbs & Mackessy 2009.
  42. ^ Lomonte et al. (a) 2014, p. 334.
  43. ^ a b c d Barlow et al. 2009, p. 2443.
  44. ^ a b Barlow et al. 2009, p. 2447.
  45. ^ Casewell et al. 2013, p. 220.
  46. ^ Gibbs & Rossiter 2008.
  47. ^ a b c Healy, Carbone & Jackson, pp. 527–537.
  48. ^ a b Casewell et al. 2013, pp. 220–221.
  49. ^ Lomonte et al. (b) 2014, pp. 112–114.
  50. ^ Heatwole & Poran 1995.
  51. ^ Biardi, Chien & Coss 2005.
  52. ^ Biardi & Coss 2011.
  53. ^ Pomento et al. 2016.
  54. ^ Sanz et al. 2006, pp. 2098–2099.
  55. ^ Zancolli et al. 2019.
  56. ^ Bohlen 2011.
  57. ^ Kazandjian et al. 2021.
  58. ^ Casewell et al. 2013, pp. 219–220.
  59. ^ Siqueira‐Silva et al. 2021, pp. 1978–1989.
  60. ^ Modahl et al. 2018.
  61. ^ Barlow et al. 2009, pp. 2444, 2447.
  62. ^ a b Barlow et al. 2009, pp. 2446–2448.
  63. ^ Casewell et al. 2013, pp. 223–225.
  64. ^ Richards et al. 2012.
  65. ^ Calvete et al. 2012, pp. 4092–4093.
  66. ^ Calvete et al. 2012, pp. 4097–4098.
  67. ^ Fry et al. 2008.

Cited sources

edit
edit