A genetic analogue of 'mark-recapture' methods for estimating population size: an approach based on molecular parentage assessments

Estimation of breeding population size using DNA-based pedigree reconstruction in brown bears

Robust estimates of demographic parameters are critical for effective wildlife conservation and management but are difficult to obtain for elusive species. We estimated the breeding and adult population sizes, as well as the minimum population size, in a high-density brown bear population on the Shiretoko Peninsula, in Hokkaido, Japan, using DNA-based pedigree reconstruction. A total of 1288 individuals, collected in and around the Shiretoko Peninsula between 1998 and 2020, were genotyped at 21 microsatellite loci. Among them, 499 individuals were identified by intensive genetic sampling conducted in two consecutive years (2019 and 2020) mainly by noninvasive methods (e.g., hair and fecal DNA). Among them, both parents were assigned for 330 bears, and either maternity or paternity was assigned to 47 and 76 individuals, respectively. The subsequent pedigree reconstruction indicated a range of breeding and adult (≥4 years old) population sizes: 128-173 for female breeders and 66-91 male breeders, and 155-200 for female adults and 84-109 male adults. The minimum population size was estimated to be 449 (252 females and 197 males) in 2019. Long-term continuous genetic sampling prior to a short-term intensive survey would enable parentage to be identified in a population with a high probability, thus enabling reliable estimates of breeding population size for elusive species.

Plastron color patterns allows for individual photo-identification in two different chelonian species

Natural marks have increasingly been used as a tool for individual identification in capture-mark-recapture techniques. Photo-identification is a non-invasive alternative to traditional marking techniques, allowing individual recognition of species through time and space. We tested the APHIS (Automatic Photo Identification) software as a software capable of identifying individuals of Testudo hermanni (Gmelin, 1789) and Emys orbicularis (Linneaus, 1758) in different populations during capture-release sessions in the field based on plastron color patterns, since they can be used as natural marks for identification. For this individual identification, SPM (Spot Pattern Matching) and ITM (Image Template Matching) procedures were tested; achieving 100 % success of individuals recognized in both procedures and visually verified comparing the images. However, ITM procedure was more efficient at recognizing recaptures than SPM because ITM allowed faster recaptures verification, since most of the matches were directly placed on the first position on the candidates list. Previous studies have used photo-identification on freshwater or sea turtles, but never with terrestrial tortoise species. Consequently, it can be corroborated that APHIS is a competent and efficient software considering photo-identification of T. hermanni and E. orbicularis and can be applied with close species with similar and unique individual color patterns in their plastron.

Estimating contemporary migration numbers of adults based on kinship relationships in iteroparous species

This paper describes the development of estimators for the contemporary migration number and rate of adults between two populations in iteroparous species. The proposed estimators are based on known half-sibling (HS) and/or parent-offspring (PO) relationships observed between populations across breeding seasons. The rationale is that HS and PO pairs exhibit information about the occurrence frequency of parental movements during the breeding interval. The proposed method allows for variance in the average number of offspring per parent within and between populations. In addition, coupled with the PO pairs found within the population, the estimators can be obtained using only genetic data. Generally, a sample size representing the square root of the population size is required to obtain meaningful migration information. We describe a detailed evaluation of the performance of the proposed estimators by running an individual-based model, and the results provide guidance regarding sample sizes to ensure the required accuracy and precision. In addition, given that there are few effective methods to estimate adult movement (especially when populations cannot be genetically distinct), we discuss the usefulness of the proposed kinship assignment method in terms of conservation biology and wildlife management.

Nearly unbiased estimator of contemporary Ne /N based on kinship relationships

This study develops a nearly unbiased estimator of the ratio of the contemporary effective mother size to the census size in a population, as a proxy of the ratio of contemporary effective size (or effective breeding size) to census size (N_e /N or N_b /N). The proposed estimator is based on both known mother-offspring (MO) and maternal-sibling (MS) relationships observed within the same cohort, in which sampled individuals in the cohort probably share MO relationships with sampled mothers. The rationale is that the frequency of MO and MS pairs contains information regarding the contemporary effective mother size and the (mature) census size, respectively. Therefore, the estimator can be obtained only from genetic data. We also evaluate the performance of the estimator by running an individual-based model. The results of this study provide the following: (a) parameter range for satisfying the unbiasedness, and (b) guidance for sample sizes to ensure the required accuracy and precision, especially when the order of the ratio is available. Furthermore, the results demonstrate the usefulness of a sibship assignment method for genetic monitoring, providing insights for interpreting environmental and/or anthropological factors fluctuating N_e /N (or N_b /N), especially in the context of conservation biology and wildlife management.

Joint estimation of growth and survival from mark-recapture data to improve estimates of senescence in wild populations

Understanding age-dependent patterns of survival is fundamental to predicting population dynamics, understanding selective pressures, and estimating rates of senescence. However, quantifying age-specific survival in wild populations poses significant logistical and statistical challenges. Recent work has helped to alleviate these constraints by demonstrating that age-specific survival can be estimated using mark-recapture data even when age is unknown for all or some individuals. However, previous approaches do not incorporate auxiliary information that can improve age estimates of individuals. We introduce a survival estimator that combines a von Bertalanffy growth model, age-specific hazard functions, and a Cormack-Jolly-Seber mark-recapture model into a single hierarchical framework. This approach allows us to obtain information about age and its uncertainty based on size and growth for individuals of unknown age when estimating age-specific survival. Using both simulated and real-world data for two painted turtle (Chrysemys picta) populations, we demonstrate that this additional information substantially reduces the bias of age-specific hazard rates, which allows for the testing of hypotheses related to aging. Estimating patterns of senescence is just one practical application of jointly estimating survival and growth; other applications include obtaining better estimates of the timing of recruitment and improved understanding of life-history trade-offs between growth and survival.

Mating System, Reproductive Success, and Sexual Selection in Bluntnose Klipfishes (Clinus cottoides)

A critical part of the sexual selection process in animals is the genetic mating system. Quantifying mating systems, especially in species with cryptic life histories can be challenging. One approach is to use genotypic markers and accurate parentage analysis, along with methods to account for bias when sampling natural populations, to calculate sexual selection metrics derived from Bateman's principles. In this study, 3 microsatellites were used to genotype 48 adults (23 female and 25 male) and 342 offspring from known mothers of live-bearing bluntnose klipfish. Parentage analysis was performed to interpret mating and reproductive success for both sexes. Metrics quantified were the opportunity for selection (I), the opportunity for sexual selection (Is), absolute (βss), and standardized (β'ss) Bateman gradients and the maximum intensity of precopulatory sexual selection (s'max). Multiple mating by both sexes were revealed by parentage analysis. However, females did not show significant Bateman gradients or a significant maximum intensity of precopulatory sexual selection (s'max), whereas male sexual selection metrics were all significantly greater than 0. These results suggest a polygynandrous mating system for this species. There is an opportunity for sexual selection to act on males but not females in this population, which is evolutionary tied to anisogamy, parental investment, and sex roles.

Inference from single occasion capture experiments using genetic markers

Accurate estimation of the size of animal populations is an important task in ecological science. Recent advances in the field of molecular genetics researches allow the use of genetic data to estimate the size of a population from a single capture occasion rather than repeated occasions as in the usual capture-recapture experiments. Estimating the population size using genetic data also has sometimes led to estimates that differ markedly from each other and also from classical capture-recapture estimates. Here, we develop a closed form estimator that uses genetic information to estimate the size of a population consisting of mothers and daughters, focusing on estimating the number of mothers, using data from a single sample. We demonstrate the estimator is consistent and propose a parametric bootstrap to estimate the standard errors. The estimator is evaluated in a simulation study and applied to real data. We also consider maximum likelihood in this setting and discover problems that preclude its general use.

BATEMANATER: a computer program to estimate and bootstrap mating system variables based on Bateman's principles

Bateman's principles continue to play a major role in the characterization of genetic mating systems in natural populations. The modern manifestations of Bateman's ideas include the opportunity for sexual selection (i.e. I(s) - the variance in relative mating success), the opportunity for selection (i.e. I - the variance in relative reproductive success) and the Bateman gradient (i.e. β(ss) - the slope of the least-squares regression of reproductive success on mating success). These variables serve as the foundation for one convenient approach for the quantification of mating systems. However, their estimation presents at least two challenges, which I address here with a new Windows-based computer software package called BATEMANATER. The first challenge is that confidence intervals for these variables are not easy to calculate. BATEMANATER solves this problem using a bootstrapping approach. The second, more serious, problem is that direct estimates of mating system variables from open populations will typically be biased if some potential progeny or adults are missing from the analysed sample. BATEMANATER addresses this problem using a maximum-likelihood approach to estimate mating system variables from incompletely sampled breeding populations. The current version of BATEMANATER addresses the problem for systems in which progeny can be collected in groups of half- or full-siblings, as would occur when eggs are laid in discrete masses or offspring occur in pregnant females. BATEMANATER has a user-friendly graphical interface and thus represents a new, convenient tool for the characterization and comparison of genetic mating systems.

Genetic fingerprinting proves cross-correlated automatic photo-identification of individuals as highly efficient in large capture-mark-recapture studies

Capture-mark-recapture (CMR) approaches are the backbone of many studies in population ecology to gain insight on the life cycle, migration, habitat use, and demography of target species. The reliable and repeatable recognition of an individual throughout its lifetime is the basic requirement of a CMR study. Although invasive techniques are available to mark individuals permanently, noninvasive methods for individual recognition mainly rest on photographic identification of external body markings, which are unique at the individual level. The re-identification of an individual based on comparing shape patterns of photographs by eye is commonly used. Automated processes for photographic re-identification have been recently established, but their performance in large datasets (i.e., > 1000 individuals) has rarely been tested thoroughly. Here, we evaluated the performance of the program AMPHIDENT, an automatic algorithm to identify individuals on the basis of ventral spot patterns in the great crested newt (Triturus cristatus) versus the genotypic fingerprint of individuals based on highly polymorphic microsatellite loci using GENECAP. Between 2008 and 2010, we captured, sampled and photographed adult newts and calculated for 1648 samples/photographs recapture rates for both approaches. Recapture rates differed slightly with 8.34% for GENECAP and 9.83% for AMPHIDENT. With an estimated rate of 2% false rejections (FRR) and 0.00% false acceptances (FAR), AMPHIDENT proved to be a highly reliable algorithm for CMR studies of large datasets. We conclude that the application of automatic recognition software of individual photographs can be a rather powerful and reliable tool in noninvasive CMR studies for a large number of individuals. Because the cross-correlation of standardized shape patterns is generally applicable to any pattern that provides enough information, this algorithm is capable of becoming a single application with broad use in CMR studies for many species.

Spatial patterns of breeding success of grizzly bears derived from hierarchical multistate models

Conservation programs often manage populations indirectly through the landscapes in which they live. Empirically, linking reproductive success with landscape structure and anthropogenic change is a first step in understanding and managing the spatial mechanisms that affect reproduction, but this link is not sufficiently informed by data. Hierarchical multistate occupancy models can forge these links by estimating spatial patterns of reproductive success across landscapes. To illustrate, we surveyed the occurrence of grizzly bears (Ursus arctos) in the Canadian Rocky Mountains Alberta, Canada. We deployed camera traps for 6 weeks at 54 surveys sites in different types of land cover. We used hierarchical multistate occupancy models to estimate probability of detection, grizzly bear occupancy, and probability of reproductive success at each site. Grizzly bear occupancy varied among cover types and was greater in herbaceous alpine ecotones than in low-elevation wetlands or mid-elevation conifer forests. The conditional probability of reproductive success given grizzly bear occupancy was 30% (SE = 0.14). Grizzly bears with cubs had a higher probability of detection than grizzly bears without cubs, but sites were correctly classified as being occupied by breeding females 49% of the time based on raw data and thus would have been underestimated by half. Repeated surveys and multistate modeling reduced the probability of misclassifying sites occupied by breeders as unoccupied to <2%. The probability of breeding grizzly bear occupancy varied across the landscape. Those patches with highest probabilities of breeding occupancy-herbaceous alpine ecotones-were small and highly dispersed and are projected to shrink as treelines advance due to climate warming. Understanding spatial correlates in breeding distribution is a key requirement for species conservation in the face of climate change and can help identify priorities for landscape management and protection.

Development of a rapid and efficient method for non-lethal DNA sampling and genotyping in scallops

Non-lethal DNA sampling has long appealed to researchers studying population and conservation genetics, as it does not necessitate removing individuals permanently from their natural environment or destroying valuable samples. However, such an approach has not yet been well established in bivalves. In this study, we demonstrate that the gill represents a good source of tissue for non-lethal sampling in scallops. Removal of a few gill filaments caused no noticeable behavioral abnormalities or increased mortality rates in Zhikong scallop (Chlamys farreri) during a three-month period of observation. To facilitate rapid gill-based DNA extraction, six methods (MA-MF) were designed and evaluated, each requiring less than one hour of processing time. The optimal method was identified as MF, in terms of maintaining DNA integrity and genotyping accuracy. Further optimization of MF method by orthogonal experimental design suggested that the utilization of gills could be limited to 2 mg of sample, which is sufficient for performing up to 20,000 PCR reactions. We also demonstrate the excellent cross-species utility of MF in two additional scallop species, Yesso scallop (Patinopecten yessoensis) and bay scallop (Argopecten irradians). Taken together, our study provides a rapid and efficient approach for applying non-lethal DNA sampling in bivalve species, which would serve as a valuable tool for maintaining bivalve populations and conservation genetics, as well as in breeding studies.

Using pedigree reconstruction to estimate population size: genotypes are more than individually unique marks

Estimates of population size are critical for conservation and management, but accurate estimates are difficult to obtain for many species. Noninvasive genetic methods are increasingly used to estimate population size, particularly in elusive species such as large carnivores, which are difficult to count by most other methods. In most such studies, genotypes are treated simply as unique individual identifiers. Here, we develop a new estimator of population size based on pedigree reconstruction. The estimator accounts for individuals that were directly sampled, individuals that were not sampled but whose genotype could be inferred by pedigree reconstruction, and individuals that were not detected by either of these methods. Monte Carlo simulations show that the population estimate is unbiased and precise if sampling is of sufficient intensity and duration. Simulations also identified sampling conditions that can cause the method to overestimate or underestimate true population size; we present and discuss methods to correct these potential biases. The method detected 2–21% more individuals than were directly sampled across a broad range of simulated sampling schemes. Genotypes are more than unique identifiers, and the information about relationships in a set of genotypes can improve estimates of population size.

Reconstruction of paternal genotypes over multiple breeding seasons reveals male green turtles do not breed annually

For species of conservation concern, knowledge of key life-history and demographic components, such as the number and sex ratio of breeding adults, is essential for accurate assessments of population viability. Species with temperature-dependent sex determination can produce heavily biased primary sex ratios, and there is concern that adult sex ratios may be similarly skewed or will become so as a result of climate warming. Prediction and mitigation of such impacts are difficult when life-history information is lacking. In marine turtles, owing to the difficultly in observing males at sea, the breeding interval of males is unknown. It has been suggested that male breeding periodicity may be shorter than that of females, which could help to compensate for generally female-biased sex ratios. Here we outline how the use of molecular-based paternity analysis has allowed us, for the first time, to assess the breeding interval of male marine turtles across multiple breeding seasons. In our study rookery of green turtles (Chelonia mydas), 97% of males were assigned offspring in only one breeding season within the 3-year study period, strongly suggesting that male breeding intervals are frequently longer than 1year at this site. Our results also reveal a sex ratio of breeding adults of at least 1.3 males to each female. This study illustrates the utility of molecular-based parentage inference using reconstruction of parental genotypes as a method for monitoring the number and sex ratio of breeders in species where direct observations or capture are difficult.

Genetic markers reveal high PIT tag retention rates in giant salamanders (Cryptobranchus alleganiensis)

Estimation of population size using mark-recapture (MRR) methods are based on the fundamental assumption that individuals retain their marks throughout the course of study. Passive Integrated Transponder (PIT) tags are useful as a cost effective, reliable marking method in many amphibian and reptile species. Few studies however, use secondary methods to evaluate tag retention rates. Failure to do so can lead to biased population estimates, erroneous conclusions, and thus poor management decisions. Surprisingly, estimates of PIT tag retention are currently lacking for the majority of amphibian species, many of which are experiencing population declines. Herein, we use genetic tagging to assess the retention of PIT tags of the eastern hellbender (Cryptobranchus alleganiensis alleganiensis). We captured and tagged 78 individuals across 35 sites. Recapture rate was 24% and genetic tagging revealed 100% tag retention across all recaptured individuals.

The genetics and genomics of Syngnathidae: pipefishes, seahorses and seadragons

The goal of this review was to provide a historical overview of how molecular techniques have increased the understanding of the ecology and evolution of the family Syngnathidae (pipefishes, seahorses and seadragons). Molecular studies based primarily on mitochondrial DNA markers have proved their worth by elucidating complex phylogenetic relationships within the family. Phylogeographic studies, which have revealed how life-history traits and past climatic events shape geographic distributions and patterns of genetic variation within syngnathid species, also provide interesting case studies for the conservation and management of threatened species. The application of microsatellite DNA markers has opened a floodgate of studies concerned with the breeding biology of these fishes, which are interesting due to their unique reproductive mode of male pregnancy. Research in this area has contributed significantly to the understanding of mating patterns and sexual selection. Molecular markers may also be employed in studies of demography, migration and local breeding population sizes. Genomic studies have identified genes that are probably involved in male pregnancy and promise additional insights into various aspects of syngnathid biology at the level of the gene. Despite these advances, much more remains to be explored. Goals for future research should include: (1) a more inclusive phylogeny to resolve outstanding issues concerning the relationships within the family and higher order taxa, (2) a broader use of molecular studies to aid management and conservation efforts, (3) the inclusion of more genera in comparative behavioural studies and (4) the continued development of genomic resources for syngnathids to facilitate comparative genomic work.

Genotyping errors in a calibrated DNA register: implications for identification of individuals

BACKGROUND

The use of DNA methods for the identification and management of natural resources is gaining importance. In the future, it is likely that DNA registers will play an increasing role in this development. Microsatellite markers have been the primary tool in ecological, medical and forensic genetics for the past two decades. However, these markers are characterized by genotyping errors, and display challenges with calibration between laboratories and genotyping platforms. The Norwegian minke whale DNA register (NMDR) contains individual genetic profiles at ten microsatellite loci for 6737 individuals captured in the period 1997-2008. These analyses have been conducted in four separate laboratories for nearly a decade, and offer a unique opportunity to examine genotyping errors and their consequences in an individual based DNA register. We re-genotyped 240 samples, and, for the first time, applied a mixed regression model to look at potentially confounding effects on genotyping errors.

RESULTS

The average genotyping error rate for the whole dataset was 0.013 per locus and 0.008 per allele. Errors were, however, not evenly distributed. A decreasing trend across time was apparent, along with a strong within-sample correlation, suggesting that error rates heavily depend on sample quality. In addition, some loci were more error prone than others. False allele size constituted 18 of 31 observed errors, and the remaining errors were ten false homozygotes (i.e., the true genotype was a heterozygote) and three false heterozygotes (i.e., the true genotype was a homozygote).

CONCLUSIONS

To our knowledge, this study represents the first investigation of genotyping error rates in a wildlife DNA register, and the first application of mixed models to examine multiple effects of different factors influencing the genotyping quality. It was demonstrated that DNA registers accumulating data over time have the ability to maintain calibration and genotyping consistency, despite analyses being conducted on different genotyping platforms and in different laboratories. Although errors were detected, it is demonstrated that if the re-genotyping of individual samples is possible, these will have a minimal effect on the database's primary purpose, i.e., to perform individual identification.

How many genetic markers to tag an individual? An empirical assessment of false matching rates among close relatives

Genetic identification of individuals is now commonplace, enabling the application of tagging methods to elusive species or species that cannot be tagged by traditional methods. A key aspect is determining the number of loci required to ensure that different individuals have non-matching multi-locus genotypes. Closely related individuals are of particular concern because of elevated matching probabilities caused by their recent co-ancestry. This issue may be addressed by increasing the number of loci to a level where full siblings (the relatedness category with the highest matching probability) are expected to have non-matching multi-locus genotypes. However, increasing the number of loci to meet this "full-sib criterion" greatly increases the laboratory effort, which in turn may increase the genotyping error rate resulting in an upward-biased mark-recapture estimate of abundance as recaptures are missed due to genotyping errors. We assessed the contribution of false matches from close relatives among 425 maternally related humpback whales, each genotyped at 20 microsatellite loci. We observed a very low (0.5-4%) contribution to falsely matching samples from pairs of first-order relatives (i.e., parent and offspring or full siblings). The main contribution to falsely matching individuals from close relatives originated from second-order relatives (e.g., half siblings), which was estimated at 9%. In our study, the total number of observed matches agreed well with expectations based upon the matching probability estimated for unrelated individuals, suggesting that the full-sib criterion is overly conservative, and would have required a 280% relative increase in effort. We suggest that, under most circumstances, the overall contribution to falsely matching samples from close relatives is likely to be low, and hence applying the full-sib criterion is unnecessary. In those cases where close relatives may present a significant issue, such as unrepresentative sampling, we propose three different genotyping strategies requiring only a modest increase in effort, which will greatly reduce the number of false matches due to the presence of related individuals.

Parentage and sibship inference from multilocus genotype data under polygamy

Likelihood methods have been developed to partition individuals in a sample into sibling clusters using genetic marker data without parental information. Most of these methods assume either both sexes are monogamous to infer full sibships only or only one sex is polygamous to infer full sibships and paternal or maternal (but not both) half sibships. We extend our previous method to the more general case of both sexes being polygamous to infer full sibships, paternal half sibships, and maternal half sibships and to the case of a two-generation sample of individuals to infer parentage jointly with sibships. The extension not only expands enormously the scope of application of the method, but also increases its statistical power. The method is implemented for both diploid and haplodiploid species and for codominant and dominant markers, with mutations and genotyping errors accommodated. The performance and robustness of the method are evaluated by analyzing both simulated and empirical data sets. Our method is shown to be much more powerful than pairwise methods in both parentage and sibship assignments because of the more efficient use of marker information. It is little affected by inbreeding in parents and is moderately robust to nonrandom mating and linkage of markers. We also show that individually much less informative markers, such as SNPs or AFLPs, can reach the same power for parentage and sibship inferences as the highly informative marker simple sequence repeats (SSRs), as long as a sufficient number of loci are employed in the analysis.

Detecting female precise natal philopatry in green turtles using assignment methods

It is well established that sea turtles return to natal rookeries to mate and lay their eggs, and that individual females are faithful to particular nesting sites within the rookery. Less certain is whether females are precisely returning to their natal beach. Attempts to demonstrate such precise natal philopatry with genetic data have had mixed success. Here we focused on the green turtles of three nesting sites in the Ascension Island rookery, separated by 5-15 km. Our approach differed from previous work in two key areas. First, we used male microsatellite data (five loci) reconstructed from samples collected from their offspring (N = 17) in addition to data for samples taken directly from females (N = 139). Second, we employed assignment methods in addition to the more traditional F-statistics. No significant genetic structure could be demonstrated with F(ST). However, when average assignment probabilities of females were examined, those for nesting populations in which they were sampled were indeed significantly higher than their probabilities for other populations (Mann-Whitney U-test: P < 0.001). Further evidence was provided by a significant result for the mAI(C) test (P < 0.001), supporting greater natal philopatry for females compared with males. The results suggest that female natal site fidelity was not sufficient for significant genetic differentiation among the nesting populations within the rookery, but detectable with assignment tests.

Review of capture-recapture methods applicable to noninvasive genetic sampling

The use of noninvasive genetic sampling to identify individual animals for capture-recapture studies has become widespread in the past decade. Strong emphasis has been placed on the field protocols and genetic analyses with fruitful results. Little attention has been paid to the capture-recapture application for this specific type of data beyond stating the effects of assumption violations. Here, we review the broad class of capture-recapture methods that are available for use with DNA-based capture-recapture data, noting the array of biologically interesting parameters such as survival, emigration rates, state transition rates and the finite rate of population change that can be estimated from such data. We highlight recent developments in capture-recapture theory specifically designed for noninvasive genetic sampling data.

Population genetic tools for pest management: a review

Population genetic tools have the potential to answer key questions in pest management including quantifying the number of genetically distinct populations represented in an invasion, the number of individuals present, whether populations are expanding or contracting, identifying the origin of invasive individuals, the number of separate introduction events that have occurred and in which order, and the rate that individuals are moving between populations. Genetic methods have only recently gained sufficient resolution to address these questions due to advances in laboratory techniques coupled with an increase in computational power. In combination, these methods may lead to a more comprehensive understanding of the dynamics of invasions. The expansion of the European starling (Sturnus vulgaris) into Western Australia is used as an applied example of how genetic methods can be integrated to provide vital information to improve pest-management strategies. Invasion events also may provide a unique opportunity to test some of these methodologies.

Genetic mating systems and reproductive natural histories of fishes: lessons for ecology and evolution

Fish species have diverse breeding behaviors that make them valuable for testing theories on genetic mating systems and reproductive tactics. Here we review genetic appraisals of paternity and maternity in wild fish populations. Behavioral phenomena quantified by genetic markers in various species include patterns of multiple mating by both sexes; frequent cuckoldry by males and rare cuckoldry by females in nest-tending species; additional routes to surrogate parentage via nest piracy and egg-thievery; egg mimicry by nest-tending males; brood parasitism by helper males in cooperative breeders; clutch mixing in oral brooders; kinship in schooling fry of broadcast spawners; sperm storage by dams in female-pregnant species; and sex-role reversal, polyandry, and strong sexual selection on females in some male-pregnant species. Additional phenomena addressed by genetic parentage analyses in fishes include clustered mutations, filial cannibalism, and local population size. All results are discussed in the context of relevant behavioral and evolutionary theory.

Intensive genetic assessment of the mating system and reproductive success in a semi-closed population of the mottled sculpin, Cottus bairdi

Most genetic surveys of parentage in nature sample only a small fraction of the breeding population. Here we apply microsatellite markers to deduce the genetic mating system and assess the reproductive success of females and males in an extensively collected, semi-closed stream population of the mottled sculpin fish, Cottus bairdi. In this species, males guard nest rocks where females deposit the eggs for fertilization. The potential exists for both males and females to mate with multiple partners and for males to provide parental care to genetically unrelated offspring. Four hundred and fifty-five adults and subadults, as well as 1,259 offspring from 23 nests, were genotyped at five polymorphic microsatellite loci. Multilocus maternal genotypes, deduced via genetic analyses of embryos, were reconstructed for more than 90% of the analysed nests, thus allowing both male and female reproductive success to be estimated accurately. There was no genetic evidence for cuckoldry, but one nest probably represents a takeover event. Successful males spawned with a mean of 2.8 partners, whereas each female apparently deposited her entire clutch of eggs in a single nest (mean fecundity = 66 eggs/female). On average, genetically deduced sires and dams were captured 1.6 and 9.3 metres from their respective nests, indicating little movement by breeders during the spawning season. Based on a 'genetic mark-recapture' estimate, the total number of potentially breeding adults (c. 570) was an order-of-magnitude larger than genetically based estimates of the effective number of breeders (c. 54). In addition, significantly fewer eggs per female were deposited in single than in multidam nests. Not only were perceived high-quality males spawning with multiple partners, but they were receiving more eggs from each female.