Gene Conversion Facilitates the Adaptive Evolution of Self-Resistance in Highly Toxic Newts

Reconstructing the histories of complex adaptations and identifying the evolutionary mechanisms underlying their origins are two of the primary goals of evolutionary biology. Taricha newts, which contain high concentrations of the deadly toxin tetrodotoxin (TTX) as an antipredator defense, have evolved resistance to self-intoxication, which is a complex adaptation requiring changes in six paralogs of the voltage-gated sodium channel (Na_v) gene family, the physiological target of TTX. Here, we reconstruct the origins of TTX self-resistance by sequencing the entire Na_v gene family in newts and related salamanders. We show that moderate TTX resistance evolved early in the salamander lineage in three of the six Na_v paralogs, preceding the proposed appearance of tetrodotoxic newts by ∼100 My. TTX-bearing newts possess additional unique substitutions across the entire Na_v gene family that provide physiological TTX resistance. These substitutions coincide with signatures of positive selection and relaxed purifying selection, as well as gene conversion events, that together likely facilitated their evolution. We also identify a novel exon duplication within Na_v1.4 encoding an expressed TTX-binding site. Two resistance-conferring changes within newts appear to have spread via nonallelic gene conversion: in one case, one codon was copied between paralogs, and in the second, multiple substitutions were homogenized between the duplicate exons of Na_v1.4. Our results demonstrate that gene conversion can accelerate the coordinated evolution of gene families in response to a common selection pressure.     

Is there more than one way to skin a newt? Convergent toxin resistance in snakes is not due to a common genetic mechanism

Convergent evolution of tetrodotoxin (TTX) resistance, at both the phenotypic and genetic levels, characterizes coevolutionary arms races between amphibians and their snake predators around the world, and reveals remarkable predictability in the process of adaptation. Here we examine the repeatability of the evolution of TTX resistance in an undescribed predator–prey relationship between TTX-bearing Eastern Newts (Notophthalmus viridescens) and Eastern Hog-nosed Snakes (Heterodon platirhinos). We found that that local newts contain levels of TTX dangerous enough to dissuade most predators, and that Eastern Hog-nosed Snakes within newt range are highly resistant to TTX. In fact, these populations of Eastern Hog-nosed Snakes are so resistant to TTX that the potential for current reciprocal selection might be limited. Unlike all other cases of TTX resistance in vertebrates, H. platirhinos lacks the adaptive amino acid substitutions in the skeletal muscle sodium channel that reduce TTX binding, suggesting that physiological resistance in Eastern Hog-nosed Snakes is conferred by an alternate genetic mechanism. Thus, phenotypic convergence in this case is not due to parallel molecular evolution, indicating that there may be more than one way for this adaptation to arise, even among closely related species.     

河豚毒素生态作用研究进展

河豚毒素(tetrodotoxin, TTX)取名于河豚鱼,最早从河豚鱼中分离纯化.自1964年河豚毒素化学结构被阐明以后,河豚毒素研究得到了生物学家、毒理学家、化学家、药理学家的广泛关注.河豚毒素具有许多天然同系物.河豚毒素及其同系物在自然界分布广泛,存在于一系列不同进化水平的海洋生物和少量的两栖动物体内.河豚毒素及其同系物可能具有防御、捕食及信息传递等生态作用,毒素在生物体内的分布与其生态作用密切相关.含有河豚毒素的生物对河豚毒素具有一定的耐受能力,其机制可能与生物体内存在河豚毒素结合蛋白或生物自身具有独特的钠离子通道结构有关.重点针对河豚毒素的生态作用及生物对河豚毒素的耐受机制进行了综述,以期为河豚毒素生态学研究及河豚毒素中毒事件的防范提供科学资料.     

Convergent adaptation to dangerous prey proceeds through the same first‐step mutation in the garter snake Thamnophis sirtalis

Convergent phenotypes often result from similar underlying genetics, but recent work suggests convergence may also occur in the historical order of substitutions en route to an adaptive outcome. We characterized convergence in the mutational steps to two independent outcomes of tetrodotoxin (TTX) resistance in separate geographic lineages of the common garter snake (Thamnophis sirtalis) that coevolved with toxic newts. Resistance is largely conferred by amino acid changes in the skeletal muscle sodium channel (Na_V1.4) that interfere with TTX‐binding. We sampled variation in Na_V1.4 throughout western North America and found clear evidence that TTX‐resistant changes in both lineages began with the same isoleucine‐valine mutation (I1561V) within the outer pore of Na_V1.4. Other point mutations in the pore, shown to confer much greater resistance, accumulate later in the evolutionary progression and always occur together with the initial I1561V change. A gene tree of Na_V1.4 suggests the I1561V mutations in each lineage are not identical‐by‐decent, but rather they arose independently. Convergence in the evolution of channel resistance is likely the result of shared biases in the two lineages of T. sirtalis—only a few mutational routes can confer TTX resistance while maintaining the conserved function of voltage‐gated sodium channels.     

Morphological variation of hypaxial musculature in salamanders (Lissamphibia: Caudata)

Despite the acknowledged importance of the locomotory and respiratory functions associated with hypaxial musculature in salamanders, variation in gross morphology of this musculature has not been documented or evaluated within a phylogenetic or ecological context. In this study, we characterize and quantify the morphological variation of lateral hypaxial muscles using phylogenetically and ecologically diverse salamander species from eight families: Ambystomatidae (Ambystoma tigrinum), Amphiumidae (Amphiuma tridactylum), Cryptobranchidae (Cryptobranchus alleganiensis), Dicamptodontidae (Dicamptodon sp.), Plethodontidae (Gyrinophilus porphyriticus), Proteidae (Necturus maculosus), Salamandridae (Pachytriton sp.), and Sirenidae (Siren lacertina). For the lateral hypaxial musculature, we document 1) the presence or absence of muscle layers, 2) the muscle fiber angles of layers at mid‐trunk, and 3) the relative dorsoventral positions and cross‐sectional areas of muscle layers. Combinations of two, three, or four layers are observed. However, all species retain at least two layers with opposing fiber angles. The number of layers and the presence or absence of layers vary within species (Necturus maculosus and Siren lacertina), within genera (e.g., Triturus), and within families. No phylogenetic pattern in the number of layers can be detected with a family‐level phylogeny. Fiber angle variation of hypaxial muscles is considerable: fiber angles of the M. obliquus externus range from 20–80°; M. obliquus internus, 14–34°; M. transversus abdominis, 58–80° (acute angles measured relative to the horizontal septum). Hypaxial musculature comprises 17–37% of total trunk cross‐sectional area. Aquatic salamanders show relatively larger total cross‐sectional hypaxial area than salamanders that are primarily terrestrial. J. Morphol. 241:153–164, 1999. © 1999 Wiley‐Liss, Inc.     

Genetic architecture of a feeding adaptation: garter snake (Thamnophis) resistance to tetrodotoxin bearing prey

Detailing the genetic basis of adaptive variation in natural populations is a first step towards understanding the process of adaptive evolution, yet few ecologically relevant traits have been characterized at the genetic level in wild populations. Traits that mediate coevolutionary interactions between species are ideal for studying adaptation because of the intensity of selection and the well-characterized ecological context. We have previously described the ecological context, evolutionary history and partial genetic basis of tetrodotoxin (TTX) resistance in garter snakes (Thamnophis). Derived mutations in a voltage-gated sodium channel gene (Na_v1.4) in three garter snake species are associated with resistance to TTX, the lethal neurotoxin found in their newt prey (Taricha). Here we evaluate the contribution of Na_v1.4 alleles to TTX resistance in two of those species from central coastal California. We measured the phenotypes (TTX resistance) and genotypes (Na_v1.4 and microsatellites) in a local sample of Thamnophis atratus and Thamnophis sirtalis. Allelic variation in Na_v1.4 explains 23 per cent of the variation in TTX resistance in T. atratus while variation in a haphazard sample of the genome (neutral microsatellite markers) shows no association with the phenotype. Similarly, allelic variation in Na_v1.4 correlates almost perfectly with TTX resistance in T. sirtalis, but neutral variation does not. These strong correlations suggest that Na_v1.4 is a major effect locus. The simple genetic architecture of TTX resistance in garter snakes may significantly impact the dynamics of phenotypic coevolution. Fixation of a few alleles of major effect in some garter snake populations may have led to the evolution of extreme phenotypes and an ‘escape’ from the arms race with newts.     

TTX accumulation in pufferfish

Tetrodotoxin (TTX) has been detected in a variety of animals. The finding of TTX in the trumpet shell Charonia sauliae strongly suggested that its origin was its food, a TTX-bearing starfish Astropecten polyacanthus. Since then, the food chain has been consistently implicated as the principal means of TTX intoxication. To identify the primary producer of TTX, intestinal bacteria isolated from several TTX-bearers were investigated for their TTX production. The results demonstrated that some of them could produce TTX. Thus the primary TTX producers in the sea are concluded to be marine bacteria. Subsequently, detritus feeders and zooplankton can be intoxicated with TTX through the food chain, or in conjunction with parasitism or symbiosis. The process followed by small carnivores, omnivores or scavengers, and by organisms higher up the food chain would result in the accumulation of higher concentrations of TTX. Finally, pufferfish at the top of the food chain are intoxicated with TTX. This hypothesis is supported by the fact that net cage and land cultures produce non-toxic pufferfish that can be made toxic by feeding with a TTX-containing diet.     

EVOLUTIONARY RESPONSE OF PREDATORS TO DANGEROUS PREY: PREADAPTATION AND THE EVOLUTION OF TETRODOTOXIN RESISTANCE IN GARTER SNAKES

Coevolutionary interactions typically involve only a few specialized taxa. The factors that cause some taxa and not others to respond evolutionarily to selection by another species are poorly understood. Preadaptation may render some species predisposed for evolutionary response to new pressures, whereas a lack of genetic variation may limit the evolutionary potential of other taxa. We evaluate these factors in the predator-prey interaction between toxic newts (Taricha granulosa) and their resistant garter snake predators (Thamnophis sirtalis). Using a bioassay of resistance to tetrodotoxin (TTX), the primary toxin in the prey, we examined phenotypic evolution in the genus Thamnophis. Reconstruction of ancestral character states suggests that the entire genus Thamnophis, and possibly natricine snakes in general, has slightly elevated TTX resistance compared to other lineages of snakes. While this suggests that T. sirtalis is indeed predisposed to evolving TTX resistance, it also indicates that the potential exists in sympatric congeners not expressing elevated levels of TTX resistance. We also detected significant family level variation for TTX resistance in a species of Thamnophis that does not exhibit elaborated levels of the trait. This finding suggests that evolutionary response in other taxa is not limited by genetic variability. In this predator-prey system, species and population differences in resistance appear to be largely determined by variation in the selective environment rather than preadaptation or constraint.     

Genetic basis of tetrodotoxin resistance in pufferfishes

Tetrodotoxin (TTX) is a highly potent neurotoxin that selectively binds to the outer vestibule of voltage-gated sodium channels. Pufferfishes accumulate extremely high concentrations of TTX without any adverse effect. A nonaromatic amino acid (Asn) residue present in domain I of the pufferfish, Takifugu pardalis, Na v1.4 channel has been implicated in the TTX resistance of pufferfishes . However, the effect of this residue on TTX sensitivity has not been investigated, and it is not known if this residue is conserved in all pufferfishes. We have investigated the genetic basis of TTX resistance in pufferfishes by comparing the sodium channels from two pufferfishes (Takifugu rubripes [fugu] and Tetraodon nigroviridis) and the TTX-sensitive zebrafish. Although all three fishes contain duplicate copies of Na v1.4 channels (Na v1.4a and Na v1.4b), several substitutions were found in the TTX binding outer vestibule of the two pufferfish channels. Electrophysiological studies showed that the nonaromatic residue (Asn in fugu and Cys in Tetraodon) in domain I of Na v1.4a channels confers TTX resistance. The Glu-to-Asp mutation in domain II of Tetraodon channel Na v1.4b is similar to that in the saxitoxin- and TTX-resistant Na+ channels of softshell clams . Besides helping to deter predators, TTX resistance enables pufferfishes to selectively feed on TTX-bearing organisms.     

Toxicity of pufferfish Takifugu rubripes cultured in netcages at sea or aquaria on land

Marine pufferfish (family Tetraodontidae) are believed to accumulate tetrodotoxin (TTX) mainly in liver and ovary through the food chain by ingesting TTX-bearing organisms such as starfish, gastropods, crustacean, flatworms, ribbonworms, etc. Consequently, it is hypothesized that non-toxic pufferfish can be produced if they are cultured with TTX-free diets in netcages at sea or aquaria on land, where the invasion of TTX-bearing organisms is completely shut off. To confirm this hypothesis, more than 5000 specimens of the pufferfish (“torafugu”, Takifugu rubripes) cultured in such manners for 1–3 years were collected from several locations in Japan during 2001–2004, and toxicity of their livers and some other parts was examined according to the Japanese official mouse assay method for TTX. In addition, typical specimens were submitted to LC/MS analysis. The results showed that all the livers and other parts tested were ‘non-toxic’ in both of the mouse assay (less than 2 MU/g) and LC/MS analysis (less than 0.1 MU/g). Thus, it is undoubtedly confirmed that pufferfish are intoxicated through the food chain, and non-toxic pufferfish can be successfully produced by netcage or land culture. The livers from these fish can be used with safety as a Japanese traditional food “fugu-kimo” (puffer liver).     

A key to totipotency: Wuschel-like homeobox 2a unlocks embryogenic culture response in maize (Zea mays L.)

The ability of plant somatic cells to dedifferentiate, form somatic embryos and regenerate whole plants in vitro has been harnessed for both clonal propagation and as a key component of plant genetic engineering systems. Embryogenic culture response is significantly limited, however, by plant genotype in most species. This impedes advancements in both plant transformation-based functional genomics research and crop improvement efforts. We utilized natural variation among maize inbred lines to genetically map somatic embryo generation potential in tissue culture and identify candidate genes underlying totipotency. Using a series of maize lines derived from crosses involving the culturable parent A188 and the non-responsive parent B73, we identified a region on chromosome 3 associated with embryogenic culture response and focused on three candidate genes within the region based on genetic position and expression pattern. Two candidate genes showed no effect when ectopically expressed in B73, but the gene Wox2a was found to induce somatic embryogenesis and embryogenic callus proliferation. Transgenic B73 cells with strong constitutive expression of the B73 and A188 coding sequences of Wox2a were found to produce somatic embryos at similar frequencies, demonstrating that sufficient expression of either allele could rescue the embryogenic culture phenotype. Transgenic B73 plants were regenerated from the somatic embryos without chemical selection and no pleiotropic effects were observed in the Wox2a overexpression lines in the regenerated T0 plants or in the two independent events which produced T1 progeny. In addition to linking natural variation in tissue culture response to Wox2a, our data support the utility of Wox2a in enabling transformation of recalcitrant genotypes.     

Myostatin rapid sequence evolution in ruminants predates domestication

Myostatin (GDF-8) is a negative regulator of skeletal muscle development. This gene has previously been implicated in the double muscling phenotype in mice and cattle. A systematic analysis of myostatin sequence evolution in ruminants was performed in a phylogenetic context. The myostatin coding sequence was determined from duiker (Sylvicapra grimmia caffra), eland (Taurotragus derbianus), gaur (Bos gaurus), ibex (Capra ibex), impala (Aepyceros melampus rednilis), pronghorn (Antilocapra americana), and tahr (Hemitragus jemlahicus). Analysis of nonsynonymous to synonymous nucleotide substitution rate ratios (K_a/K_s) indicates that positive selection may have been operating on this gene during the time of divergence of Bovinae and Antilopinae, starting from approximately 23 million years ago, a period that appears to account for most of the sequence difference between myostatin in these groups. These periods of positive selective pressure on myostatin may correlate with changes in skeletal muscle mass during the same period.     

Isolation of a <Emphasis Type="Italic">Ve</Emphasis> homolog, <Emphasis Type="Italic">mVe1</Emphasis>, and its relationship to verticillium wilt resistance in <Emphasis Type="Italic">Mentha longifolia</Emphasis> (L.) Huds

As a step toward greater understanding of the genetics of verticillium wilt resistance in plants, we report the sequencing of a candidate wilt resistance gene, mVe1, from the mint diploid model species, Mentha longifolia (Lamiaceae). mVe1 is a putative homolog of tomato (Solanum lycopersicum L.) verticillium wilt (Ve) resistance genes. The mVe1 gene has a coding region of 3,051 bp. The predicted mVe1 protein contains a leucine-rich repeat domain, a common feature of plant disease resistance proteins. We compared 13 mVe1 alleles from three mint species. These alleles shared 96.2–99.6% nucleotide identity. We analyzed four M. longifolia populations segregating with respect to mVe1 alleles and wilt resistance versus susceptibility and found one association between mVe1 genotype and wilt phenotype. We conclude that mVe1 may play a role in mint verticillium wilt resistance, but variation for resistance in our segregating progenies is likely polygenic. Therefore, further investigations of mVe1 and identification of additional candidate genes are both warranted.     

Pronounced differences in genetic structure despite overall ecological similarity for two Ambystoma salamanders in the same landscape

Studies linking genetic structure in amphibian species with ecological characteristics have focused on large differences in dispersal capabilities. Here, we test whether two species with similar dispersal potential but subtle differences in other ecological characteristics also exhibit strong differences in genetic structure in the same landscape. We examined eight microsatellites in marbled salamanders (Ambystoma opacum) from 29 seasonal ponds and spotted salamanders (Ambystoma maculatum) from 19 seasonal ponds in a single geographic region in west-central Massachusetts. Despite overall similarity in ecological characteristics of spotted and marbled salamanders, we observed clear differences in the genetic structure of these two species. For marbled salamanders, we observed strong overall genetic differentiation (F _ST = 0.091, F′ _ST = 0.375), three population-level clusters of populations (K = 3), a strong pattern of isolation by distance (r = 0.58), and marked variation in family-level structure (from 1 to 23 full-sibling families per site). For spotted salamanders, overall genetic differentiation was weaker (F _ST = 0.025, F′ _ST = 0.102), there was no evidence of population-level clustering (K = 1), the pattern of isolation by distance (r = 0.17) was much weaker compared to marbled salamanders, and there was less variation in family-level structure (from 10 to 36 full-sibling families per site). We suspect that a combination of breeding site fidelity, effective population size, and generation interval is responsible for these marked differences. Our results suggest that marbled salamanders, compared to spotted salamanders, are more sensitive to fragmentation from various land-use activities and would be less likely to recolonize extirpated sites on an ecologically and conservation-relevant time frame.     

Toxicity and population structure of the Rough‐Skinned Newt (Taricha granulosa) outside the range of an arms race with resistant predators

Species interactions, and their fitness consequences, vary across the geographic range of a coevolutionary relationship. This spatial heterogeneity in reciprocal selection is predicted to generate a geographic mosaic of local adaptation, wherein coevolutionary traits are phenotypically variable from one location to the next. Under this framework, allopatric populations should lack variation in coevolutionary traits due to the absence of reciprocal selection. We examine phenotypic variation in tetrodotoxin (TTX) toxicity of the Rough‐Skinned Newt (Taricha granulosa) in regions of allopatry with its TTX‐resistant predator, the Common Garter Snake (Thamnophis sirtalis). In sympatry, geographic patterns of phenotypic exaggeration in toxicity and toxin‐resistance are closely correlated in prey and predator, implying that reciprocal selection drives phenotypic variation in coevolutionary traits. Therefore, in allopatry with TTX‐resistant predators, we expect to find uniformly low levels of newt toxicity. We characterized TTX toxicity in northwestern North America, including the Alaskan panhandle where Ta. granulosa occur in allopatry with Th. sirtalis. First, we used microsatellite markers to estimate population genetic structure and determine if any phenotypic variation in toxicity might be explained by historical divergence. We found northern populations of Ta. granulosa generally lacked population structure in a pattern consistent with northern range expansion after the Pleistocene. Next, we chose a cluster of sites in Alaska, which uniformly lacked genetic divergence, to test for phenotypic divergence in toxicity. As predicted, overall levels of newt toxicity were low; however, we also detected unexpected among‐ and within‐population variation in toxicity. Most notably, a small number of individuals contained large doses of TTX that rival means of toxic populations in sympatry with Th. sirtalis. Phenotypic variation in toxicity, despite limited neutral genetic divergence, suggests that factors other than reciprocal selection with Th. sirtalis likely contribute to geographic patterns of toxicity in Ta. granulosa.     

Pigment pattern formation in the larval salamander Ambystoma maculatum

As part of an ongoing comparative study of pigment patterns and their formation in embryos and larvae of ambystomatid salamanders, Ambystoma maculatum from two differnt populations, one in the northern (New York) and one in the central (Tennessee) United States, were investigated. Scanning electron microscopy was used to study early neural crest development. Light microscopy in combination with markers for the two pigment cell types (xanthophores and melanophores) made it possible to follow pigment cell migration before the pigment cells were fully differentiated. A bilateral pigment pattern consisting of two horizontal melanophore stripes surrounding an interstripe area populated by xanthophores formed in the larvae. In both populations, some variation was present in the form of a continuum ranging from clear horizontal stripes to extreme cases with a random pattern. Unlike the other ambystomatids that have been investigated, the neural crest cells in A. maculatum do not form aggregates and no vertical bars are formed. Instead, both the pattern and its formation are very similar to what has been reported for salamandrids. If pattern formation mechanisms can act as developmental constraints we would expect the A. maculatum pattern to be the primitive condition in the Ambystomatidae, using the Salamandridae as the outgroup. There is no strong support for this when aggregate formation is used as a character and mapped onto phylogenies for the group. The aggregate formation mechanism, and the pigment pattern that it leads to, have most likely been secondarily lost in A. maculatum. © 1993 Wiley-Liss, Inc.     

Larval growth in polyphenic salamanders: making the best of a bad lot

Polyphenisms are excellent models for studying phenotypic variation, yet few studies have focused on natural populations. Facultative paedomorphosis is a polyphenism in which salamanders either metamorphose or retain their larval morphology and eventually become paedomorphic. Paedomorphosis can result from selection for capitalizing on favorable aquatic habitats (paedomorph advantage), but could also be a default strategy under poor aquatic conditions (best of a bad lot). We tested these alternatives by quantifying how the developmental environment influences the ontogeny of wild Arizona tiger salamanders (Ambystoma tigrinum nebulosum). Most paedomorphs in our study population arose from slow-growing larvae that developed under high density and size-structured conditions (best of a bad lot), although a few faster-growing larvae also became paedomorphic (paedomorph advantage). Males were more likely to become paedomorphs than females and did so under a greater range of body sizes than females, signifying a critical role for gender in this polyphenism. Our results emphasize that the same phenotype can be adaptive under different environmental and genetic contexts and that studies of phenotypic variation should consider multiple mechanisms of morph production.