How Metamorphosis Is Different in Plethodontids: Larval Life History Perspectives on Life-Cycle Evolution

Plethodontid salamanders exhibit biphasic, larval form paedomorphic, and direct developing life cycles. This diversity of developmental strategies exceeds that of any other family of terrestrial vertebrate. Here we compare patterns of larval development among the three divergent lineages of biphasic plethodontids and other salamanders. We discuss how patterns of life-cycle evolution and larval ecology might have produced a wide array of larval life histories. Compared with many other salamanders, most larval plethodontids have relatively slow growth rates and sometimes exceptionally long larval periods (up to 60 mo). Recent phylogenetic analyses of life-cycle evolution indicate that ancestral plethodontids were likely direct developers. If true, then biphasic and paedomorphic lineages might have been independently derived through different developmental mechanisms. Furthermore, biphasic plethodontids largely colonized stream habitats, which tend to have lower productivity than seasonally ephemeral ponds. Consistent with this, plethodontid larvae grow very slowly, and metamorphic timing does not appear to be strongly affected by growth history. On the basis of this, we speculate that feeding schedules and stress hormones might play a comparatively reduced role in governing the timing of metamorphosis of stream-dwelling salamanders, particularly plethodontids.

Phylogeny of Mental Glands, Revisited

Mental glands and their associated delivery behaviors during courtship are unique to the plethodontid salamanders. Because previous interpretations of the evolution of these features were conducted using older phylogenetic hypotheses, we reanalyzed these traits with newer courtship descriptions and contemporary phylogenetic methods. Using Bayesian ancestral state reconstruction methods that have been developed since the first phylogenetic analyses were conducted in the mid-1990s, we reconstructed mental gland and courtship behavior evolution on a Bayesian phylogeny of the nuclear gene Rag1. The most probable ancestral condition for plethodontids was resolved as presence of a mental gland. Loss of a mental gland occurred in each subfamily and was recovered as the most probable ancestral condition for the Spelerpinae. In contrast, parsimony reconstruction recovered the presence of a mental gland in the ancestor to Spelerpinae with multiple secondary losses. We hypothesize that that absence of a mental gland is possibly ancestral in some clades (i.e., Spelerpinae) and secondary in others (e.g., paedomorphic Eurycea). The most probable ancestral form of the mental gland is likely to be the large pad-type distributed extensively in Plethodontinae and Bolitoglossinae. Desmognathans have the most unique mental glands, occurring in an anterior protrusion or bifurcated form (in Desmognathus wrighti). Fan-shaped mental glands evolved independently in Eurycea and Oedipina. Small pads arose independently in Bolitoglossinae, Plethodontinae, and Spelerpinae. Head-rubbing behavior for mental gland delivery mode was recovered as the most probable and parsimonious ancestral state for the Plethodontidae, with independent losses of this behavior in Plethodontinae and Spelerpinae. Because head-rubbing was observed in outgroups, we hypothesize that head-rubbing behavior predated mental gland evolution. Pulling, snapping, slapping, and biting behaviors evolved independently in the Plethodontinae and Spelerpinae and are not homologous with head-rubbing. All hypotheses of mental gland and courtship evolution invoke homoplasy.

Experimental Analysis of Allocation during Larval Development in Ambystomatid Salamanders

We examined how variation in food availability and temperature influences patterns of energy allocation to growth, storage, metamorphosis, and reproduction in larval Ambystoma maculatum and A. tigrinum. In each species, an allocation vector consisting of metamorphic size, duration of larval period, fat body mass, and gonad mass was influenced by both food and temperature. Larvae grew rapidly at high food levels and high temperature, and delayed metamorphosis when grown at low temperature. High food levels resulted in larger fat bodies and decreased time to metamorphosis in A. tigrinum, and larger gonads in both species. The allocation vector collectively representing growth rate, storage, reproduction and metamorphosis differed between species: in A. maculatum, only the correlation between metamorphic size and fat body mass was significant, whereas all four variables were influenced in A. tigrinum. An interaction of food and temperature on the allocation vector was only detected for A. tigrinum. Our results indicate that interaction among allocation variables is stronger in some species than others and that these associations might be important in understanding life history variation between these species.

Transcriptional response of Mexican axolotls to Ambystoma tigrinum virus (ATV) infection

Background

Very little is known about the immunological responses of amphibians to pathogens that are causing global population declines. We used a custom microarray gene chip to characterize gene expression responses of axolotls (Ambystoma mexicanum) to an emerging viral pathogen, Ambystoma tigrinum virus (ATV).

Result

At 0, 24, 72, and 144 hours post-infection, spleen and lung samples were removed for estimation of host mRNA abundance and viral load. A total of 158 up-regulated and 105 down-regulated genes were identified across all time points using statistical and fold level criteria. The presumptive functions of these genes suggest a robust innate immune and antiviral gene expression response is initiated by A. mexicanum as early as 24 hours after ATV infection. At 24 hours, we observed transcript abundance changes for genes that are associated with phagocytosis and cytokine signaling, complement, and other general immune and defense responses. By 144 hours, we observed gene expression changes indicating host-mediated cell death, inflammation, and cytotoxicity.

Conclusion

Although A. mexicanum appears to mount a robust innate immune response, we did not observe gene expression changes indicative of lymphocyte proliferation in the spleen, which is associated with clearance of Frog 3 iridovirus in adult Xenopus. We speculate that ATV may be especially lethal to A. mexicanum and related tiger salamanders because they lack proliferative lymphocyte responses that are needed to clear highly virulent iridoviruses. Genes identified from this study provide important new resources to investigate ATV disease pathology and host-pathogen dynamics in natural populations.

Effect of thyroid hormone concentration on the transcriptional response underlying induced metamorphosis in the Mexican axolotl (Ambystoma)

Background

Thyroid hormones (TH) induce gene expression programs that orchestrate amphibian metamorphosis. In contrast to anurans, many salamanders do not undergo metamorphosis in nature. However, they can be induced to undergo metamorphosis via exposure to thyroxine (T₄). We induced metamorphosis in juvenile Mexican axolotls (Ambystoma mexicanum) using 5 and 50 nM T₄, collected epidermal tissue from the head at four time points (Days 0, 2, 12, 28), and used microarray analysis to quantify mRNA abundances.

Results

Individuals reared in the higher T₄ concentration initiated morphological and transcriptional changes earlier and completed metamorphosis by Day 28. In contrast, initiation of metamorphosis was delayed in the lower T₄ concentration and none of the individuals completed metamorphosis by Day 28. We identified 402 genes that were statistically differentially expressed by ≥ two-fold between T₄ treatments at one or more non-Day 0 sampling times. To complement this analysis, we used linear and quadratic regression to identify 542 and 709 genes that were differentially expressed by ≥ two-fold in the 5 and 50 nM T₄ treatments, respectively.

Conclusion

We found that T₄ concentration affected the timing of gene expression and the shape of temporal gene expression profiles. However, essentially all of the identified genes were similarly affected by 5 and 50 nM T₄. We discuss genes and biological processes that appear to be common to salamander and anuran metamorphosis, and also highlight clear transcriptional differences. Our results show that gene expression in axolotls is diverse and precise, and that axolotls provide new insights about amphibian metamorphosis.

Microarray analysis identifies keratin loci as sensitive biomarkers for thyroid hormone disruption in the salamander Ambystoma mexicanum

Ambystomatid salamanders offer several advantages for endocrine disruption research, including genomic and bioinformatics resources, an accessible laboratory model (Ambystoma mexicanum), and natural lineages that are broadly distributed among North American habitats. We used microarray analysis to measure the relative abundance of transcripts isolated from A. mexicanum epidermis (skin) after exogenous application of thyroid hormone (TH). Only one gene had a >2-fold change in transcript abundance after 2 days of TH treatment. However, hundreds of genes showed significantly different transcript levels at days 12 and 28 in comparison to day 0. A list of 123 TH-responsive genes was identified using statistical, BLAST, and fold level criteria. Cluster analysis identified two groups of genes with similar transcription patterns: up-regulated versus down-regulated. Most notably, several keratins exhibited dramatic (1000 fold) increases or decreases in transcript abundance. Keratin gene expression changes coincided with morphological remodeling of epithelial tissues. This suggests that keratin loci can be developed as sensitive biomarkers to assay temporal disruptions of larval-to-adult gene expression programs. Our study has identified the first collection of loci that are regulated during TH-induced metamorphosis in a salamander, thus setting the stage for future investigations of TH disruption in the Mexican axolotl and other salamanders of the genus Ambystoma.

Early gene expression during natural spinal cord regeneration in the salamander Ambystoma mexicanum

In contrast to mammals, salamanders have a remarkable ability to regenerate their spinal cord and recover full movement and function after tail amputation. To identify genes that may be associated with this greater regenerative ability, we designed an oligonucleotide microarray and profiled early gene expression during natural spinal cord regeneration in Ambystoma mexicanum. We sampled tissue at five early time points after tail amputation and identified genes that registered significant changes in mRNA abundance during the first 7 days of regeneration. A list of 1036 statistically significant genes was identified. Additional statistical and fold change criteria were applied to identify a smaller list of 360 genes that were used to describe predominant expression patterns and gene functions. Our results show that a diverse injury response is activated in concert with extracellular matrix remodeling mechanisms during the early acute phase of natural spinal cord regeneration. We also report gene expression similarities and differences between our study and studies that have profiled gene expression after spinal cord injury in rat. Our study illustrates the utility of a salamander model for identifying genes and gene functions that may enhance regenerative ability in mammals.

Effects of lung reduction on rheotactic performance in amphibian larvae

In this study we test the hypothesis that lung reduction in amphibians is a larval adaptation to reduce buoyancy and facilitate rheotaxis in lotic environments. Larvae of the spotted salamander Ambystoma maculatum (family Ambystomatidae) were raised from eggs with (= control) and without (= experimental) access to the air-water interface. The experimental larvae developed significantly smaller and more solidified lungs than the controls. Although mean snout-vent length did not differ between the control and experimental groups, the latter had a significantly greater mean body mass, indicating greater specific gravity in experimental versus controls. To evaluate rheotactic behavior, larvae were subjected to a water current in an artificial stream. The distance and time required to establish a purchase on the stream bottom were recorded. Although the times did not differ, the distances displaced downstream were significantly less in experimental than in control larvae. We conclude that lung reduction and solidification reduce buoyancy and enhance rheotactic ability in salamander larvae. These findings suggest that lunglessness in plethodontid salamanders evolved as an adaptation to stream environments.

Effects of developmental and growth history on metamorphosis in the gray treefrog, Hyla versicolor (Amphibia, Anura)

In ecological models, the timing of amphibian metamorphosis is dependent upon rate of larval growth, e.g., tadpoles that experience a decrease in growth rate can initiate metamorphosis early. Recent authors have suggested that this plasticity may be lost at some point during the larval period. We tested this hypothesis by exposing groups of tadpoles of the gray treefrog, Hyla versicolor, to different growth schedules. In endocrine models, metamorphosis is dependent on thyroxine levels and thyroxine is antagonized by prolactin (amphibian larval growth hormone), consistent with the idea that a rapidly growing tadpole can delay metamorphosis. Thus, we also manipulated the rate of development by supplementing or maintaining natural thyroxine levels for half of the tadpoles in each growth treatment. All tadpoles that received thyroxine supplements metamorphosed at the same time regardless of growth history. They also metamorphosed earlier than tadpoles not treated with thyroxine. Tadpoles not given thyroxine supplements metamorphosed at different times: those growing rapidly during day 15-34 metamorphosed earlier than tadpoles growing slowly. Growth rate before day 15 and after day 34 had no effect on metamorphic timing. The difference in larval period between these rapidly growing tadpoles and their sisters given thyroxine treatments was less than the same comparison for tadpoles that grew slowly during the same period. This apparent prolactin/thyroxine antagonism did not exist after day 34. These results are consistent with the hypothesis of a loss of plasticity in metamorphic timing.