Genetic Consequences of Multiple Translocations of the Banded Hare-Wallaby in Western Australia

Many Australian mammal species now only occur on islands and fenced mainland havens free from invasive predators. The range of one species, the banded hare-wallaby (Lagostrophus fasciatus), had contracted to two offshore islands in Western Australia. To improve survival, four conservation translocations have been attempted with mixed success, and all occurred in the absence of genetic information. Here, we genotyped seven polymorphic microsatellite markers in two source (Bernier Island and Dorre Island), two historic captive, and two translocated L. fasciatus populations to determine the impact of multiple translocations on genetic diversity. Subsequently, we used population viability analysis (PVA) and gene retention modelling to determine scenarios that will maximise demographic resilience and genetic richness of two new populations that are currently being established. One translocated population (Wadderin) has undergone a genetic bottleneck and lost 8.1% of its source population’s allelic diversity, while the other (Faure Island) may be inbred. We show that founder number is a key parameter when establishing new L. fasciatus populations and 100 founders should lead to high survival probabilities. Our modelling predicts that during periodic droughts, the recovery of source populations will be slower post-harvest, while 75% more animals—about 60 individuals—are required to retain adequate allelic diversity in the translocated population. Our approach demonstrates how genetic data coupled with simulations of stochastic environmental events can address central questions in translocation programmes.


Introduction
Translocation, the anthropogenic movement of a group of organisms from one location to another, is an increasingly necessary tool for conservation management [1]. Conservation translocations take many forms (reintroductions, reinforcements, genetic rescue, assisted colonisation) and in all cases should be carefully planned, as they are costly and their success is hard to predict [2][3][4][5][6][7]. For example, in Since 1974, four translocations of L. fasciatus have been attempted, and two captive breeding colonies have been established but since discontinued (Figure 1d). The captive colony at Peron Captive Breeding Centre (hereafter Peron CBC) was managed by the Western Australian Department of Biodiversity, Conservation and Attractions (DBCA), and was sourced from Bernier Island. In contrast, a captive breeding colony at Dryandra was sourced from Dorre Island. Two translocation attempts were unsuccessful-specifically the movement of 21 animals from Dorre Island to southern Dirk Hartog Island from 1974 to 1978, and a movement of 18 animals from the Peron CBC to Francois Peron National Park in 2001. The failure of these translocations was attributed to predation by feral cats, drought (Dirk Hartog Island translocation), and the impact of livestock [35][36][37]. The two successful translocations were each founded from the Peron CBC, with 91 individuals released onto Faure Island between 2004 and 2013, and 12 individuals into Wadderin Sanctuary in 2013. Faure Island is a 4561-hectare island in Shark Bay managed by the Australian Wildlife Conservancy (AWC) where all feral predators and livestock have been removed and was reported to support a population size of approximately 300 in 2017 [38]. Wadderin Sanctuary on the mainland (Figure 1d) is a 430-hectare predator-free enclosure managed by a community group [39] and now holds approximately 30 animals [40].
Currently, two new populations of L. fasciatus are being established. From 2017 to 2018, a population of L. fasciatus was reintroduced to a 7832 ha safe haven (a fenced area from which introduced predators have been removed) at Mount Gibson Wildlife Sanctuary. A total of 119 individuals were translocated to Mt. Gibson from Bernier (n = 40), Faure (n = 19) and Dorre (n = 60) Islands. In 2017, a translocation programme began for L. fasciatus to Dirk Hartog Island (58,640 ha), forming part of DBCA's "Return to 1616" project following the removal of exotic predators and herbivores [41]. Fifty-six individuals from both Bernier and Dorre Islands were proposed to be moved over 3 years (total n = 112). In this study, we aim to determine how the translocation history of L. fasciatus has affected the genetic health of all extant populations, particularly as serial translocations via intermediary captive populations have led to the possibility of genetic bottlenecks. We use genetic and population viability analysis models to explore the impact on source populations of concurrent harvesting scenarios and how the retention of genetic diversity can be maximised, and we predict growth rates and genetic diversity of the new Dirk Hartog Island and Mount Gibson populations. Most analyses are conducted under a baseline assumption of regular drought cycles that reduce reproductive success and survivorship, and we test the sensitivity of our results to both reduced and increased drought frequencies, as both scenarios are possible under climate change [42].

DNA Extraction and Nuclear Microsatellite Amplification
Genomic DNA was extracted from ear tissue using a standard "salting out" protocol [43]. Ten microsatellite loci were initially amplified using primers derived from the tammar wallaby (Macropus eugenii; Me14, Me17; [44]), yellow-footed rock-wallaby (Petrogale xanthopus; Y105, Y175, Y151, Y148; [45]), and allied rock-wallaby (P. assimilis; Pa593, Pa297, Pa385, Pa55; [46]). Briefly, for samples prior to 2016, polymerase chain reactions (PCRs) were carried out for individual microsatellite primers in a total volume of 30 µL with~100 ng DNA, 1X PCR buffer, 400 µM of dNTPs, 2 mM MgCl 2 , 0.2 µM of each primer, and 0.825 U Taq or, after 2016, in three PCR multiplexes using the Qiagen Multiplex PCR kit for a total volume of 7.5 µL per multiplex with~5-10 ng DNA, 1X Qiagen buffer, 0.2 µM primer mix, and water (Table S1). Fluorescently labelled DNA fragments were separated using an ABI373xl capillary sequencer (Applied Biosystems, Foster City, CA, USA) and scored manually with the aid of GENEMARKER software (v1.5, Soft Genetics, State College, PA, USA). Allele size was determined by co-running a Genescan500 standard (Applied Biosystems, Melbourne, Australia). Data were checked for input errors and duplicate genotypes using the Excel Microsatellite Toolkit add-in [46]. Deviations from Hardy-Weinberg equilibrium and linkage disequilibria were tested using GENEPOP v4.1.4 [47]. For markers with fewer than five alleles, a complete enumeration algorithm was used to estimate the exact p value; and for markers with five or more alleles, the Markov chain algorithm of Guo and Thompson [48] was used to generate an unbiased estimate of the exact p value. The presence of any null alleles was tested with MICROCHECKER v2.2.3 [49].

Statistical Analyses of Genetic Data
The programme STRUCTURE v2.3.4 [50] was used to estimate the number of genetically distinct populations and to assign individuals to populations. STRUCTURE uses a Bayesian clustering method to assign individuals to one of k populations and to estimate the degree of inter-population admixture. While some assumptions made by the software will likely not be met by the island populations in our system, such as the distribution of genotypes under Hardy Weinberg Equilibrium (HWE) and marker linkage equilibria, this approach still provides a useful assessment of population genetic divergence. Our models assumed no admixture between locations, but since in some cases locations represent recently sub-sampled source populations, and therefore could have shared allele frequencies, we ran models that assumed either correlated or uncorrelated allele frequencies between groups for comparison. After preliminary assessment of convergence times for the Monte Carlo Markov chain, a burn-in period of 100,000 steps was chosen, followed by 1,000,000 steps of the chain. To estimate k, four replicate runs at each value of k from 1 to 8 were performed, and the most likely value was estimated from the plot of ln Pr (X|k) vs. k, as well as from Evanno's method [51] which plots ∆ k (a second order rate of change of ln Pr (X|k)) vs. k, using CLUMPAK beta version [52]. STRUCTURE figures were generated using DISTRUCT v1.1 [53] in CLUMPAK [52]. To assess genetic similarity between populations using an approach that does not rely on the same biological assumptions, a principal coordinates analysis was conducted based on genetic distance, run in GENALEX v6.503 [54,55]. The level of genetic differentiation among all populations was determined by estimating pairwise F ST and Jost's Dest, also in GENALEX v6.503. Probabilities were calculated as the proportion of times the observed data generated values greater than values generated from 999 random permutations.
Standard genetic diversity metrics, including mean number of individuals per marker, mean number of alleles per marker, and observed and expected heterozygosities, were estimated for each sampling site, in GENALEX v6.503. Allelic richness, a measure of allelic diversity that controls for variable sample size, was estimated in HP-RARE [56] based on rarefaction to ten individuals. As marker Me17 was monomorphic except for one individual in the Faure population, in which it was homozygous for the alternate allele, allelic richness was calculated with and without this marker. Population inbreeding coefficients were measured using Wright's F IS in GENEPOP v4.1.4 [47]. Expected and observed heterozygosities, allelic richness, and FIS measurements were made for source populations separated into year cohorts. We estimated the genetic effective population size (Ne) using the linkage disequilibrium method on populations with a sample size of 25 or more [57][58][59], excluding singleton alleles (those that occur in one copy in one heterozygote) to prevent an upward bias in Ne estimation, as implemented in NeESTIMATOR v2.1 [60].
To assess whether any of the managed or natural populations underwent a genetic bottleneck, we compared the heterozygote distribution for each marker with the number of alleles for each population in BOTTLENECK v1.2 [61], using the two-phase substitution model with default proportions of the infinite alleles and single stepwise mutation models. A one-tailed Wilcoxon signed-rank test was used to estimate the likelihood that the observed data deviate from what is expected under mutation-drift equilibrium, the most appropriate test with fewer than 20 loci [62,63]. We also reported whether the Diversity 2020, 12, 448 6 of 22 allele frequency distribution deviated from an L-shaped mode, which can be a qualitative assessment of whether a population has passed through a genetic bottleneck [64].
Queller and Goodnight's [65] estimator was used to measure the mean pairwise relatedness within each population, as an indicator of genotype diversity and identity by descent within source populations. Whether the mean pairwise relatedness of a population was statistically different from 0 (defined as mean pairwise relatedness across the entire data set) was estimated by randomly permuting the data in the population pairwise matrix 999 times, and determining the proportion of permutations that gave a relatedness value greater than the observed value, as calculated in GENALEX v6.503.

Population Modelling
To estimate the number of founders needed to retain low-frequency alleles within newly established populations, scenarios were simulated in ALLELERETAIN v1.1 [66,67] implemented in R v3.6.2 [68]. ALLELERETAIN is an individual-based model that simulates population growth using a user-defined suite of parameters based on life-history traits. It estimates the probability of retention of selectively neutral alleles in founder populations over generations, with the starting frequency of these alleles set by the user. We were interested in assessing the impacts of translocating individuals of varying founder group size (N) on allele retention. To do this, we tested founder N values from 20 to 200 in increments of 20 individuals with allele frequency set to 0.05. We then repeated each scenario but considered the translocation to have occurred during a drought year. This was simulated by setting the initial survival after translocation to 0.4 and excluding reproduction for the first breeding cycle. For all simulations, we assumed that one L. fasciatus breeding cycle lasts 9 months [69]. Details and justification of the demographic parameters used in the model are provided in Table 1 and Table S2. Table 1. Life history parameters used for L. fasciatus population modelling in VORTEX and ALLELERETAIN. Parameter values were obtained from data in [33,34,41,70] and refined by C.S. and J.S. Where a parameter is not listed, defaults are used. V: VORTEX; A: ALLELERETAIN; EV: environmental variation; BHW: banded hare wallaby; DDR: density-dependent reproduction function.

Life History Parameter Value Used in
Species description Inbreeding depression V Lethal equivalents 3.14 V % due to recessive lethal 50 V EV concordance of reproduction and survival 0.5 V EV correlation among populations 0.5-0.8 V To estimate the impact of current translocation programmes on extinction probabilities and heterozygosities of the newly established L. fasciatus populations and the source populations, population viability analysis (PVA) simulations were run in VORTEX v10.17.2 [71,72]. VORTEX is an individual-based model that uses Monte Carlo simulations to estimate how factors intrinsic to individuals within populations alter growth rates, birth rates and extinction probabilities [72]. Initially, we developed a baseline PVA for a closed population with a starting N of 2000 individuals and a carrying capacity of 3000, where droughts at an average frequency of 1 in 6.25 calendar years reduced survival and reproduction by 50%, similar to empirical observations [33,34]. Although there is no empirical evidence for inbreeding depression, we included it in our models as there has been as yet no formal effort to detect it and populations do regularly pass through extreme population bottlenecks [33,34]. However, as populations appear to persist through these bottlenecks without obvious signs of depression, we halved the default number of lethal equivalents. Other island populations of marsupials are known to survive small effective population sizes and high levels of inbreeding [73]. Baseline model parameters are provided in Table 1 and Table S3.

Reproductive system
Various translocation scenarios (each of 1000 replicates) were then simulated that tested the number of founders and source populations needed to maximise genetic diversity and population growth for a translocation programme that, due to the logistical cost to the management agency, runs for a maximum of two years ( Table 2). Bernier Island was chosen for single source population translocations for two reasons. Firstly, as this population is smaller in size and has lower genetic diversity, it represents a conservative, worst case scenario, and secondly, this island is easier to access due to protected landing beaches. To quantify the impact of drought on translocation success, the best performing of these scenarios was re-run with both increased and decreased drought frequencies, and with no drought. Finally, a case study (scenario 8) was simulated based on the recent movements of L. fasciatus from Bernier and Dorre Islands to Dirk Hartog Island, and from Bernier, Dorre, and Faure Islands to the Mount Gibson Wildlife Sanctuary. Although the carrying capacity of Faure Island could theoretically be similar to Bernier and Dorre Islands (predicted here to be 3000) based on land area, there is currently no empirical evidence for this and the census population size remains well below 1000. We therefore tested two different carrying capacities for Faure Island, 3000 and a more conservative 1000. To further assess the impact of drought cycles in translocation planning, both the current and post-drought estimates of source population sizes were used as inputs in case study models. Table 2. Hypothetical and case study translocation scenarios * examined using population viability analyses.

Target Population Scenario Description
Dirk Hartog Island

Marker Performance
Of the ten microsatellite markers amplified, three were monomorphic in this species. The informativeness of the seven polymorphic markers varied, with the number of alleles ranging from two (Y175 and Me17) to eight (Pa593) ( Table S4). Expected heterozygosity (mean ± SE) ranged from 0.03 ± 0.03 (Me17), which was polymorphic in only one population, to 0.66 ± 0.02 (Pa593). The mean percentage of samples with missing genotype data across all populations ranged from 0.0% (markers Pa593, Me14 and Y105) to 5.1% (marker Y148). Marker Y175 violated the Hardy-Weinberg equilibrium (p < 0.05) in three populations out of six. Analysis in MICROCHECKER suggested the unexpectedly high homozygote frequency was caused by a null allele in all three populations. The percentage of samples with missing genotype data for this marker ranged from 0.0% (Dryandra) to 12.3% (Peron CBC), with a mean of 3.8%. After assessing the impact of null alleles on the estimation of global and pairwise FSTs in FreeNA [74], no significant differences were observed between estimates calculated with and without correcting for null alleles. Considering this, and the low number of markers, Y175 was retained for further analysis.

Genetic Diversity
Structure analysis revealed that the most likely number of genetically distinct clusters of L. fasciatus is two (k = 2; Figure 1e), representing the two remnant wild populations on Bernier and Dorre Islands. There was very good agreement in results whether allele frequencies were assumed to be correlated or uncorrelated between locations, and the absolute probability for k = 2 was slightly greater for correlated allele frequencies. Of the contemporary translocated populations, Faure Island and Wadderin are predominantly of Bernier Island origin, with negligible genomic contribution from Dorre Island. Of the extinct captive populations, Peron CBC was an extension of Bernier Island, whereas Dryandra was predominantly of Dorre Island origin. A small amount of potential admixture was detected between Bernier and Dorre Islands. In total there are ten individuals (four from Bernier and six from Dorre) that have an 80% or more probability of assignment to the other source population. Principal coordinate analysis agreed with the STRUCTURE results ( Figure 1f). Both pairwise FST and Jost's Dest measurements were similar (F ST = 0.026, p < 0.001; Dest = 0.029, p < 0.001), reflective of weak, significant divergence between the two parental populations. A summary of all pairwise FSTs is provided in Table S5, which shows that greatest distance exists between Bernier Island-derived populations and Dryandra (the sole Dorre Island-derived population). Pairwise relatedness analysis revealed that one source population (Bernier Island, p < 0.005), one extinct captive population (Peron CBC, p < 0.002), and one contemporary translocated population (Wadderin, p < 0.014) have relatedness values statistically greater than zero ( Figure S1).
Overall, genetic diversity was low across all populations (mean expected heterozygosity, HE = 0.34-0.45). Of the source populations, Dorre Island exhibited higher heterozygosity and lower inbreeding than Bernier Island, although allelic richness was comparable. Dorre Island had one private allele (alleles present in one population only) for marker Y105 (q = 0.02), Faure Island had two (Pa593, q = 0.05; Me17, q = 0.10), and Peron CBC had one (Y148, q = 0.01). When the two parental populations were compared only to each other, Dorre Island had four private alleles and Bernier Island had one. Diversity estimates within year cohorts for both Bernier and Dorre Islands showed a similar trend of higher than expected heterozygosities in the late 1990s, and lower than expected heterozygosities in more recent cohorts. Both expected heterozygosity and allelic richness have declined over the sampling period (Supplementary Figure S2). There were no significant differences in genetic variation between translocated and source populations except for Faure Island which showed a significant decrease in observed heterozygosity compared to Peron CBC using a one-tailed Wilcoxon signed-rank test (p < 0.05), and this island also indicated significant inbreeding due to non-random mating. Considering the two parental populations as one remnant L. fasciatus metapopulation, Faure has retained 93% of allelic diversity, while Wadderin has retained 87%. Calculating allelic richness without marker Me17 increases its value in all populations, but Faure no longer has the highest allelic diversity (Table S6). Caution is needed when interpreting all diversity results, however, as sample sizes for diversity measurements were often very low and varied between sites and temporal data.
Estimates of effective population sizes in extant populations reflected overall diversity measurements (Table 3), as Dorre Island had the highest Ne in the 2016-2017 cohort (Ne = 140, 95% CI 29-∞), which was 71% higher than for Bernier Island L. fasciatus sampled in the same period (Ne = 82, 95% CI 12-∞). Values of infinity for upper confidence limits arise from a lack of evidence for variation in the genetic characteristic which, in this case, is likely due to the small, homogeneous marker panel [59,60]. Our estimates of effective population size should therefore be considered cautiously and await confirmation with a more powerful dataset. Combining samples from all years, only one source population, Dorre Island, had statistical support (p = 0.039) for passing through a genetic bottleneck (Table 4). However, when Bernier and Dorre Islands were separated into year cohorts, both populations showed evidence of genetic bottlenecks in the earlier cohorts (1990s), and for Dorre Island also in 2013. Of the managed populations, Wadderin was the only population showing evidence of a genetic bottleneck (Table 4). Thirty individuals and 10 polymorphic loci are recommended to achieve power of at least 0.8 with the Wilcoxon signed-rank test however, so these results should be considered with caution. Both Wadderin and Peron CBC had an allele frequency distribution shifted towards more common alleles, consistent with a bottleneck, although this qualitative test requires at least 30 individuals and 10 to 20 polymorphic loci to be reliable [61].

Modelling Conservation Translocations
In the absence of periodic droughts, ALLELERETAIN demonstrated the strongest increases in the probability of retaining low-frequency alleles when the founder population size approached 20, and began to plateau after 40 ( Figure 2). Sixty founders are needed to retain 90% of alleles at a frequency of 0.05, and 80 founders are needed to retain 95% of alleles (recommended thresholds) at a frequency of 0.05. These numbers rose to 120 and 140 founders, respectively, when realistic impacts of drought were simulated. of current estimates, whereas growth on Faure Island was limited to 29% (Nc ~386). Interestingly, when Faure Island was assumed to have a K of 3000, there was little improvement in survival probability post-harvest if a drought census size was assumed (63% with K of 3000 vs. 60% with K of 1000). However, there was a dramatic increase in predicted Nc as 1400 was reached when harvesting occurred with a non-drought Nc and 900 with a drought impacted Nc, compared to 495 and 386 with a K of 1000, respectively.  Table 5. Viability of a translocated L. fasciatus population after 50 calendar years, testing four scenarios. One hundred individuals translocated in one or over two years, including a drought with a mean frequency of one every 6.25 calendar years. Carrying capacity was set at 10,000 and scenarios were run for 1000 replicates.  Population viability analysis modelling showed that when contrasting founder population sizes between 60 and 140, survival probabilities increase with increasing founder number until this relationship begins to asymptote from around 100 individuals (Table S7a-d). This level of harvesting also had limited impact on the source populations. Mixing the two source populations of Bernier and Dorre Islands increased heterozygosity in the translocated population, relative to only using founders from Bernier Island, and there was minimal impact of moving all individuals in one year versus splitting the translocation over two years (Table 5). Considering the additional logistical cost of a second year of translocation activities, the best scenario in terms of maximising survival probability, genetic diversity (heterozygosity and allelic diversity), and the final population size of the translocated population was to harvest the two source populations in one year (Scenario 4, Table 5). Table 5. Viability of a translocated L. fasciatus population after 50 calendar years, testing four scenarios. One hundred individuals translocated in one or over two years, including a drought with a mean frequency of one every 6.25 calendar years. Carrying capacity was set at 10,000 and scenarios were run for 1000 replicates. Simulations that varied drought frequency indicated the likelihood of an appreciable impact on the survival probability of a translocated population after 50 calendar years, dropping from 100% under the assumption of no drought, to 79% with a realistic drought frequency of 1 in 6.25 years, and 58% with a drought frequency of 1 in 5 years (Figure 3). There was also predicted to be a profound effect on the final N, as predicted carrying capacity (set to 10,000) was reached with no drought, whereas final N was estimated at 2277 individuals with a 1 in 6.25-year drought and 1064 at 1 in 5 years (Figure 3).

HE
0.369/0.374 0.344/0.343 P (surv): mean probability of survival; N: final mean population size of extant populations; HE = mean expected heterozygosity. * Numbers to the left of the slash are for a carrying capacity of 1000, and to the right are for a carrying capacity of 3000. Scenarios were run for 1000 replicates. Simulation of translocations recently conducted by management agencies between 2017 and 2018-harvesting Bernier and Dorre Islands for translocation to Dirk Hartog Island, and concurrent harvesting of Bernier, Dorre, and Faure Islands for translocation to Mount Gibson-including the impacts of 1 in 6.25-year droughts, predicted good survival probabilities over 50 years for the translocated populations of 83% and 84%, respectively. This result was also achieved even if harvesting occurred with a reduction in the size of source populations following drought (Table 6). In contrast, the timing of harvesting had a more significant effect on the source populations. Harvesting at current census population sizes (Nc) resulted in high survival probabilities (99%, 98%, and 93% for Bernier, Dorre, and Faure Islands, respectively; Table 6), assuming Faure Island to have a carrying capacity (K) of 1000. Bernier and Dorre Island population sizes stabilise at~82% of current estimates, whereas the Faure Island population stabilises after~65% growth from current estimates to a census size of around 500. If harvesting were to occur immediately after a drought, predicted survival probabilities reduced after 50 calendar years, most markedly to 60% on Faure Island, and there was a further reduction in predicted Nc in Bernier and Dorre Island populations to around 74% of current estimates, whereas growth on Faure Island was limited to 29% (Nc~386). Interestingly, when Faure Island was assumed to have a K of 3000, there was little improvement in survival probability post-harvest if a drought census size was assumed (63% with K of 3000 vs. 60% with K of 1000). However, there was a dramatic increase in predicted Nc as 1400 was reached when harvesting occurred with a non-drought Nc and 900 with a drought impacted Nc, compared to 495 and 386 with a K of 1000, respectively. Table 6. Comparison of population viability analysis (PVA) outputs after 50 calendar years when harvesting source populations using either current conservative source census sizes, or likely source census sizes following droughts. Models are based on scenario 8, describing recent and ongoing movement of L. fasciatus, and assuming an average drought frequency of 1 in 6.25 years.

Discussion
As the remnant range of L. fasciatus consists of just two adjacent offshore islands, the species-level impact of detrimental stochastic events affecting one or both islands could be profound. Establishing insurance populations and restoring the species to other suitable habitats free from introduced predators is a fundamental management objective, as is the retention of remnant genetic diversity to maintain adaptive capacity. While it is difficult to make conclusive summaries from the limited genetic data we provide here, we nonetheless attempt to explain observed trends that justify a more extensive genetic study. For example, our data suggest that both remnant island populations underwent genetic bottlenecks in the mid to late 1990s, potentially linked to the boom and bust cycles of L. fasciatus that occur in response to periodic droughts and thus a marked reduction in primary productivity [33,34]. Further, we show that the impact of past conservation actions (captive breeding and translocations) has possibly manifested in different ways in the two surviving translocated populations (Faure Island and Wadderin). Based on the limited samples available, Faure Island may be inbred due to non-random breeding but has relatively high allelic diversity, whereas Wadderin does not show signs of genetic inbreeding due to non-random breeding but may have passed through two selectively stronger bottlenecks (Bernier Island to Peron CBC, and Peron CBC to Wadderin). Our population viability and genetic modelling has revealed the importance of considering periodic fluctuations in population size when planning translocations and informed on the number of founders needed to avoid genetic bottlenecks. Taken together, these lines of evidence point to a critical need for genetic management and guidance for future translocations.

Genetic Diversity in Remnant Populations
Island populations are expected to lose genetic diversity as a result of genetic drift, especially in the absence of migration (i.e., gene flow) that acts to increase allelic diversity [75,76]. Here, we found that genetic diversity (allelic richness and heterozygosity) of L. fasciatus on Bernier and Dorre Islands was approximately two-fold lower relative to values reported for mainland populations of other Australian marsupials [77,78], and reviewed in [79], which is unsurprising considering the low number of polymorphic markers used in this study. Interestingly, values were similar to the Shark Bay island population of another hare-wallaby, the rufous hare-wallaby (Lagorchestes hirsutus), which was measured with shared markers and also had lower diversity than its mainland counterpart [80]. We also identified extremely low effective population sizes for L. fasciatus (Table 3) given census sizes of~2800 for Bernier Island and~2500 for Dorre Island [32,33], leading to good agreement in Ne/Nc ratios of 0.05 and 0.07, respectively. These estimates are several times smaller than those for other mammals, e.g., 0.2 to 0.8 reported in [81,82] and a mean of 0.75 reported in [83]. However, due to our low power to estimate Ne from genetic data, these trends require confirmation.
Bernier and Dorre Islands have been separated from mainland Australia for 8000 years (around 4000 generations), and from each other for around 5000 years [84,85], providing substantial opportunity for stochastic loss of allelic variation through genetic drift. The finding that around 10 remnant source population individuals (out of 130) were assigned with 80% or more probability to the other source population is therefore surprising. This observation could be consistent with recent gene flow, leading to migrant populations that have remained relatively isolated from the resident populations. However, this is unlikely, especially as it would have needed to have happened reciprocally. A more likely explanation is either a lack of discriminatory power in our markers, perhaps due to being developed in sister species [86], or there has been a technical issue, such as mis-labelling. Adequate sampling to properly determine genetic structure and the use of more powerful markers (e.g., single nucleotide polymorphic markers generated through next-generation sequencing), would resolve these alternative hypotheses.

Impact of Past Translocations on Genetic Diversity
The expected heterozygosities of either of the established translocated populations were not less than their source population (Peron CBC) or the parental population (Bernier Island). Despite describing two examples of serial translocations occurring from a single source population via a captive breeding programme, and hence providing opportunity for two bottleneck events, in only one case did we find evidence for a bottleneck-the Wadderin population translocated in 2013 with 12 founders. Wadderin was founded from the Peron CBC population, which itself was founded from Bernier Island (Figure 1d), and each translocation contributed a 2.4% and 5.7% loss of allelic diversity, culminating in a total loss from the Bernier source population of 8.1%. Faure, on the other hand, did not undergo a genetic bottleneck according to our analysis, and had a relatively high allelic diversity. The high allelic diversity is surprising considering its demographic history, and we found it to be driven by the discovery of an additional allele (marker Me17), present in one individual in homozygous form. While this allele could have arisen in the translocated Faure Island population, its independent verification would benefit from expanding sampling of Faure and Bernier Islands. While we were not able to detect a genetic signal of a bottleneck in the Peron CBC population, founded by 25 individuals, this does not preclude the possibility that a bottleneck has occurred that we were not able to detect or from which the population later recovered.
Since translocations began to Faure Island, however, observed heterozygosity appears to have decreased below equilibrium expectations and inbreeding coefficients have increased, possibly due to limited dispersal and increased breeding between related animals post-translocation. Another possibility is an increase in mating success heterogeneity [78]. These factors could also have contributed to the island's notably low effective population size estimate. This is of concern, as small populations are prone to higher rates of loss of genetic variation via genetic drift, especially if inbreeding depression leads to lower rates of reproduction and further reductions in population size [17]. However, the apparently high inbreeding does not appear to be associated with poor health of the Faure Island population, as it is well established with ongoing recruitment [87]. One explanation for this apparent paradox is that the boom and bust population cycles characteristic of Shark Bay L. fasciatus populations have led to a purging of lethal alleles, negating any detrimental effect of inbreeding [34]. A more prosaic explanation is that the Faure Island sampling cohort was both small (n = 10) and consists of a disproportionately high number of related individuals, leading to an overestimate of the inbreeding coefficient. A more comprehensive genetic analysis and a future study of the social dynamics of the Faure Island population could resolve this uncertainty.
While our bottleneck results require independent replication, a lack of bottleneck signal for the Peron CBC population with 25 founders supports results from our allele retention modelling. We show that a 70% probability of retaining a low frequency allele is maintained with a founder size of 25, but that retention probabilities drop dramatically with founder sizes less than 20. Other populations that have passed through bottlenecks tend to have greater founder numbers. For example, in birds, a bottleneck occurred when the Hawaiian gallinule (Gallinula galeata sandvicensis) was reduced to approximately 60 individuals [88], and reductions to 100 individuals may be sufficient to observe bottleneck effects in other avian species [89]. In the Maud Island frog (Leiopelma pakeka), simulations indicate that a sustained increase in irreversible allele loss begins when populations are reduced to 140 [90]. The apparent low power to detect genetic bottlenecks in translocated L. fasciatus populations could reflect relatively small changes in allelic diversity between source and founder populations (due to low effective population sizes and genetic variation in the source), compounded by our small panel of only moderately polymorphic markers. Further, the relatively rapid demographic recovery inherent to L. fasciatus during the boom phase of natural population cycles could mask any bottleneck signal that depends on heterozygosity excess [62]. Therefore, ongoing monitoring, including genetic monitoring, of translocated populations during their establishment would provide useful data with which to detect any detrimental genetic effects that might be initially masked.

Towards an Optimal Translocation Protocol for L. fasciatus
A clear trade-off exists in conservation translocation programmes between maximising viability of a new population and minimising the negative impact on critical source populations [90]. For example, sufficient founders are required to ensure a translocated population is buffered from some post-establishment mortality, as well as to retain rare alleles that bolster evolutionary potential [67,90]. Smaller founder numbers, on the other hand, reduce any negative impacts on source populations, but translocated populations will be more sensitive to mortality and the loss of rare alleles. Our PVA models for L. fasciatus show that for newly translocated populations, survival probabilities increase with increasing founder number until this relationship asymptotes at around 120 individuals, and the release of 100 individuals ensures good growth and survival in a hypothetical haven (island or fenced sanctuary with a carrying capacity of 10,000). That carrying capacity did not significantly impact 50-year predicted survival probabilities for the Faure Island population when using reduced drought-influenced census sizes, indicates the sensitivity of the successful establishment of L. fasciatus translocated populations to the founder number, and suggests that sufficient numbers are required to withstand regular population reduction caused by stochastic events. Importantly, harvesting the numbers recommended here has no apparent impact on the critical source populations, either in terms of population growth or survival probability over a 50-year period.
Our allele retention models predict that 80 founders are needed to retain 95% of low-frequency alleles, noting that retaining 90 or 95% genetic variation are recommended minimum thresholds for maintaining the evolutionary potential of populations [17,91]. Hence, aggregating the results of both modelling processes suggests that at least 100 individual founder animals are needed to maximise viability and retention of allelic diversity, a comparable figure to other taxa (e.g., 120 in frogs [90], 60 to 120 in passerines [19]). Further, by simulating the effects of periodic drought on population growth and survival, we have provided predictions about when best to translocate. For example, if harvesting were to occur immediately after a drought, we predict a 10-22% reduction in the census size of the source population after 50 years, relative to harvesting when populations are at their peak. In turn, we predicted that 140 founders would be required to retain 95% of alleles with a population frequency of 0.05 in a new population under drought conditions (as opposed to 80 founders under no drought); a 75% increase in the number harvested from the critical source populations. Therefore, climatic cycles appear to be an important consideration in translocation programmes, and we advise that the movement of animals should occur outside periods of low rainfall and when populations have recovered demographically from their impacts.
One way to increase genetic diversity of species that exist as only small and genetically depauperate remnants is to manage the species as a metapopulation [92]. However, whether to combine genetically distinct populations is contentious [1,7,78,[93][94][95]. While there is an immediate and obvious benefit of maximising genetic diversity and heterozygosity, thereby reducing the chance of inbreeding depression, mixing comes with risks, including outbreeding depression [1]. For example, if source populations derive from different habitats or climates, and have been long isolated, local adaptation can cause the rise of locally adapted alleles and for gene complexes to be in gametic disequilibrium. Locally adapted alleles can be diluted and gametic associations can be broken by recombination after mixing between differently adapted individuals, resulting in hybrid individuals with low fitness [96][97][98][99]. Other risks include competition for resources when there are morphological and size differences between source populations, and behavioural differences could affect mate choice, potentially reducing Ne and leading to a bias in representation of future generations (e.g., as occurs in the boodie Bettongia lesueur, [78]). Here, however, the habitats of the two adjacent L. fasciatus source populations have been similar across the 4000-5000 years (approx. 2000 generations) of their separation [79]. Further, there are no known differences in size, morphology, or mating behavior between islands, suggesting a low probability of outbreeding depression [95]. Considering the overall low genetic diversity maintained in this species, and by mixing we achieve a 10-12% increase in HE of the recipient population compared to using only Bernier Island as a source population, we advocate mixing the two source island populations in future translocations and predict a low risk of any adverse consequences.

Conclusions
Our work shows how an integrated analysis of genetics and population modelling can be used to inform management planning by simulating sensitive trade-offs involved in translocating threatened species. While much work has been done on the numbers of founders needed to establish new populations, less is known about the impact of harvesting for translocations on source populations [89], which is highly relevant for species with few remnant populations remaining. We suggest that the optimal translocation protocol for L. fasciatus is to mix the Bernier and Dorre Island populations in the same year, harvesting 60 individuals from each island source population (total n = 120). This gives a good probability of retaining low-frequency private alleles from each island, maximises survival probability and heterozygosity of translocated populations, limits logistical cost, and, crucially, has no major impact on the source populations. As the translocations currently underway to Dirk Hartog Island and Mt. Gibson approximate these recommendations, it will be valuable to monitor these populations over time to compare performance to model predictions. It is also worth considering Faure Island as an alternative source population to Bernier Island for future translocations, as was done with the Mt. Gibson translocation, to relieve pressure on one of the two remnant wild populations. However, in addition to possible inbreeding, our modelling shows Faure Island is more sensitive to stochastic events due to a smaller census population, further compounded if its carrying capacity is less than Bernier and Dorre Islands, and so the PVA models developed here should be used to assess the impact of any future harvesting. Wherever feasible, and after consideration of potential over-harvesting, future supplementations of the established translocated populations on Faure Island and Wadderin Sanctuary should derive from Dorre Island which, as a genetic mixing event, would be expected to increase gene diversity of the recipient population. Further, Dorre Island has a higher heterozygosity, more private alleles and lower inbreeding than Bernier Island. Finally, while we show that L. fasciatus seems resilient to harvesting and is well suited to translocation programmes, including multiple harvests from source populations, regular droughts, and limited carrying capacity, which have a substantial impact on population viability, particularly for populations with census sizes less than 500. We recommend that translocations are avoided during extended periods of drought and subsequent demographic recovery, due to a notably lower growth rate of the source populations after harvesting, and reduced capacity of new populations to retain allelic diversity.
Supplementary Materials: The following are available online at http://www.mdpi.com/1424-2818/12/12/448/s1, Spreadsheet S1: Genetic samples; Table S1: PCR multiplexes; Table S2: Details and justification of the demographic parameters used in ALLELERETAIN; Table S3: Details and justification of the demographic parameters used in VORTEX; Table S4: Performance of seven polymorphic microsatellite markers; Table S5: Population pairwise FSTs and probability values; Table S6: Comparing allelic richness with and without marker Me17; Table S7a-d: Summary of PVA outputs for translocation scenarios; Figure S1: Pairwise relatedness; Figure S2: Temporal change in genetic diversity for the two remnant source populations.