Molecular systematics and phylogeography of the Plethodon elongatus species group: combining phylogenetic and population genetic methods to investigate species history

Plethodon elongatus and P. stormi (Caudata: Plethodontidae) are Pacific Northwest endemic species which occur in northwestern California and southwestern Oregon. Studies on these salamanders have resulted in differing taxonomic conclusions, but the underlying historical hypotheses, at both inter- and intraspecific levels, have never been examined in a molecular framework. Here, representatives of 81 populations from throughout the range of both taxa are sequenced. Portions of three mitochondrial protein-coding genes (cytochrome b, NADH dehydrogenase subunit 4, and ATPase 6) were sequenced. Four haplotype groups with nonoverlapping geographical ranges were recovered in separate and combined analyses of the data. One clade corresponds to the distribution of P. stormi, while the remaining three comprise P. elongatus. Phylogenetic relationships among haplotype groups differ in separate analyses of the genes but converge on a well-supported topology, with P. elongatus and P. stormi as monophyletic sister taxa, in combined Maximum Parsimony and Maximum Likelihood analyses. Population genetic analyses of mismatch distributions and Tajima's D-statistic are consistent with range expansion for the largest clade within P. elongatus, covering the northern two-thirds of the species range. In contrast, the P. stormi haplotype clade and the P. elongatus clade from the southern third of the species range may have been relatively stable. Morphological boundaries between P. elongatus and P. stormi are largely congruent with mitochondrial DNA breaks and continued treatment as sister taxa is supported. Although mitochondrial DNA haplotype groups may reflect historical separation within P. elongatus, genetic barriers are incongruent with intraspecific patterns of morphological variation.     

Isolation of dinucleotide microsatellite loci from red‐backed salamander (Plethodon cinereus)

We isolated 13 variable dinucleotide microsatellites from red‐backed salamanders (Plethodon cinereus). After generating fragments using degenerate oligonucleotide primer‐polymerase chain reaction (DOP‐PCR), AC repeats were captured using biotinylated probes and streptavidin‐coated magnetic particles. Captured fragments were cloned into plasmids, screened for microsatellites with a simple PCR reaction, and select plasmids then sequenced. PCR primers were designed and optimized for robust amplification, and nine primers have been further optimized for multiplex reactions with fluorescent primers. These nine loci are variable with an average of 6.11 alleles per locus and an average heterozygosity of 0.61.     

Novel tetranucleotide microsatellite markers from the Del Norte salamander (Plethodon elongatus) with application to its sister species the Siskiyou Mountain salamander (P. stormi)

Eleven tetranucleotide microsatellite loci were developed for the Del Norte salamander (Plethodon elongatus). The loci were variably polymorphic, ranging from two to 20 alleles per locus, with expected heterozygosities ranging from 0.07 to 0.86. The loci also amplified in a congener, the Siskiyou Mountain salamander (P. stormi). The microsatellite loci will be used to assess the utility of highly polymorphic markers to assay within‐ and between‐species differentiation between these two closely related species.     

Salamander morph frequencies do not evolve as predicted in response to 40 years of climate change

The global climate is changing rapidly, yet biotic responses remain uncertain. Most studies focus on changes in species ranges or plastic responses like phenology, but adaptive evolution could be equally important. Studying evolutionary responses is challenging given limited historical data and a poor understanding of genetically variable traits under selection. We take advantage of a historical dataset to test for an adaptive response to climate change in a widespread, polymorphic amphibian, the eastern red‐backed salamander Plethodon cinereus. We resurveyed color morph frequencies across New England to test for an adaptive shift in response to climate change. We modeled historical and present‐day morph proportions as a function of climate and tested the accuracy of predictions both within and across different time periods. Our models showed moderate accuracy when predicting morph frequencies within time periods, but poor accuracy across time periods. Despite substantial changes in climate and significant relationships between morph frequency and climate variables within periods, we found no evidence for the predicted shift in morph frequencies across New England. The relationship between climate and color morph frequencies is likely more complex than originally suggested, potentially involving the interplay of additional factors such as microclimate variation, land use changes, and frequency‐dependent selection. Model extrapolation and changes in the correlation structure of climate variables also likely contributed to poor predictive ability. Evolution could provide a means to moderate the effects of climate change on many species. However, we often do not understand the direct links between climate variation, traits, and fitness. Therefore, forecasting climate‐mediated evolution remains an ongoing and important challenge for understanding climate change threats to species.     

The role of climate in the dynamics of a hybrid zone in Appalachian salamanders

I examined the potential influence of climate change on the dynamics of a previously studied hybrid zone between a pair of terrestrial salamanders at the Coweeta Hydrologic Laboratory, U.S. Forest Service, in the Nantahala Mountains of North Carolina, USA. A 16‐year study led by Nelson G. Hairston, Sr. revealed that Plethodon teyahalee and Plethodon shermani hybridized at intermediate elevations, forming a cline between ‘pure’ parental P. teyahalee at lower elevations and ‘pure’ parental P. shermani at higher elevations. From 1974 to 1990 the proportion of salamanders at the higher elevation scored as ‘pure’P. shermani declined significantly, indicating that the hybrid zone was spreading upward. To date there have been no rigorous tests of hypotheses for the movement of this hybrid zone. Using temperature and precipitation data from Coweeta, I re‐analyzed Hairston's data to examine whether the observed elevational shift was correlated with variation in either air temperature or precipitation from the same time period. For temperature, my analysis tracked the results of the original study: the proportion of ‘pure’P. shermani at the higher elevation declined significantly with increasing mean annual temperature, whereas the proportion of ‘pure’P. teyahalee at lower elevations did not. There was no discernable relationship between proportions of ‘pure’ individuals of either species with variation in precipitation. From 1974 to 1990, low‐elevation air temperatures at the Coweeta Laboratory ranged from annual means of 11.8 to 14.2 °C, compared with a 55‐year average (1936–1990) of 12.6 °C. My re‐analyses indicate that the upward spread of the hybrid zone is correlated with increasing air temperatures, but not precipitation, and provide an empirical test of a hypothesis for one factor that may have influenced this movement. My results aid in understanding the potential impact that climate change may have on the ecology and evolution of terrestrial salamanders in montane regions.     

PREDATOR PERCEPTION OF BATESIAN MIMICRY AND CONSPICUOUSNESS IN A SALAMANDER

In Batesian mimicry a palatable mimic deceives predators by resembling an unpalatable model. The evolution of Batesian mimicry relies on the visual capabilities of the potential predators, as prey detection provides the selective force driving evolutionary change. We compared the visual capabilities of several potential predators to test predictions stemming from the hypothesis of Batesian mimicry between two salamanders: the model species Notophthalmus viridescens, and polymorphic mimic, Plethodon cinereus. First, we found mimicry to be restricted to coloration, but not brightness. Second, only bird predators appeared able to discriminate between the colors of models and nonmimic P. cinereus. Third, estimates of salamander conspicuousness were background dependent, corresponding to predictions only for backgrounds against which salamanders are most active. These results support the hypothesis that birds influence the evolution of Batesian mimicry in P. cinereus, as they are the only group examined capable of differentiating N. viridescens and nonmimetic P. cinereus. Additionally, patterns of conspicuousness suggest that selection from predators may drive the evolution of conspicuousness in this system. This study confirms the expectation that the visual abilities of predators may influence the evolution of Batesian mimicry, but the role of conspicuousness may be more complex than previously thought.     

A PARAMETRIC METHOD FOR ASSESSING DIVERSIFICATION‐RATE VARIATION IN PHYLOGENETIC TREES

Phylogenetic hypotheses are frequently used to examine variation in rates of diversification across the history of a group. Patterns of diversification‐rate variation can be used to infer underlying ecological and evolutionary processes responsible for patterns of cladogenesis. Most existing methods examine rate variation through time. Methods for examining differences in diversification among groups are more limited. Here, we present a new method, parametric rate comparison (PRC), that explicitly compares diversification rates among lineages in a tree using a variety of standard statistical distributions. PRC can identify subclades of the tree where diversification rates are at variance with the remainder of the tree. A randomization test can be used to evaluate how often such variance would appear by chance alone. The method also allows for comparison of diversification rate among a priori defined groups. Further, the application of the PRC method is not restricted to monophyletic groups. We examined the performance of PRC using simulated data, which showed that PRC has acceptable false‐positive rates and statistical power to detect rate variation. We apply the PRC method to the well‐studied radiation of North American Plethodon salamanders, and support the inference that the large‐bodied Plethodon glutinosus clade has a higher historical rate of diversification compared to other Plethodon salamanders.     

Evaluating within‐population variability in behavior and demography for the adaptive potential of a dispersal‐limited species to climate change

Multiple pathways exist for species to respond to changing climates. However, responses of dispersal‐limited species will be more strongly tied to ability to adapt within existing populations as rates of environmental change will likely exceed movement rates. Here, we assess adaptive capacity in Plethodon cinereus, a dispersal‐limited woodland salamander. We quantify plasticity in behavior and variation in demography to observed variation in environmental variables over a 5‐year period. We found strong evidence that temperature and rainfall influence P. cinereus surface presence, indicating changes in climate are likely to affect seasonal activity patterns. We also found that warmer summer temperatures reduced individual growth rates into the autumn, which is likely to have negative demographic consequences. Reduced growth rates may delay reproductive maturity and lead to reductions in size‐specific fecundity, potentially reducing population‐level persistence. To better understand within‐population variability in responses, we examined differences between two common color morphs. Previous evidence suggests that the color polymorphism may be linked to physiological differences in heat and moisture tolerance. We found only moderate support for morph‐specific differences for the relationship between individual growth and temperature. Measuring environmental sensitivity to climatic variability is the first step in predicting species' responses to climate change. Our results suggest phenological shifts and changes in growth rates are likely responses under scenarios where further warming occurs, and we discuss possible adaptive strategies for resulting selective pressures.     

Effect of search method and age class on mark-recapture parameter estimation in a population of red-backed salamanders

Mark-recapture methods are generally considered to more accurately reflect population trends than count data, which is especially important for indicator species. Terrestrial salamanders are often used as indicators of forest ecosystem health and may be monitored through diurnal cover object searches or nocturnal activity searches. Our goal was to determine whether search method affected encounter probabilities, whether these probabilities differed between age classes, and whether the inclusion of search method in mark-recapture models affected abundance estimates. We used program MARK to analyze 3 years of red-backed salamander (Plethodon cinereus) mark-recapture data using Pollock’s robust design taken from a 144 m² plot. Initial encounter probabilities during night searches were consistently greater than during diurnal cover object searches. As a result, inclusion of search method in models sometimes affected abundance estimates. There was no difference between adult and juvenile encounter probabilities nor were abundance estimates affected by inclusion of age class, but there was yearly variation in the juvenile abundance estimates. For these reasons, we recommend that sampling of terrestrial salamanders include nocturnal activity searches and be conducted over multiple years. Monitoring programs of other species should take into account whether the selected search method(s) may be more likely to sample different subsets of the population of interest and how this may restrict their inferences.     

Female philopatry and male-biased dispersal in a direct-developing salamander, Plethodon cinereus

The local resource competition hypothesis and the local mate competition hypothesis were developed based on avian and mammalian systems to explain sex-biased dispersal. Most avian species show a female bias in dispersal, ostensibly due to resource defence, and most mammals show a male bias, ostensibly due to male-male competition. These findings confound phylogeny with mating strategy; little is known about sex-biased dispersal in other taxa. Resource defence and male-male competition are both intense in Plethodon cinereus, a direct-developing salamander, so we tested whether sex-biased dispersal in this amphibian is consistent with the local resource competition hypothesis (female-biased) or the local mate competition hypothesis (male-biased). Using fine-scale genetic spatial autocorrelation analyses, we found that females were philopatric, showing significant positive genetic structure in the shortest distance classes, with stronger patterns apparent when only territorial females were tested. Males showed no spatial genetic structure over the shortest distances. Mark-recapture observations of P. cinereus over 5 years were consistent with the genetic data: males dispersed farther than females during natal dispersal and 44% of females were recaptured within 1 m of their juvenile locations. We conclude that, in this population of a direct-developing amphibian, females are philopatric and dispersal is male-biased, consistent with the local mate competition hypothesis.     

Primary and secondary somatosensory projections in direct-developing Plethodontid Salamanders

In three species of plethodontid salamanders (Plethodon jordani, Hydromantes italicus, and Bolitoglossa subpalmata), primary and secondary somatosensory pathways were investigated by means of tract-tracing in vivo and in vitro using biocytin, horseradish peroxidase, and neurobiotin. Afferent sensory fibers of cranial nerves V, VII, and X and the brachial nerve run in the dorsal funiculus of the medulla oblongata and spinal cord. Fibers ascend to the level of, but do not enter, the cerebellum. In the caudal medulla oblongata, sensory tracts of the cranial nerves descend in a dorsal and a dorsolateral bundle and reach the level of the fourth spinal nerve. Two bundles are likewise formed by spinal afferent fibers, which descend to the level of the seventh spinal nerve. Secondary somatosensory projections ascend in contralateral ventral, contralateral lateral, and ipsilateral lateral tracts, the latter two corresponding to the spinal lemniscal tracts of Herrick. These tracts reach the cerebellum, mesencephalic, and diencephalic targets (tegmentum, torus, tectum, tuberculum posterius, pretectum, and ventral thalamus) ipsi- and contra-laterally. The projection to the tectum is confined to fiber layer 4. Fibers of the ascending tracts cross in the cerebellar and tectal commissure. Our study demonstrates that the ascending secondary somatosensory pathways of plethodontid salamanders differ remarkably from those of other amphibians. J. Morphol. 238:307–326, 1998. © 1998 Wiley-Liss, Inc.     

Genetic variation in Plethodon cinereus and Plethodon hubrichti from in and around a contact zone

Climate change poses several challenges to biological communities including changes in the frequency of encounters between closely related congeners as a result of range shifts. When climate change leads to increased hybridization, hybrid dysfunction or genetic swamping may increase extinction risk—particularly in range‐restricted species with low vagility. The Peaks of Otter Salamander, Plethodon hubrichti, is a fully terrestrial woodland salamander that is restricted to ~18 km of ridgeline in the mountains of southwestern Virginia, and its range is surrounded by the abundant and widespread Eastern Red‐backed Salamander, Plethodon cinereus. In order to determine whether these two species are hybridizing and how their range limits may be shifting, we assessed variation at eight microsatellite loci and a 1,008 bp region of Cytochrome B in both species at allopatric reference sites and within a contact zone. Our results show that hybridization between P. hubrichti and P. cinereus either does not occur or is very rare. However, we find that diversity and differentiation are substantially higher in the mountaintop endemic P. hubrichti than in the widespread P. cinereus, despite similar movement ability for the two species as assessed by a homing experiment. Furthermore, estimation of divergence times between reference and contact zone populations via approximate Bayesian computation is consistent with the idea that P. cinereus has expanded into the range of P. hubrichti. Given the apparent recent colonization of the contact zone by P. cinereus, future monitoring of P. cinereus range limits should be a priority for the management of P. hubrichti populations.     

Sperm aggregations in the spermatheca of the red back salamander (Plethodon cinereus)

The spermatheca of Plethodon cinereus is a compound tubular gland that stores sperm from mating in early spring (March–April) to oviposition in summer (June–July). The seasonal variation of sperm storage in this species has previously been studied by light and transmission electron microscopy. In this paper, sperm aggregations, interaction of sperm with the spermathecal epithelium, and spermathecal secretions are studied using scanning electron microscopy. Within spermathecal tubules, relatively small groups of sperm are aligned along their entire lengths in parallel arrays. This pattern is similar to other plethdontids with complex spermathecae. Lumina of spermathecal tubules are filled with secretory material in April prior to the arrival of sperm, and after sperm appear, a coating of secretory material persists on the apices of the spermathecal epithelium. Sperm peripheral to the central luminal mass can become embedded in the secretory matrix or pushed deeper into the spermathecal epithelium. The spermathecal secretions may serve to attract and prolong the viability of sperm, but sperm that become enmeshed in the secretions or epithelium are phagocytized. Sperm and spermathecal secretions are largely absent after ovulation and in summer months, and new secretory vacuoles are formed in fall, although mating does not occur until spring.     

Lead-phase and red-stripe color morphs of red-backed salamanders Plethodon cinereus differ in hematological stress indices: A consequence of differential predation pressure?

Throughout the animal kingdom there are species that have two or more phenotypic forms or 'morphs', and many of these are amphibians. In North America, the red-backed salamander Plethodon cinereus can have either a red dorsal stripe or no dorsal stripe (lead-phase form), and evidence to date indicates the lead-phase form incurs a greater number of attacks from predators. In a recent collection of 51 P. Cinereus, blood smears of both color morphs (35 red-stripe, 16 lead-phase) were examined to obtain numbers of circulating leukocytes (via light microscopy), which can be used to indirectly estimate levels of stress hormones in vertebrates via a 'hematological stress index', which is the ratio between the number of two leukocyte types (neutrophils and lymphoeytes). Our results showed that lead-phase salamanders tended to have greater numbers of circulating neutrophils and lower numbers of circulating iymphocytes than red-stripe morphs, leading to higher average neutrophil-lymphocyte ratios in lead-phase individuals. Since the salamanders were held (refrigerated) for 7 days before sampling, we cannot be certain if this effect is a stress reaction to the captivity or the normal level for this morph. However comparison with two sets of related salamanders that were captured and sampled immediately indicates the red-stripe salamanders were either not stressed from the captivity at all, or their white blood cell distributions had returned to normal after 7 days of captivity. Taken together, our results indicate that lead-phase forms of P. Cinereus have higher stress levels than the red-stripe forms, which may be a consequence of their higher exposure to, and/or attacks from, predators. They may also indicate that the lead-phase form is less-suited to captivity than the red-stripe form of this species.     

Behavioral interactions between salamanders and centipedes: competition in divergent taxa

Red-backed salamanders, Plethodon cinereus, use territorialadvertisement in the form of agonistic displays and pheromonalscent marking as a mechanism for intraspecific interferencecompetition. Although ecological and behavioral interactionsamong species of salamanders have been well studied, littleis known about the interactions between territorial P. cinereusand other ecologically similar species, such as large predatoryinvertebrates. Our field data indicate that P. cinereus anda large syntopic centipede, Scolopocryptops sexspinosus, exhibitnegative spatial associations in natural habitats, possiblyindicating interspecific territoriality. Only seven instancesof salamander/centipede co-occurrence were recorded from a fieldsample of 247 occupied cover objects. Cover object size waspositively correlated with salamander SVL (tip of the snoutto the anterior end of the cloaca), but there was no correlationof cover object size to centipede length. Data on the abilityof P. cinereus to differentiate among chemicals on the substratesuggest that visual cues are not necessary to elicit a territorialresponse from intruding salamanders. Although in laboratorytrials salamanders behaved similarly toward intruders of bothspecies, biting was directed only toward centipedes. Salamandersspent significantly more time approaching centipedes than theydid approaching other salamanders. Approach behavior was oftenassociated with nose tapping and may be an investigative, ratherthan aggressive, behavior. We suggest that territorial P. cinereusrespond similarly to intruding salamanders and centipedes, butthat they escalate more readily to biting centipedes becauseS. sexspinosus is sightless and thus unable to respond to visualsignals.     

Heterozygosity–behavior and heterozygosity–fitness correlations in a salamander with limited dispersal

Inbreeding and genetic drift can decrease genetic heterozygosity, and this low heterozygosity can depress fitness, resulting in heterozygosity–fitness correlations (HFCs). HFCs are typically small in magnitude, a result often attributed to power of the analyses. Animal behaviors often affect fitness and are often heritable to some degree. We hypothesized that heterozygosity influences behavior, which, in turn, potentially influences fitness. Specifically, in red-backed salamanders, Plethodon cinereus, which have limited dispersal and the potential for inbreeding, we tested whether heterozygosity, as estimated from six microsatellite loci, affected home range size, juvenile growth, and survival. We found that salamanders with higher heterozygosity had larger home ranges and grew faster, which is indicative of reproductive success. However, we found no effects of heterozygosity on survival. We conclude that, because activity in P. cinereus is tightly linked to food uptake and mass gain, heterozygosity influences growth via effects on foraging behavior. Future research should investigate how the relationship between heterozygosity, behavior, and fitness may be affected or mediated by endocrine or immune systems.     

Salamander diversity of reforested abandoned surface coal mines in the Appalachian Region,U.S.A.

Surface coal mining is a pervasive form of land use and leading cause of deforestation in the Appalachian Region of the United States. Contemporary reclamation practices tend to produce strongly altered landscapes that inhibit forest succession, so there is a need for the development of alternative methods. Surface coal mines that preceded modern reclamation laws, known as abandoned mine lands, have become reforested in many locations and may serve as a model for improving the restoration of forests in postmining landscapes. I compared salamander diversity between 5 forests growing atop abandoned surface coal mines to 5 similar aged reference forests in southeastern Ohio using time constrained transect searches in autumn 2012 and spring 2013. The survey revealed 1,480 salamander individuals, representing 11 different species from 3 different families. The vast majority of salamanders (93%) were from the family Plethodontidae. Salamander species richness was significantly higher in reference forests than in mined forests, but this difference was contingent on the presence of stream‐dependent species in the reference forests. Salamander relative abundance was not significantly different between reference and mined forests. This survey demonstrates that reforested abandoned surface coal mines can provide suitable habitat for many forest‐dependent, terrestrial salamander species, but impaired water quality of headwater streams may limit the colonization of these forests by stream‐dependent species.     

Multiple paternity in a salamander with socially monogamous behaviour

In the majority of birds and mammals, social monogamy is not congruent with genetic monogamy. No research to date has compared social and genetic monogamy in amphibians. We analysed paternity in clutches of red-backed salamanders (Plethodon cinereus), a species in which social monogamy has been demonstrated in the laboratory, and 28% of individuals in the forest are found in male-female pairs in the noncourtship season. We collected 16 clutches of eggs of P. cinereus in the southern Appalachian Mountains of Virginia and collected tail clippings from attending mothers. We genotyped embryos and adults at five microsatellite loci in order to analyse paternity of clutches. Most clutches (84.6%) had multiple sires, with two to three sires per clutch. In this study, 25% of clutches had males in addition to females attending eggs. None of the mothers of these clutches were genetically monogamous. All attending males sired some of the offspring in the clutch that they attended (between 9% and 50%) but never sired a majority in that clutch. We conclude that, at least in this population, social monogamy in P. cinereus is not concomitant with genetic monogamy.     

Demographic consequences of climate variation along an elevational gradient for a montane terrestrial salamander

Climate change represents a significant threat to amphibians. However, for many species, the relationship between demography and climate is unknown, which limits predictive models. Here, we describe the life history variation of Plethodon montanus using capture–recapture data over a period of 4 years, along an elevational gradient to determine how survival and growth vary with climate, and how these relationships vary with elevation. We used a hierarchical model to estimate asymptotic size and growth rate and a spatial Cormack-Jolly-Seber model to estimate the probability of capture and survival and dispersal variance. We found that during the active season, growth and survival rates are both positively affected by precipitation; however, while survival was positively affected by temperature at all elevations, higher temperatures led to a decrease in growth at lower elevations, while at higher elevations the opposite was true. During the inactive season, we found reduced growth rates, whereas survival was lower compared with the active at lower elevations but was higher at higher elevations. Increased inactive season temperatures resulted in decreased survival while we found that temperature, amount of snow, and elevation interacted to influence survival. At low elevations, which were warmer, survival decreased with increasing snowfall but at higher elevations, survival generally increased with increasing snowfall. Our results demonstrate that understanding how the environment can affect salamander demography to develop mechanistic models will require knowledge of the actual environmental conditions experienced by a given population as well as an understanding of the overall differences in climate at a given site.     

Pheromonal markers as indicators of parasite load: parasite-mediated behavior in salamanders (Plethodon angusticlavius)

Assessment of parasite load of conspecifics may be important during social interactions such as courtship and aggressive encounters. We used a correlational study to test whether pheromonal markers can be used to assess parasite load of conspecifics, and whether the parasite load of the pheromone receivers affected their responses. We tested the responses of parasitized and nonparasitized Ozark zigzag salamanders, Plethodon angusticlavius, to territorial markers (fecal pellets) from conspecific males. Males and females were simultaneously exposed to fecal pellets placed in front of two artificial burrows located at the opposite ends of their chambers. The treatments were (1) fecal pellet of male with low parasite load versus fecal pellet of male with high parasite load, (2) fecal pellet of male with low parasite load versus control pellet (chemical blank), and (3) fecal pellet of male with high parasite load versus control pellet. Nonparasitized females spent significantly more time near fecal pellets of males with low parasite load in test condition 1, whereas the behavior of parasitized females was not significantly different from random in any test condition. Males responded differently to treatments only in condition 3; males with low parasite loads spent significantly more time near control pellets, whereas males with high parasite loads spent significantly more time near fecal pellets of males with high parasite loads. This study demonstrates that pheromonal markers may be used for assessment of parasite load of conspecifics and that responses by both males and females may be influenced by their own level of infection. Received: 21 January 1999 / Received in revised form: 17 March 2000 / Accepted: 13 November 2000