Insights into the mating habits of the tiger salamander (Ambystoma tigrinum tigrinum) as revealed by genetic parentage analyses

Does the potential strength of sexual selection differ between mating systems with and without defensive behaviours? A meta-analysis

The Darwin-Bateman paradigm predicts that females enhance their fitness by being choosy and mating with high-quality males, while males should compete to mate with as many females as possible. In many species, males enhance their fitness by defending females and/or resources used by females. That is, males directly defend access to mating opportunities. However, paternity analyses have repeatedly shown that females in most species mate polyandrously, which contradicts traditional expectations that male defensive behaviours lead to monandry. Here, in an extensive meta-analysis, encompassing 109 species and 1026 effect sizes from across the animal kingdom, we tested if the occurrence of defensive behaviours modulates sexual selection on females and males. If so, we can illuminate the extent to which males really succeed in defending access to mating and fertilisation opportunities. We used four different indices of the opportunity for sexual selection that comprise pre-mating and/or post-mating episodes of selection. We found, for both sexes, that the occurrence of defensive behaviours does not modulate the potential strength of sexual selection. This implies that male defensive behaviours do not predict the true intensity of sexual selection. While the most extreme levels of sexual selection on males are in species with male defensive behaviours, which indicates that males do sometimes succeed in restricting females' re-mating ability (e.g. elephant seals, Mirounga leonina), estimates of the opportunity for sexual selection vary greatly across species, regardless of whether or not defensive behaviours occur. Indeed, widespread polyandry shows that females are usually not restricted by male defensive behaviours. In addition, our results indicate that post-mating episodes of selection, such as cryptic female choice and sperm competition, might be important factors modulating the opportunity for sexual selection. We discuss: (i) why male defensive behaviours fail to lower the opportunity for sexual selection among females or fail to elevate it for males; (ii) how post-mating events might influence sexual selection; and (iii) the role of females as active participants in sexual selection. We also highlight that inadequate data reporting in the literature prevented us from extracting effect sizes from many studies that had presumably collected the relevant data.

Do female amphibians and reptiles have greater reproductive output if they have more mates?

Abstract

In general, males mate with multiple females to increase individual reproductive success. Whether or not, and under what circumstances, females benefit from multiple mating has been less clear. Our review of 154 studies covering 184 populations of amphibians and reptiles showed that polyandry was widespread and variable among and within taxonomic groups. We investigated whether amphibian and reptile females had greater reproductive output as the number of sires for offspring increased. Meta-analysis revealed significant heterogeneity in the dataset of all taxa. Expected heterozygosity was a significant moderator (covariate) of positive relationships between female reproductive output and the number of sires, but a sensitivity test showed the result was tenuous. Significant heterogeneity remained despite controlling for expected heterozygosity and other variables but was resolved for most taxonomic groups with subgroup meta-analyses. Subgroup meta-analyses showed that only female salamanders (Caudata) had significantly greater reproductive output with an increased number of sires. For many species of Caudata, males cannot coerce females into accepting spermatophores. We therefore suggest that if females control the number of matings, they can use polyandry to increase their fitness. Caudata offers ideal models with which to test this hypothesis and to explore factors enabling and maintaining the evolution of female choice. Outstanding problems may be addressed by expanding taxonomic coverage and data collection and improving data reporting.

Significance Statement

Many factors and combinations of factors drive polyandry. Whether or not females benefit from mating with more than one male remains equivocal. Focusing on amphibians and reptiles, our analyses demonstrate that female salamanders produced more offspring when mated with multiple males, whereas this was not the case for reptiles. Unlike many other species in our dataset, the polyandrous female salamanders fully control sperm intake and have chosen to mate multiple times. We further highlight problems and key directions for future research in the field.

The evolution of pueriparity maintains multiple paternity in a polymorphic viviparous salamander

The reduction in fecundity associated with the evolution of viviparity may have far-reaching implications for the ecology, demography, and evolution of populations. The evolution of a polygamous behaviour (e.g. polyandry) may counteract some of the effects underlying a lower fecundity, such as the reduction in genetic diversity. Comparing patterns of multiple paternity between reproductive modes allows us to understand how viviparity accounts for the trade-off between offspring quality and quantity. We analysed genetic patterns of paternity and offspring genetic diversity across 42 families from two modes of viviparity in a reproductive polymorphic species, Salamandra salamandra. This species shows an ancestral (larviparity: large clutches of free aquatic larvae), and a derived reproductive mode (pueriparity: smaller clutches of larger terrestrial juveniles). Our results confirm the existence of multiple paternity in pueriparous salamanders. Furthermore, we show the evolution of pueriparity maintains, and even increases, the occurrence of multiple paternity and the number of sires compared to larviparity, though we did not find a clear effect on genetic diversity. High incidence of multiple paternity in pueriparous populations might arise as a mechanism to avoid fertilization failures and to ensure reproductive success, and thus has important implications in highly isolated populations with small broods.

Darwinian sex roles confirmed across the animal kingdom

Since Darwin's conception of sexual selection theory, scientists have struggled to identify the evolutionary forces underlying the pervasive differences between male and female behavior, morphology, and physiology. The Darwin-Bateman paradigm predicts that anisogamy imposes stronger sexual selection on males, which, in turn, drives the evolution of conventional sex roles in terms of female-biased parental care and male-biased sexual dimorphism. Although this paradigm forms the cornerstone of modern sexual selection theory, it still remains untested across the animal tree of life. This lack of evidence has promoted the rise of alternative hypotheses arguing that sex differences are entirely driven by environmental factors or chance. We demonstrate that, across the animal kingdom, sexual selection, as captured by standard Bateman metrics, is indeed stronger in males than in females and that it is evolutionarily tied to sex biases in parental care and sexual dimorphism. Our findings provide the first comprehensive evidence that Darwin's concept of conventional sex roles is accurate and refute recent criticism of sexual selection theory.

The effects of contig length and depth on the estimation of SNP frequencies, and the relative abundance of SNPs in protein-coding and non-coding transcripts of tiger salamanders (Ambystoma tigrinum)

BACKGROUND

Next-generation sequencing methods have contributed to rapid progress in the fields of genomics and population genetics. Using this high-throughput and cost-effective technology, a number of studies have estimated single nucleotide polymorphism (SNP) frequency by calculating the mean number of SNPs per unit sequence length (e.g., mean SNPs/kb). However, both read length and contig depth are highly variable and thus raise doubt about simple methods of SNP frequency estimation.

RESULTS

We used 454 pyrosequencing to identify 2,980 putative SNPs in the eastern tiger salamander (Ambystoma tigrinum tigrinum) transcriptome, then constructed analytical models to estimate SNP frequency. The model which considered only contig length (i.e., the method employed in most published papers) was evaluated with very poor likelihood. Our most robust model considered read depth as well as contig length, and was 7.5 × 10⁵⁵ times more likely than the length-only model. Using this novel modeling approach, we estimated SNP frequency in protein-coding (mRNA) and non-coding transcripts (e.g., small RNAs). We found little difference in SNP frequency in the contigs, but we found a trend of a higher frequency of SNPs in long contigs representing non-coding transcripts relative to protein-coding transcripts. These results support the hypothesis that long non-coding transcripts are less conserved than long protein-coding transcripts.

CONCLUSIONS

A modeling approach (i.e., using multiple model construction and model selection approaches) can be a powerful tool for identifying selection on specific functional sequence groups by comparing the frequency and distribution of polymorphisms.

Comparative transcriptomics and gene expression in larval tiger salamander (Ambystoma tigrinum) gill and lung tissues as revealed by pyrosequencing

Biologists are beginning to unravel the complexities of gene expression in model organisms by studying the transcriptome, the complement of genes that are transcribed in a given tissue. It is unclear, however, if findings from model systems apply to non-model organisms because of environmental effects on gene expression. Furthermore, there have been few efforts to quantify how transcriptome or gene expression varies across individuals and across tissues in natural environments. Herein, we describe transcriptomic profiling of gene expression in lung and gill tissue of three larval tiger salamanders. We do so with a hierarchical experimental design that captures variation in expression among genes, among tissues, and among individuals. Using 454 pyrosequencing, we produced high-quality sequence data of 59 megabases and assembled ~200,000 reads into 19,501 contigs. These contigs BLASTed to 3,599 transcripts, of which 721 were expressed in both tissues, 1,668 were unique to gill, and 1,210 unique to lung. Our data showed tissue-specific patterns in gene expression level with variation among transcripts and individuals. We identified genes and gene ontology terms related to respiration and compared their relative expression levels between gill and lung tissues. We also found evidence of exogenous genes associated with larval salamanders, and we identified ~1400 potential molecular markers (microsatellites and single nucleotide polymorphisms) that are associated with expressed genes. Given the tissue-specific differences we observed in transcriptomes, these data reinforce the idea that changes in gene expression serve as a primary mechanism underlying phenotypic plasticity.

A model for estimating the minimum number of offspring to sample in studies of reproductive success

Molecular parentage permits studies of selection and evolution in fecund species with cryptic mating systems, such as fish, amphibians, and insects. However, there exists no method for estimating the number of offspring that must be assigned parentage to achieve robust estimates of reproductive success when only a fraction of offspring can be sampled. We constructed a 2-stage model that first estimated the mean (μ) and variance (v) in reproductive success from published studies on salmonid fishes and then sampled offspring from reproductive success distributions simulated from the μ and v estimates. Results provided strong support for modeling salmonid reproductive success via the negative binomial distribution and suggested that few offspring samples are needed to reject the null hypothesis of uniform offspring production. However, the sampled reproductive success distributions deviated significantly (χ(2) goodness-of-fit test p value < 0.05) from the known simulated reproductive success distribution at rates often >0.05 and as high as 0.24, even when hundreds of offspring were assigned parentage. In general, reproductive success patterns were less accurate when offspring were sampled from cohorts with larger numbers of parents and greater variance in reproductive success. Our model can be reparameterized with data from other species and will aid researchers in planning reproductive success studies by providing explicit sampling targets required to accurately assess reproductive success.

Comparing the performance of analytical techniques for genetic PARENTAGE of half-sib progeny arrays

The prevalence of female multiple mating in natural populations is important for many questions in mating system evolution. Several statistical techniques use genetic data to estimate the number of fathers that contribute gametes to broods, but they have not been widely compared to assess the magnitude of differences in their performance. With a combination of new data and reanalysis of previously published data, we compared five analytical approaches: (1) allele-counting, (2) parental reconstruction in GERUD, (3) a Bayesian probability model to estimate the frequency of multiple mating (FMM), (4) computer simulations based on population allele frequencies in HAPLOTYPES and (5) Bayesian parental reconstruction in PARENTAGE. The results show that choice of analysis technique can significantly affect estimates of sire number. Estimates from GERUD conformed exactly to results obtained from strict exclusion of potential sires in an experimental context. However, estimates yielded by HAPLOTYPES and PARENTAGE sometimes exceeded the numbers from GERUD by as much as 120 and 55%, respectively. We recommend GERUD over these other approaches for most purposes because of its accuracy and consistency in this analysis. Our novel genetic data set allowed us to investigate the extent and frequency of multiple paternity in a marbled salamander (Ambystoma opacum) population in South Carolina, USA. A. opacum contrasted with other salamander species by having relatively low levels of multiple paternity (only 31-54% compared with 71-96%). Although A. opacum had the lowest level of multiple paternity under all analytical approaches used here, the magnitude of differences among species varied.

Condition-dependent mate choice and a reproductive disadvantage for MHC-divergent male tiger salamanders

Major histocompatibility complex (MHC) alleles likely have adaptive value because of overdominance, in which case MHC heterozygous individuals have increased fitness relative to homozygotes. Because of this potential benefit, the evolution of sexual reproduction between MHC-divergent individuals (i.e. negative assortative mating, NAM) may be favoured. However, the strongest evidence for MHC-based NAM comes from inbred animals, and context-dependent mating preferences have rarely been evaluated although they often occur in nature. We assessed the extent MHC-based mating preferences among wild tiger salamanders (Ambystoma tigrinum) using multiple molecular approaches. We genotyped 102 adults and 864 larvae from 36 breeding trials at both microsatellite and MHC loci. Parentage analysis revealed that reproductive success among males was positively associated with increased tail length and that with respect to the focal female, MHC-similar males sired a significantly higher number of offspring than more dissimilar males. This trend was consistent, even under context-dependent scenarios that favour traditional MHC-based NAM. These results suggest that the most MHC-divergent males may be at a reproductive disadvantage in pairwise breeding trials. Our data add to a growing body of evidence that suggests where it exists, MHC-based choice is probably dynamic and mediated by many factors that vary in the wild, notably signals from other indicator traits and by the quality and quantity of potential mates.

Inferring population history and demography using microsatellites, mitochondrial DNA, and major histocompatibility complex (MHC) genes

Microsatellites and mitochondrial DNA (mtDNA) have traditionally been used in population genetics because of their variability and presumed neutrality, whereas genes of the major histocompatibility complex (MHC) are increasingly of interest because strong selective pressures shape their standing variation. Despite the potential for MHC genes, microsatellites, and mtDNA sequences to complement one another in deciphering population history and demography, the three are rarely used in tandem. Here we report on MHC, microsatellite, and mtDNA variability in a single large population of the eastern tiger salamander (Ambystoma tigrinum tigrinum). We use the mtDNA mismatch distribution and, on microsatellite data, the imbalance index and bottleneck tests to infer aspects of population history and demography. Haplotype and allelic variation was high at all loci surveyed, and heterozygosity was high at the nuclear loci. We find concordance among neutral molecular markers that suggests our study population originated from post-Pleistocene expansions of multiple, fragmented sources that shared few migrants. Differences in N(e) estimates derived from haploid and diploid genetic markers are potentially attributable to secondary contact among source populations that experienced rapid mtDNA divergence and comparatively low levels of nuclear DNA divergence. We find strong evidence of natural selection acting on MHC genes and estimate long-term effective population sizes (N(e)) that are very large, making small selection intensities significant evolutionary forces in this population.

Variation in multiple paternity in natural populations of a free-spawning marine invertebrate

For free-spawning marine invertebrates, fertilization processes control the genetic diversity of offspring. Each egg can potentially be fertilized by a sperm from a different male, and hence genetic diversity within a brood varies with levels of multiple paternity. Yet, few studies have characterized the frequency of multiple paternity in natural spawns. We analysed patterns of multiple paternity in two populations of the colonial ascidian Botryllus schlosseri using microsatellites. Because previous studies have shown that at moderate to high population densities, competition among male-phase B. schlosseri colonies results in the nearest male dominating the paternity of a brood, we specifically tested the effect of population density on patterns of paternity. Paternity was estimated using three multilocus indices: minimum number of fathers, counts of sperm haplotypes, and effective paternity (K(E)). Multiple paternity was evident in more than 92% of the broods analysed, but highly variable, with a few broods displaying unequal contributions of different males. We found no effect of population density on multiple paternity, suggesting that other factors may control paternity levels. Indirect benefits from increasing the genetic diversity of broods are a possible explanation for the high level of multiple paternity in this 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.