Introduction

The last recorded presence of the Eurasian otter (Lutra lutra) in the Netherlands dates from 1989: it was a dead individual trapped in a fyke. Otter numbers had plummeted from a flourishing population in the 1930s (van Wijngaarden and van de Peppel 1970). This rapid extinction of the Dutch otter population is typical of what happened in many industrialised and densely populated Western-European countries. The factors contributing to the decline of this aquatic top predator were environmental pollution (Roos et al. 2001), an increase in road kill due to habitat fragmentation (Hauer et al. 2002a; Kruuk and Conroy 1991; Sommer et al. 2005) and incremented pressure from fisheries with fyke nets.

After the extinction of the otter in the Netherlands, there was a call for the return of this iconic species to the Dutch wetlands. In the mid-1990s measures were taken to restore natural otter habitat, improve water quality, reconstruct interconnecting corridors and build road underpasses in a specific area in the north of the Netherlands. It was made mandatory to use stop grids on fyke nets in this area, to reduce the risk of young otters drowning.

In 2002 a reintroduction programme was started, since it seemed unlikely that otters would recolonise spontaneously from the nearest populations in Germany in the short term. Reintroduction projects attempt to re-establish species within their historical ranges by releasing wild or captive-bred individuals following extirpation or extinction in the wild (Seddon et al. 2007). They are often carried out to fulfill a biodiversity preservation or restoration objective. According to the IUCN Reintroduction Guidelines (1998), the ultimate goal of reintroduction projects should be “the establishment of a self-sustaining population that requires minimal long-term management”. Several authors have pointed out that in order to achieve this goal, projects should be accompanied by “focused monitoring” to identify the possible cause of success or failure (Nichols and Williams 2006; Seddon et al. 2007; Armstrong and Seddon 2008) and that the resulting information should be used to improve future project designs.

Currently, most inferences in reintroduction biology are made by induction, gained from post hoc interpretation of monitoring results or by exploratory comparative analyses (Seddon et al. 2007). The main information needed for a long-term perspective on reintroduced populations is on the vital statistics, e.g. mortality and recruitment. This entails extensive monitoring, which is an obstacle for many projects because of its expense (Morell 2008). Yet from the many anecdotal and contradictory results (Seddon et al. 2007), it is clear that post-release monitoring is badly needed in reintroduction efforts.

Genetic monitoring, i.e. quantifying temporal changes in population metrics or other population data generated using molecular markers, is becoming increasingly important for monitoring the consequences of anthropogenic change on wild species (Allendorf and Luikart 2006; Schwartz et al. 2007; Wayne and Morin 2004). Molecular markers can be used either as a diagnostic tool for individual identification and traditional population monitoring, or for estimating changes in population genetic parameters, providing information on genetic diversity, inbreeding, effective population size or migration (Schwartz et al. 2007). Genetic data is invaluable for ecologists and wildlife managers, especially when it is combined with behavioural, demographic, or spatial information (DeYoung and Honeycutt 2005). In contrast to traditional techniques like telemetry, genetic monitoring offers advanced possibilities for the long-term monitoring of population structure. The development of noninvasive genetic procedures, i.e. extracting DNA from biological samples such as hair, faeces, or urine without handling, capturing or even observing the animals (Kohn and Wayne 1997; Taberlet et al. 1999; Piggott and Taylor 2003; Waits and Paetkau 2005) yields details on elusive and nocturnal animals like bears (Kendall et al. 2009), coyotes (Kohn et al. 1999; Prugh et al. 2005), otters (Ferrando et al. 2008; Hung et al. 2004), lynx (Palomares et al. 2002) and wolves (Creel et al. 2003). Combining genetic monitoring and noninvasive sampling into noninvasive genetic monitoring offers an excellent supplementary technique for studying the structure of reintroduced populations, as this enables time-dependent processes such as trends in inbreeding and recruitment to be monitored, especially if DNA is sampled from the outset of the reintroduction.

Despite being very promising, the technique, and thereby the quality of the results, is constrained by the quantity and quality of DNA extracted from the biological samples (Miller et al. 2002; Taberlet et al. 1996). Genotyping from secondary material is prone to several problems, e.g. allelic dropout caused by the random amplification of only one of two alleles at a heterozygous locus due to the scarcity of template DNA. Another type of error, a “false allele”, is an artefact generated during the amplification process, which often reveals a spurious or third allele. To ensure that estimates of abundance and individual identification are reliable, these errors need to be detected, e.g. by repeating the DNA amplifications independently several times, in order to obtain trustworthy consensus genotypes (Taberlet et al. 1996).

Though otters are elusive and nocturnal animals and are therefore hard to trace (Kruuk 2006; Mason and Macdonald 1986) they leave scent markers of their home range at prominent sites of aquatic habitats, such as bridges, logs, and sandy beaches. In our survey, which was started concomitantly with the reintroduction of the first otters into the Netherlands, we used DNA isolated from freshly deposited spraints (faeces) and anal secretions (jellies) to assess (a) the presence and distribution of the founders, (b) the recruitment of new individuals (offspring) into the population, (c) to determine the ancestry of these offspring and infer the social and genetic structure of the population, and (d) to estimate the abundance of otters in the study area during successive years. We did so in the belief that only noninvasive genetic monitoring could give us the answers to the probing questions surrounding the reintroduction of elusive animals (see Armstrong and Seddon 2008; Morell 2008; Schwartz et al. 2007). Using DNA as a tag to identify individuals we were able to build a detailed life table and a pedigree of the otter population, to evaluate the first phase of the reintroduction programme and to suggest how the reintroduction programme should be continued.

Materials and methods

Study area

The study area comprised lowland peat marshes in the north of the Netherlands, ca. 200 km2 in extent and located at 52°38′–52°50′ N, 5°53′–6°09′ E. The area consists of a mosaic of peat grasslands, reed beds (Phragmites australis) and swamp woodland, mainly composed of Willow (Salix spp.), Black alder (Alnus glutinosa) and Birch species (Betula spp.) intersected by small ditches and pools. Most of the surrounding landscape is intensively farmed. The mean annual temperature is 9.5°C and on average there are 10 frost days per year. The area can be divided into three zones: the Weerribben and the Rottige Meenthe, which comprise peat grassland, woodland and ditches, and the Wieden, where larger lakes are also present (Fig. 1).

Fig. 1
figure 1

Geographical location of the study area. The black border surrounds the area where 30 otters were released until 1 April 2008. The expanding population was monitored using noninvasive genetic sampling

Origin and release

The otters used in the Dutch reintroduction project had either been captured in the wild in Belarus, Latvia, or Poland, or originated from captivity or rehabilitation programmes in Sweden (Finnish origin), the Czech Republic, or Germany. The otters caught in the wild came from areas with high otter densities. Between July 2002 and April 2008, 30 otters were sequentially released in the study area (Table 1). The first set of otters was released in the Weerribben, the central location, and in subsequent years otters were gradually added to the Wieden and Rottige Meenthe areas (Table 1). At the time of release the animals varied in age between 1 and 5 years old. Before their release, tissue and blood samples were taken for DNA fingerprinting and all the otters were tagged with a transponder and were fitted with a radio transmitter, implanted intraperitoneally, for monitoring initial post-release movements and survival. The implantations were carried out by a veterinarian at Burgers Zoo (Arnhem, The Netherlands), where the otters were kept in captivity for 3–18 days and observed before being released. The study was conducted in accordance with Dutch legislation on the protection and welfare of vertebrate animals used for experimental and other scientific purposes.

Table 1 Number of otters released until 1 April 2008

Sample collection

Sprainting behaviour often shows seasonality and there is evidence that winter is the best period to collect otter spraints (Kruuk 1992) and that the microsatellite DNA analysis of otter faeces is most successful in cold months, when spraints are collected in the early morning (Hájková et al. 2006). During summer it is difficult to collect fresh spraints in our area because the tall grass and reeds obstruct visibility. Therefore, our surveys were carried out in the winter half year (October to the end of March) of consecutive years from 2002 to 2008.

Each winter period we checked the whole release area for otter activity (spraints, spoor, landing sites, trails). In November and December each year we did a first survey of the whole release area covering about 10–15 km2 per day by foot, bike and boat. From January to March we conducted a second survey. We tried to cover the whole release area at least twice and revisited promising sites suggested by the field managers. On average we spent 48 days in the field each winter, searching for otter marks. We tried to revisit each location from which we had collected a spraint, with the aim of obtaining at least 3 fresh spraints from that location. The constraints were the accessibility of the terrain and the availability of waterways. We used GPS to record the locations of spraints and thus ascertain the spatial organisation of the population. If possible, the area was visited the day before collection and old spraints were marked to increase the chance of identifying fresh spraints the next day. When this strategy could not be employed, we collected all the spraints that seemed to be fresh. To minimise DNA degradation, spraints were collected in the morning. The samples were immediately put into 10 ml plastic phials containing 99% ethanol and taken to the lab, where they were stored at −20°C until DNA extraction and analysis. A total of 1,265 spraints were collected for genetic analysis.

We conducted a vigorous publicity campaign to encourage people to contact us when dead otters were found, so that we would have a as complete picture of the population as possible. Dead otters were delivered to Alterra throughout the year and subjected to post mortem analysis to determine the most likely cause of death. Tissue samples from the cadavers were stored in phials containing 99% ethanol. Most of the otters brought in were road kills from outside the release area.

Rationale of the monitoring design

Since our otter population is small (n = 30 founders) and isolated, and immigration from Germany seems highly unlikely, it is effectively a closed population. At the start of the project we constructed a reference database of the genetic profiles of the founders, so that in subsequent years we could infer successful mating and recruitment from new genetic profiles obtained from spraints or dead animals. We updated the database of genetic profiles yearly. Since all potential parents were known, we applied complete exclusion as our method of parentage analysis (Blouin 2003; Jones and Ardren 2003). A spraint was classified as from an offspring if (a) the genetic profile did not match existing profiles from previous years, and (b) the profile could unequivocally be assigned to a known male and female. The offspring–mother–father combination was subsequently checked by comparing the distribution of GPS-coordinate observations of the animals involved.

DNA extraction

Faecal DNA was extracted using a modification of the hexadecyltrimethylammonium bromide (CTAB)-based extraction (Parsons et al. 1999; Hung et al. 2004). This entailed removing a spraint from the phial of ethanol and briefly putting it on filter paper to remove most of the ethanol. Next, a raisin-sized piece was put in a 2 ml Eppendorf tube together with 1 ml of CTAB buffer (100 mM Tris–HCl pH 8, 20 mM EDTA, 1.4 M NaCl, 2% CTAB) and the spraint was homogenised by stirring with a small rod. After adding additional CTAB buffer to bring the total volume to 2 ml, the mixture was briefly vortexed and left on a shaker for 15 min. This mixture was centrifuged for 5 min and 1.5 ml of the supernatant was transferred into a new tube together with 0.5 ml chloroform. After two rounds of extraction, DNA was precipitated by adding 0.67 ml isopropanol to 1 ml of the cleared suspension. The resulting pellet was resuspended in 0.18 ml of ATL buffer. Spraint pellets and tissue from released and dead individuals were further processed following the protocol of the DNeasy Blood and Tissue Kit (Qiagen) for DNA purification.

Genetic typing

During the first three winters seven microsatellites were sufficient for individual typing and parentage assessments: Lut701, Lut715, Lut717, Lut733, Lut818, Lut832 and Lut833 (Dallas and Piertney 1998). Subsequently because of the loss of released individuals, the occurrence of offspring and the increasing relatedness among individuals, we had to gradually increase the number of microsatellite loci: OT04, OT05, OT07, OT14, OT17, OT19 and OT22 (Huang et al. 2005) and RI18 (Beheler et al. 2005). We only used tetranucleotide microsatellite loci to reduce the occurrence of stutter bands and ambiguity in scoring that often happens with dinucleotide loci.

PCR reactions were performed in a total volume of 10 μl containing 0.3 Units of Taq (Invitrogen Taq DNA polymerase (18038-034), amounts of PCR buffer and W-1 according to the Invitrogen protocol, 130 nM of each primer, 200 μM of each dNTP, 4.25 mM MgCl2 and 320 μg/ml BSA. Forward primers were labelled with either an IRD-700 or an IRD-800. The PCR programme used was 95°C/3 min and (90°C/30 s, T a/30 s, 72°C/1 min.) × 39 cycles. For most primers T a was 60°C, except for locus Lut715 (T a = 58°C), loci Lut733, Lut782 and Lut818 (T a = 59°C), locus Lut717 (T a = 61°C), and locus OT07 (T a = 62°C). The same protocol was used for tissue extracts, except that a dilution factor of 10 was applied and 2 μl of this diluted extract was used. For sex identification we used the DBY7Ggu primer following the protocol of Hedmark et al. (2004).

PCR products of microsatellite loci and sexual typing were genotyped on a 6.5% polyacrylamide gel containing 7 M Urea and 1× TBE on a Li-Cor 4300 platform.

Data analysis

We did not use the same set of microsatellite loci every year because initially the seven loci of Dallas and Piertney (1998) had sufficient power to distinguish individuals and assess parentage. In 2006/07 these seven loci were still sufficient for identifying individuals, but not for assessing parentage. After optimising and adding the second set of loci we made two new sets of loci: (i) a set for distinguishing individuals (Lut715, Lut717, Lut733, Lut833, OT07, OT19 and OT22); and (ii) a set with the remainder of the loci, which we used solely to confirm parentage assessment or in cases of doubtful identification of an individual because of the failure of a locus. Budget constraints prevented us from using all loci. The criterion for compiling our first set of loci was a P (ID)sib of < 0.01. The probability of identity, P (ID), is the probability that two individuals drawn at random from a population have the same genotype at multiple loci (Creel et al. 2003; Taberlet and Luikart 1999; Waits et al. 2001). It is considered to be the most common statistic used to quantify the power of molecular markers in distinguishing two individuals. The P (ID)sib, the P (ID) among a population composed solely of siblings, gives an upper limit to the possible range of P (ID) in a population. At the beginning of the project, in 2002/03, the P (ID) among the released animals was 1.9 × 10−7 and the P (ID)sib was 1.9 × 10−3. Because of the changing population composition, e.g. loss of released individuals, presence of offspring and increased relatedness, in 2006/07 P (ID) was 4.1 × 10−6 and P (ID)sib was 4.4 × 10−3. During the 2007/08 season the P (ID)sib of our first set of loci was 2.1 × 10−3, but when all 15 loci were used it increased to 1.4 × 10−5.

To reduce the chance of mistyping, we applied a modified multiple tube approach (Gagneux et al. 1997; Hung et al. 2004; Taberlet et al. 1996). The constraint on the modified approach was that with our current protocol we could only run ca. 50 PCRs from one faecal extract and therefore had to adjust the number of replicates when using all 15 microsatellite loci. Our approach was as follows: (1) Each sample was amplified three times for locus LUT715. This locus was chosen because of good, repeatable results in previous experiments. (2) The sample was discarded from the subsequent analyses if there were less than three PCR products. In case of three PCR products a sample was still discarded if, after scoring the results, it resulted in three different typings. (3) Selected samples were amplified three times for the remaining loci from the first set. (4) Three independent typings with the same single allele at a locus confirmed a homozygote. Three independent typings with the same two alleles at a locus confirmed a heterozygote. Samples with two typings of a heterozygotes and one homozygous typing, were scored as heterozygous with the two alleles appearing in these typings. (5) For loci that were typed twice as homozygous and once as heterozygous, or for loci that were scored as homozygous for different alleles, three additional independent typings were performed. When among the six typings an allele was recorded at least twice, the sample was accepted as heterozygote. If an allele appeared only once among the six typings the sample was accepted as a possible homozygote. (6) Those samples that could not be appropriately typed up after six typing attempts were discarded. (7) When the genotypes of two samples were the same at six loci and the only mismatch at the seventh locus may have been due to allelic dropout, we considered the two samples to be the same multilocus genotype if they came from geographical locations close to each other. As safeguard, these samples were also typed for the second set of loci. If the only mismatch was an unambiguous different typing at the seventh locus, the samples were always typed for the second set of loci. (8) The consensus genotypes obtained were compared with the reference database and, if possible, assigned to known individuals. (9) In case of new profiles we tried to assess parentage and completed the profile for the remaining loci. Final parentage was assessed on 15 loci, with complete exclusion as the criterion (Blouin 2003).

Since we had complete genetic profiles for the otters released, we calculated allelic dropouts and false alleles as deviations from the expected profile after assigning spraint samples to known individuals. The allelic dropout for a locus was computed from the number of homozygotes typed for the heterozygous individuals, divided by the total number of heterozygous samples. Allelic dropout occurred in 15.2% of the heterozygous samples and varied among loci (range: 8.1% (LUT833)–22.4% (LUT717)). When using three independent amplifications, the probability of obtaining false homozygotes was 0.0018, using the equation P = K × (K/2)n − 1 (Gagneux et al. 1997), where K is the observed frequency of false homozygotes averaged over all individuals and loci and n is the number of repeated amplifications. Strictly speaking, false alleles are undefined amplification products that show up as a spurious or third allele. These occurred in only 2.1% of the PCR reactions. More broadly speaking, a heterozygous typing of a homozygote individual could also be considered as a false allele. The observed frequency of false heterozygotes was 14.9% and varied among loci (range: 4.0% (LUT733)–26.2% (LUT717)). As expected, these values are comparable to the allelic dropout results.

Results

Spraint surveys

Of a total of 1,265 spraint samples, 582 were successfully genotyped, resulting in an overall success rate of 46% (Table 2). Our success rate gradually improved during the years as a result of refining the laboratory protocols. During the 2007/08 winter we achieved a success rate of 59% (Table 2). The number of spraints per individual per census varied from 1 to 27, indicating considerable capture heterogeneity among individuals.

Table 2 Summarised genetic data on otter spraints collected between October 2002 and April 2008 (a) and otters recovered dead from the same period and from which no spraint observations were known (b)

Of the 30 individuals released, only 15 (50%) were re-identified from spraints (Table 2). We identified a total of 54 offspring during the period 2002–2008: 23 females and 31 males (Tables 2, 3). Of these offspring, 41 (76%) were identified from spraints and 13 (24%) as dead animals without previous spraint signature (Table 2). Among the offspring, the likelihood of being detected through spraints differed significantly between the sexes: 91% for the females and 65% for the males (Χ2 = 5.18, P = 0.023; Table 2).

Table 3 (a) Reproductive patterns and number of offspring observed during the period July 2002–April 2008. (b) The number of new offspring identified during the successive winter periods, grouped according to the origin of the mother and the origin of the father. The parent otter could be a founding animal (founder) or an offspring from previous years (native-born); (c) the number of reproducing females and males during the successive winter periods. Individuals are classified according to their origin: either a founding animal (founder) or an individual born in the release area (native-born). Also shown is the range in the number of matings, i.e. different females, per male

Mating and recruitment

We inferred a total of 30 matings and 54 offspring (Tables 3, 4). In the consecutive winter periods, an increasing number of females were involved in mating; offspring females started to mate after 2 years (Table 3a, c). The number of males involved in reproduction was much lower and initially declined, but it increased during winter 2006/07 when the first generation of offspring males became sexually active (Table 3c). However, offspring males took twice as long as offspring females to become involved in population recruitment (Table 3a, c). The number of females a male mated with varied greatly: from one to seven (Table 3c). During winter 2005/06 only a single male (A08) contributed to the reproduction by mating with seven different females (Tables 3a, 4).

Table 4 Summarised data on female and male reproductive success during the period July 2002–April 2008

The number of offspring identified increased steadily during consecutive winters (Table 3b). The winter of 2007/08 saw the first third-generation offspring resulting from matings between second-generation males and females (the offspring of the founder, i.e. first generation, population) the previous year (Table 3b, c).

Reproductive success

Throughout the observation period we found offspring of 16 founders (9 females and 7 males: Table 3). These founders contributed to reproduction predominantly during the early years; later they were replaced by their descendants, especially during the winter of 2006/07 (Table 3b, c). Despite the involvement of many individuals in the mating process (Table 4), reproductive success was very skewed, especially for the males (Fig. 2). When considering only those individuals that could have reached maturity by winter 2007/08, 34% of the males observed in our population produced offspring; for the females this figure was more balanced, but still only 51% of the females became involved in successful reproduction (Fig. 2). The average number of offspring per female was 1.67 ± 4.92 (mean ± variance; n = 33), for the males this was 1.90 ± 24.71 (n = 29). The high variances are indicative of the skewed distribution of reproductive success.

Fig. 2
figure 2

Frequency distribution of the reproductive success of females and males during the period 2002–2008

A large number of offspring (25 out of 54) shared the same father (A08; Table 4). This male was solely responsible for the 15 offspring observed during winter 2006/07 (Table 3b, c). After the disappearance of A08 two other males (NB11 and NB15) took over his dominant position. The more dominant one of the two was NB15, a son of A08 (Table 4), so A08 still continued to put his stamp on the population.

Fate of the individuals

By 1 April 2008 though we still had records on the presence of 30% (6 out of 20) of the females released, we no longer had evidence that any of the ten males released were still alive. The comparable figures for the offspring females and males were 78 and 29%, respectively (Table 5). The data shown in the alive and dead columns of Table 5 show that offspring males had a significantly higher mortality than offspring females (Χ2 = 12.41, P < 0.001; Table 5). The majority of the dead otters delivered to Alterra were subadult males that had been found outside the release area beyond a radius of ca. 50 km. In most instances they were traffic victims.

Table 5 Fate of the individuals since the start of the reintroduction programme in 2002, grouped by origin (founder or native-born) and sex

The fate of 43% of the otters released is unknown, the comparable figure for the observed offspring is only 12% (Table 5).

State of the population

Counting the identified genetic profiles and dead animals reveals a steady growth in the population, with individuals designated as “native-born” taking over (Fig. 3). Initially, the population hardly grew and it took 4 years before the actual population size exceeded the number of animals released. The total number of individuals identified during winter 2007/08 was 47.

Fig. 3
figure 3

The development and composition of the reintroduced Dutch otter population: the number of animals identified and their origin, either a founding animal or offspring born in the release area (native-born) during subsequent winters. The dotted line indicates the cumulative number of animals released during the period 2002–2008

Discussion

Noninvasive genetic methodology

After the initial euphoria concerning the use of noninvasively collected biological material for genetic research, the shortcomings arising from the small amounts of DNA, or the poor quality of the DNA or of the extract became clear, and a more cautious approach was advocated (Schwartz et al. 1999; Taberlet et al. 1999). Since then, much emphasis has been on the development and optimisation of lab and statistical protocols to safeguard against unjustified interpretations of the results (Bonin et al. 2004; Mills et al. 2000; Waits and Paetkau 2005). If mistakes occur in the genetic profiles due to genotyping errors caused by allelic dropouts, there is a risk of overestimating population size (Mills et al. 2000). The risk can be reduced by including fewer loci in the genotype or allowing one or two mismatches among loci, thereby reducing the probability of creating spurious genotypes. This, however, reduces the potential to distinguish closely related genotypes. In contrast to most other noninvasive genetic studies, we had the advantage of working with a closed population originating from a limited number of founders that had been genotyped beforehand. We could therefore even permit uncertain or wrong scoring when interpreting our data. The fact that our P (ID)sib values ranged from 1.9 × 10−3 (2002/03) to 2.1 × 10−3 (based on 7 loci) or 1.4 × 10−5 (based on 15 loci) (2007/08) indicates that the number of loci we used was sufficient to differentiate between individuals, juveniles and even siblings of the second and third generations. We followed the so-called matching approach (Creel et al. 2003), allowing a maximum of one mismatch to be scored as an identical genotype when we did the initial typing with seven microsatellites. Subsequently, we located the geographical position for the samples for which we allowed a mismatch and checked whether it fell within the area of the assigned known individual or assigned parents in case of native-born individuals. When there was a big discrepancy between the two geographical locations, the sample was processed for another set of microsatellite loci, as control. We are therefore confident that we did not overestimate population abundance through spraint genotyping. On the contrary, we most likely underestimated the current population size because we only detected 76% of the offspring through spraint genotyping (Table 2a). The remaining 24%, mostly subadult males, only became apparent when discovered dead and delivered to Alterra for post mortem analysis (Table 2b). Based on the placental scars of dead adult females delivered to Alterra we estimated the number of offspring per female to be 2.2 (n = 4). The latter value is close to the 2.3 newborn per female reported by Hauer et al. (2002b) for a German population. Since we inferred 30 matings (Table 4), this would give rise to 66–69 native-born offspring while we observed only 54 offspring (Table 4). Therefore, we suspect that we have underestimated the current population size by at least 20%.

Simulations with the CAPWIRE program (Miller et al. 2005), that takes into account capture heterogeneity among individuals based on spraint genotyping, indicated that the expected number of observed genotypes during winter 2007/08 was 56 (confidence limits 49–67); on the basis of the genetic profiles we identified 47 genotypes (Fig. 3). This also suggests that we are underestimating the current density.

Reproductive success and dominance hierarchies

Faecal genotyping has been used frequently in recent years for studying otter distribution and abundance (Arrendal et al. 2007; Dallas et al. 2003; Ferrando et al. 2008; Hájková et al. 2008; Hung et al. 2004; Lanszki et al. 2008). Here we also report on the recruitment of offspring into the population. We were able to detect 54 offspring, and deduced that founder females were reproducing successfully within 4 months of their release. We also found evidence for second-generation breeding within 4 years (males NB11 and NB15 mating with female offspring; Table 4). The involvement of native-born individuals in reproduction is one of the criteria of the short-term success of reintroduction efforts (IUCN 1998; Morell 2008; Seddon et al. 2007).

Descriptive field studies suggest that the social system of Eurasian otters consists of female territories and independent, larger male territories (Erlinge 1968; Kruuk 2006). Our genetic data on the reproductive success of males and females confirm these observations, since (i) in contrast to the females, only a limited number of males contributed to reproduction (Table 3), (ii) among these males there was a hierarchy of dominance, with a few highly productive males (A08 and NB15; Table 4; Fig. 2), and (iii) males mated with several females during one season (Table 3c). As a consequence, male reproductive success is highly skewed and will cause increased relatedness in our reintroduced population.

The dominance of a few males is also reflected in the distribution patterns of individual animals (data not shown) and the number of spraints collected per individual. Males with large territories were represented more frequently in our spraint surveys—a finding that agrees with previous studies (see Hájková et al. 2008). At the start, in 2002, A12 was the dominant male; he was succeeded by A08 in the Weerribben (winter 2005/06) and by A18 in the Wieden (winter 2004/05) (Table 2a). After A18 died, A08 took control over the whole release area during winter 2005/06, having relinquished the Weerribben in 2006/07 to his son NB15 who then became the most dominant male (Table 2a). The number of dominant males dividing up the release area is much lower than expected. Currently only 2–3 adult males occupy the area, whereas it was expected that the area was big enough for 5–6 males.

The dominant behaviour of adult males caused the exodus of juvenile and subadult males—a pattern frequently observed in otters (Kruuk and Moorehouse 1991) and other carnivores (Creel and MacDonald 1995). Juvenile males were usually observed at the borders of the release area, occasionally roaming the area and waiting for an opportunity to enter. When A08 disappeared from the Wieden his place was taken over by NB30 (winter 2007/08; Table 2a). However, while roaming through his newly acquired territory this male was killed in a road traffic accident before he had the chance to reproduce. To our knowledge we are the first to document the social behaviour of the otter in such detail and this was only feasible because of the detailed results obtained through noninvasive genetic analyses.

Effective population size

A consequence of the social system and the skewed reproductive success is that the average number of individuals (N) counted in a population does not necessarily reflect the true population size. To estimate the effective number of breeding individuals Ne, i.e. the number of individuals in an ideal population (random mating, equal sex contributions to the next generation) that gives rise to the same observed variation, corrections must be made to the number counted for the juvenile/adult ratio, the variation in numbers in time, the sex ratio and the non-random distribution of family size (Falconer and MacKay 1996). Ne is the most reliable indicator of the effective number of breeding individuals in a population. When correcting for variation in family size Ne is estimated by Ne ≈ 8 N/(V f + V m + 4), where N is the observed number of adults and V f and V m are the observed variances in family size for females and males. In an ideal population both V f and V m are equal to 2, thereby resembling a poisson distribution of reproductive success with a mean of 2 offspring per mating (Falconer and MacKay 1996). Using the variance in reproductive success as a measure for variance in family size, Ne was estimated as 8 N/(4 + 4.92 + 24.71); thus, about 24% of the observed number. Since the main cause of this reduction is the high variance in male reproductive success (Fig. 2), which is a consequence of the social system of otters, it is difficult to take measures to reduce the large variance in male reproductive success. The real effective number will be even lower after correcting for the other factors. Thus, despite the increase observed in the number of animals in the population, the population is still vulnerable because the effective population is small.

Our findings are in line with the observations on other species. Frankham (1995) estimated the Ne/N ratio in natural populations of animals to be 0.11 and argued that the main determinants of the reduction were variation in numbers between years and the variance in effective family size. The variation in family size reduced effective population sizes to an average of 54% of the census sizes. Heywood (1986), though working with plant populations, estimated the Ne/N ratio to be 0.34 and, interestingly, demonstrated that the main contributor to this reduction was the large variance in seed production among individual plants.

Sex-biased dispersal, mortality and the modification of the release plan

The genetic results revealed that mortality among young males was much higher than their female siblings (Table 5). Road traffic accidents turned out to be responsible for 87% of the otters found dead and most of the victims were males found outside the release area (within a radius of ca. 50 km around the release area). Apparently, juvenile females remain in the release area, while juvenile males leave the area. This pattern is confirmed by the geographical distribution of spraint observations (data not shown): juvenile females stay close to the territories of their mothers and have a high incidence of breeding next to the mother’s territory in subsequent years. The juvenile males, on the other hand, are forced out of the area by the dominant males (cf. Kruuk and Moorehouse 1991).

Our results show that traffic is the main threat to the reintroduced otters, as has been found in other studies throughout Europe (Hauer et al. 2002a; Kruuk and Conroy 1991; Sommer et al. 2005). This emphasises the need for increased habitat connectivity and should receive highest priority in the near future, to ensure that otters are encouraged to migrate to suitable areas outside the study area, as otherwise a new isolated otter population will have been created.

Because of the sex-biased dispersal it was decided to change the release strategy during the course of the project. Since the juvenile females stayed in the release area and the juvenile males were forced out of the area, it was decided that new releases in the adjacent satellite areas like the Oldematen and the Lindevallei (Fig. 1) should be only of females, as this might discourage the males from travelling further away from the release area.

The need for post-release monitoring in reintroduction studies

As pointed out by Seddon et al. (2007), most of the information on reintroduction programmes consists of descriptive accounts and retrospective analyses. The research questions addressed have largely been driven by the monitoring data available, instead of the monitoring being driven by the questions. Before starting the monitoring we discussed which information would be essential to field managers and how this could best be obtained. We decided to use noninvasive genetic monitoring as our main strategy, and to use traditional telemetry only to detect the initial establishment of the founders. Typing all founding animals before release provided us with a base reference population from the outset. This individual-oriented approach is feasible, especially in reintroduction studies which are usually based on a limited number of individuals; it will provide information on processes related to small populations. During the first 3 years of the reintroduction project the project had a bad press because the general public did not see the otters and the media reported traffic incidents and the low number of identified individuals (Fig. 3). It was not until we had demonstrated that (i) otters were reproducing (Table 3), (ii) the number of individual otters identified exceeded the number of otters released (Fig. 3), and (iii) we were on course for a new Dutch population (Fig. 3), that attitudes changed. The otter population did not start to grow until after the offspring of the founders took part in the reproduction (Fig. 3). Only genetics could reveal these details. The lesson is that, when reintroducing animals, one should be prepared for a slow start and accept losses due to mortality. These processes are inevitable.

The future of the Dutch otter population

Despite the initial success in achieving a growing population and detecting mating between native-born individuals, both of which are short-term criteria for successful reintroductions (Morell 2008; Seddon et al. 2007), the future of the Dutch otter population remains uncertain. We have created a new, small and isolated population with a low effective population size. This population is currently approaching the maximum carrying capacity of the release area. The original plan aimed at a connection with a second settled population in nearby wetlands. This has not yet been achieved, and consequently many animals are killed in traffic incidents when they move away from the current population. Moreover, the otters have their own strategy: only a few males dominate the reproductive process, thereby lowering the effective population size and increasing relatedness and the occurrence of inbreeding. Whether inbreeding will have an effect on the current population is yet unknown, since no data are available on the occurrence of inbreeding depression in otters. However, data from other large mammal species (Laikre 1999; Liberg et al. 2005) indicate that inbreeding effects can be substantial. To counteract possible negative effects of inbreeding, the preferred option would be to establish connection with other or newly founded nearby populations. If this is not feasible, the only option might be to regularly introduce a new unrelated male into the current population.

Conclusions

Noninvasive genetic monitoring has given us insight into the hidden life of otters and has proved to be a valuable tool for conservation purposes. Traditional monitoring techniques would not have provided us with an as detailed picture of the reintroduced otter population. The results on the differences in social and territorial behaviour between males and females forced us to reconsider our initial release plan and also indicate that however well the introduction methods are prepared for by developing a scientific and socio-economic strategy, there will be surprises because the key players have a strategy of their own (cf. Morell 2008). Our studies confirm the need for intensive post-release monitoring of reintroduced species (cf. Armstrong and Seddon 2008).