Admixture and reproductive skew shape the conservation value of ex situ populations of the Critically Endangered eastern black rhino

Small populations of endangered species risk losing already eroded genetic diversity, important for adaptive potential, through the effects of genetic drift. The magnitude of drift can be mitigated by maximising the effective population size, as is the goal of genetic management strategies. Different mating systems, specifically those leading to reproductive skew, exacerbate genetic drift by distorting contributions. In the absence of an active management strategy, reproductive skew will have long-term effects on the genetic composition of a population, particularly where admixture is present. Here we examine the contrasting effects of conservation management strategies in two ex situ populations of the Critically Endangered eastern black rhino (Diceros bicornis michaeli) , one managed as a semi-wild population in South Africa (SAx), and one managed under a mean-kinship breeding strategy in European zoos. We use molecular data to reconstruct pedigrees for both populations and validate the method using the zoo studbook. Using the reconstructed pedigree and studbook we show there is male sex-specific skew in both populations. However, the zoo’s mean-kinship breeding strategy effectively reduces reproductive skew in comparison to a semi-wild population with little genetic management. We also show that strong male reproductive skew in SAx has resulted in extensive admixture, which may require a re-evaluation of the population’s original intended role in the black rhino meta-population. With a high potential for admixture in many ex situ populations of endangered species, molecular and pedigree data remain vital tools for populations needing to balance drift and selection


Introduction
In the face of genetic erosion through defaunation (Bogoni et al. 2020), and range contraction (Britnell et al. 2023;Pacifici et al. 2020), ex situ populations of endangered species may serve as valuable reservoirs of genetic diversity and individuals for reintroduction (Farré et al. 2022).These, typically small, populations are vulnerable to inherent risks, including low resilience to perturbation and stochastic effects, and risk being driven into an extinction vortex through genetic drift and inbreeding depression (Williams et al. 2021).Genetic management of ex situ populations to minimize drift and preserve valuable diversity is therefore a major challenge for endangered species conservation and the restoration of resilient populations (DeWoody et al. 2021;García-Dorado & Caballero 2021;Hohenlohe et al. 2021).
One approach to minimize the effects of drift in ex situ populations is to employ active management strategies to reduce reproductive skew (Gooley et al. 2018).These may include male rotation strategies for group-living animals (Leus et al. 2011;Mucha & Komen 2016) or direct interventions to bypass competition, such as artificial insemination (Howard et al. 2016).In captive breeding programmes, current best practice to minimize drift is to equalise individual genetic contributions using a mean-kinship breeding strategy informed by a multi-generational pedigree.Under this strategy individuals are paired to minimize the average relatedness within the population (Ivy & Lacy 2012;Willoughby et al. 2015).An alternative management approach is to maintain animals in semi-wild conditions with natural mate choice.Doing so has the advantage of maintaining selection pressure to retain adaptive genetic variation, but at the potential expense of loss of genetic diversity due to high variance in mating success (Cain et al. 2014).Another factor that can impact the value of ex situ populations for reintroduction is if they have been founded or supplemented using animals of diverse or unknown origins, whether opportunistically or intentionally (Modesto et al. 2018;Senn et al. 2019).In this study we evaluate the conservation value and reintroduction potential of two ex situ populations of the Critically Endangered black rhino (Diceros bicornis, Emslie 2020) that exemplify the management strategies and challenges outlined above.
The black rhino has undergone a precipitous bottleneck from over 850,000 animals to only ~ 6000 today due to the poaching crisis (Moodley et al. 2017), and the ongoing threat of poaching necessitates the protection of a network of small and isolated populations (Emslie 2020).Ex situ populations are a key component of the meta-population, serving as sources of animals for translocations.But while small populations are easier to protect and manage, they suffer the compounding costs of genetic drift and inbreeding which are amplified with decreasing population size.Due to the severe bottleneck, the black rhino has suffered substantial genetic erosion (Moodley et al. 2017), so minimising further genetic drift in the extant small populations is important to maintain what little diversity remains.Reproductive skew has been identified as a priority concern for the species (Kenya Wildlife Service 2022).Female reproductive skew has been demonstrated among black rhino in both captive (Edwards et al. 2015a, b) and free-ranging populations, with variation in lifetime reproductive success (Harvey Sky et al. 2022), inter-calving interval (Patton et al. 2008) and age of first reproduction (Law et al. 2013).Male reproductive skew has been found to be even more extreme among freely mixing black rhino in semi-wild conditions (Garnier et al. 2001) with dominant males securing the highest number of matings with females (Cain et al. 2014).
Previous management of the species as four 'ecotypes', that became associated with historic subspecies classifications (eastern D. b. michaeli, southern D. b. minor, south-western D. b. bicornis and western D. b. longipes, Rookmaaker 1995), has left a high potential for admixture in ex situ populations.Based on recent genomic analyses, the eastern black rhino is considered to comprise multiple evolutionary significant units (ESUs, Sánchez-Barreiro et al. 2023), which have likely been mixed both in situ (Muya et al. 2011) and ex situ.The southern ESU (now reclassified as D. b. bicornis, Sánchez-Barreiro et al. 2023) is considered less diverse, but shows significant divergence from the eastern clade, estimated to have diverged 641Ka (Moodley et al. 2020).These two deepest clades are divided by the Zambezi river which serves as a significant barrier to animal movement, although rare instances of admixture along this boundary are still possible (Sánchez-Barreiro et al. 2023).Admixture has the potential to release the burden of inbreeding depression, but with increasing divergence may result in outbreeding depression (Frankham et al. 2011).Where admixture and reproductive skew occur simultaneously, the introgressed alleles may be eliminated, or conversely, may spread disproportionately throughout the population.In its most extreme form this results in a genetic sweep and the erosion of 'native' alleles (Adams et al. 2011;Gibson et al. 2019).Therefore, reproductive skew not only accelerates genetic drift but can dramatically alter the evolutionary trajectory of a population in the context of admixture.
Pedigrees form a critical management tool for the black rhino artificial meta-population, as it is almost entirely reliant on human-mediated translocations to maintain gene flow (Emslie 2020).Accurately constructing and maintaining a multi-generational pedigree requires intense sampling and long-term monitoring of a population (Galla et al. 2022), which is challenging for free-ranging animals.Furthermore, the accuracy of primary assumptions dictates the subsequent accuracy of a pedigree.For example, founders are often assumed to be unrelated which, if untrue, may result in highly inbred populations with negative long-term implications for breeding programmes (Hogg et al. 2019).Conversely, if founders are incorrectly assumed to belong to the same evolutionary group, this can result in genetic admixture in populations of conservation concern which may compromise their conservation value (Hvilsom et al. 2013).By combining observational data with molecular information, many of the assumptions and pitfalls of an incomplete pedigree can be addressed.Molecular data has been used to correct studbook records for captive African painted dogs (Lycaon pictus, Miller-Butterworth et al. 2021), captive and reintroduced Californian condors (Gymnogyps californianus, Moran et al. 2021), and has revealed relatedness amongst founders in populations of the Critically Endangered Kākāpo (Strigops habroptilus, Bergner et al. 2014).
In this study we incorporate molecular data with studbook records for the management of two key ex situ populations of the Critically Endangered eastern black rhino: the population of eastern black rhino in European Zoos managed under the European Association of Zoos and Aquaria (EAZA) Endangered Species Programme (EEP), and an ex situ population of eastern black rhino that is extralimital in South Africa (SAx).Admixture with the southern black rhino is known to be present in the SAx population, although the full extent of this has not been established (Kretzschmar, pers. comm.).Understanding admixture in conservation populations is key to inform both conservation and management strategy, and will dictate the suitability of individual animals for different reintroduction plans.
Here we use the genotypes obtained from a panel of 15 autosomal microsatellite loci combined with mitochondrial d-loop haplotypes to: 1) validate the marker panel for pedigree reconstruction using the studbook of European zoos and reconstruct the pedigree of the SAx population of admixed black rhino; 2) compare the genetic diversity maintained in the two ex situ populations under different management strategies (mean-kinship breeding in the European zoos and free-ranging in SAx); 3) evaluate the effectiveness of the management strategies at moderating reproductive skew and therefore genetic drift; and, 4) quantify the contribution of southern ancestry to the zoo and SAx populations.Based on these analyses we discuss the implications for conservation strategies utilizing animals from these two major ex situ populations of eastern black rhino.

Methods
Ethical approval for sampling of European zoo animals was granted by the EEP Taxon Advisory Group for rhino.Ethical approval for the use of veterinary by-product samples from South Africa was granted by the University of Manchester, UK.Tissue samples from SAx were imported from South Africa under CITES export permit number 279652 and CITES import permit number 563301.Details pertaining to the locations of animals used in this study are not disclosed due to the confidentiality of black rhino data which is a policy of the IUCN SSC African Rhino Specialist Group.

Study populations
The EAZA Endangered Species Programme (EEP) for the eastern black rhino is maintained in European zoos and managed using a detailed studbook.The studbook contains comprehensive records dating back to the early 1960s, while some lineages can be traced back to 1948.The population currently numbers around ~ 90 individuals (Pilgrim, pers.comm.) and breeding recommendations are made using a mean-kinship strategy (Lacy et al. 2011).
The second ex situ population (SAx) is extralimital and now held on a private reserve in South Africa.The founding animals were captured in Kiboko, Kenya, in 1961, and translocated to Addo National Park.This population was intended as an insurance population at a time when Kenya's black rhino were threatened by steep decline (Emslie & Brooks 1999;Muya et al. 2011;Okita-Ouma et al. 2007).However, three males of a divergent ESU, the southern black rhino (D. b. bicornis), were transferred into the ex situ population while in Addo National Park in 1977(Hall-Martin 1984).As a result, the population contains an unknown proportion of southern admixture (Kretzschmar, pers. comm.).The original intention of this insurance population was to return SAx descendants to Tsavo National Park to supplement the small in situ population which had declined to only 15 animals in 1991 (Hall-Martin 1991).However, admixed animals may not be favourable for translocation into native populations in East Africa.

Sampling and genotyping
Samples were collected from zoo animals, between 2017 and 2021, using non-invasive nasal swabs (N = 75), or whole blood (N = 10) where this could be obtained opportunistically.One muscle tissue sample from a deceased animal, preserved in 2ml absolute ethanol at -20°C, was also included.Nasal swabs were preserved using silica gel at room temperature and blood samples collected in EDTA tubes were frozen at -20°C.From the SAx population we extracted DNA from pinna offcuts preserved in sodium chloride.These were obtained during routine ear-notching operations performed to facilitate individual animal identification on the reserve.Samples used in this study dated up to 2018, but samples from historic assignments and some adults were not available.We sampled approximately 80% of the SAx and 85% of the EEP population.
DNA was extracted using a QIAamp DNA mini kit (QIA-GEN GMBH, Hilden, Germany).Mitochondrial DNA was amplified using D-loop primers mt15996L and mt16502H following Moodley et al., (2017).Individuals were Sanger sequenced in the forward and reverse direction (Eurofins Genomics Europe) and consensus sequences generated using Geneious software (Geneious Prime 2023.0.1, https:// www.genei ous.com).These sequences were aligned with the corresponding homologous sequenced region from the Moodley et al., (2017) metapopulation database and a haplotype network was constructed using PopArt (Version 1.7, Leigh & Bryant 2015).
Using the 15 markers we calculated mean allelic richness (AR), and observed and expected heterozygosity (H O and H E , respectively) for each population using R packages PopGenReport and Poppr (Adamack and Gruber 2014;Kamvar et al. 2014;Adamack & Gruber 2014;Kamvar et al. 2014) Population level differentiation, F ST , and the inbreeding coefficient, F IS , were calculated using a Weir and Cockerham estimator in R package hierfstat (Weir & Goudet 2017).We estimated the effective population sizes for the EEP and SAx populations using the linkage disequilibrium method, under random mating and with a critical minor allele frequency of 0.05, implemented in N e Estimator v2 (Do et al. 2014).

Pedigree reconstruction
Pedigree reconstruction for the EEP and SAx was carried out on the microsatellite dataset using FRANz (version 2.0.0,Riester et al. 2009).Pedigree reconstruction was implemented for each population separately by adding observed mothers in a stepwise fashion to the FRANz input file.For the EEP, we ran FRANz without the full-sib heuristic, with a female reproductive span of 5 to 30 years of age, a male reproductive span of 5 to 40, and ran the simulation for 1,000,000 iterations.We first ran FRANz without any known relationships included and retained any consensus assignments to known mothers in the second model, to which we then also added mother-offspring pairs with the rarest mtDNA genotypes.After each run, we revised the input file of known relationships to remove any parents with a posterior probability of < 0.8 and added further known mothers based on the frequency of the mtDNA genotype in the population.We averaged the posterior probability for each parent over 20 repeat runs once all mothers were added and removed any assignment with a probability < 0.8.We compared the resulting pedigree to studbook records and counted the number of erroneous or missing assignments to benchmark this method.Marker DB49 showed the highest number of mismatches with the studbook, but removing this marker did not improve the accuracy of pedigree reconstruction in comparison with the studbook data.
For the SAx population we repeated the above mentioned process.A mother-offspring studbook has been curated from observational data collected during ear-notching operations.We first used mitochondrial (mtDNA) genotypes to remove any erroneous mother-offspring relationships.Where a different mother to the observed mother was assigned, the microsatellite genotypes of the candidate parents were compared with the offspring and any assigned father to check for correct assignment.Finally, we added additional parents for unsampled individuals to increase network connectivity.These were based on known relationships from observational data (mother-calf relationships, (Hall-Martin 1986), or previous genetic relationships (Harper pers. comm., unpublished data, SI Table 2).We ran the final configuration 20 times to generate consensus probabilities for all parent offspring pairs and eliminate spurious assignments.The resulting pedigree network was visualised using Helium (Shaw et al. 2014).

Quantifying admixture
In order to quantify the proportion of southern ancestry across the pedigrees we joined the reconstructed SAx pedigree produced by FRANz, with historic relationships assigned from a combination of historic records (Hall-Martin 1986) and molecular data (Harper pers. comm.).Based on the early assignments among population founders we used the R package OptiSelect (Wellmann et al. 2012) to quantify the contribution of each male founder to both the zoo and SAx pedigrees.

Evaluating reproductive skew
Here we consider that reproductive skew is a result of the intrinsic (individual fitness) and the extrinsic (management strategy) components.To evaluate variation in reproductive skew and genetic drift between the two populations, we used the breeding histories from the studbook and the reconstructed pedigree.For the EEP we restricted our analyses to the 32-year time period between 1990 and 2022 to reflect the current management strategy.We excluded any offspring that were stillborn or did not reach one year of age.We compiled a list of candidate parents including all animals that were alive in the EEP between 1990 and 2022 and were > 5 years of age (Brodie et al. 2011).Animals were excluded as candidate parents on death or on transfer out of the zoo breeding population.For each candidate parent we recorded the years they were assigned as sire or dam in the list of offspring born between 1990 and 2022.
For the SAx population we restricted our analyses to the period between 2003 and 2018 as the last males were translocated into the population in 2001 and black rhino have a gestation period of a minimum of 15 to 18 months (Le Roex & Ferreira 2020;Okita-Ouma et al. 2021).To compare measures of skew to the zoo population we used a minimum age of five years for candidate parents, even though freeranging males tend not to reproduce successfully until they are competitive at nine or ten years of age (Law & Linklater 2014).In the SAx population some males are located on an adjoining reserve to reduce conflict.Males were occasionally transferred between these two populations or at times held in an enclosure with females.We included all animals aged > 5 as candidate parents to reflect the management strategy.We used the reconstructed pedigree to list the births assigned to each candidate parent.For each population, we calculated the multinomial index of reproductive skew, using the R package SkewCalc (Ross et al. 2020), which accounts for a non-linear relationship between age and reproductive success.If < 5% of the posterior distribution of M falls within the Region of Practical Equivalence (ROPE) of − 0.1 to 0.1, the null hypothesis of no skew can be rejected.
To evaluate how well each population has preserved founder genetic diversity we quantified the contribution of each founder to the current population (2022) for males and females from zoos, and the contributions of male founders to the current and total SAx population.We were unable to make accurate estimations for female founders from the SAx population as early relationships in the pedigree were unassigned for animals that were not sampled in this study.

Mitochondrial diversity
A total of 177 (SAx = 103, EEP = 74) mitochondrial d-loop consensus sequences of length 456 bp were obtained, featuring 28 polymorphic sites (Genbank OR709675-OR709685).We identified 10 haplotypes, of which 3 were shared between the EEP and SAx, one was unique to SAx, and 6 haplotypes were unique to the EEP.Two of the EEP haplotypes fall within the southern ESU (Fig. 1, SN; Moodley et al. 2017), the first of which is from a male 'Kalusho' from the discontinued D. b. minor EEP, and the second from a D. b. michaeli EEP male 'Sammy' indicating southern maternal ancestry.
All of the four SAx haplotypes had an eastern origin, confirmed by comparisons with the Moodley et al., (2017) historic mtDNA dataset (Fig. 1).This low number reflects the small number of female founders.Based on their relationships in the pedigree we were able to infer that the original female founders 'Ida' and 'Brunni' share an eastern mtDNA haplotype.This is likely due to the close proximity of their capture sites in Kiboko, Kenya and indicates a common ancestor.A second mitochondrial haplotype was associated with the lineage of female 'Kate'.Of the final two haplotypes found in the SAx population, one was found only in a male founder from Tanzania, 'Richard' and one originated from the EEP population through female 'Tana', one of two captive-born females translocated from Port Lympne Safari Park in the UK to the SAx population in 2004 (King & Beer 2018).
Both populations had similar expected and observed heterozygosity (EEP: H E = 0.665, H O = 0.706, SAx: H E = 0.673, H O = 0.699) (Table 1).Private alleles were higher in the EEP (28) than in SAx (6) consistent with the larger number of founders.Population F IS was not significantly different from zero in both ex situ populations (EEP -0.062, SAx-0.039),indicating no excess of inbreeding.N e was estimated at 23.3 and 12.8 for the EEP and SAx populations respectively using the linkage disequilibrium method (N e estimator v2, Do et al. 2014).The two populations showed weak differentiation with an F ST value of 0.029.

Validating pedigree reconstruction using the EEP studbook
Based on 48 known maternal relationships, 30 known paternal relationships were correctly generated by FRANz with reference to the studbook data.78 relationships were correctly assigned as missing, as the parent was not sampled in our dataset.Eight paternal relationships were unassigned despite the father being present in the dataset with a posterior probability of < 0.8.Four parent-offspring relationships were incorrectly assigned by FRANz with two being assigned to full-siblings and two to half-siblings.

Reconstructing the pedigree of the SAx population
The proposed pedigree includes a total of 168 relationships within the population.77 mother-offspring relationships were assigned based on our molecular data, of which 69 confirmed observational assignments.FRANz highlighted eight false assignments in the observational pedigree based on > 2 mismatched markers between mother and offspring genotypes.In two of these cases the correct mother was able to be assigned with molecular data, whereas for the other mismatches the mothers were presumably not sampled.We were able to assign paternity for 68 animals based on microsatellite data.Five impossible assignments resulting from an unsampled father were removed due to either spatial considerations (animals held on different parts of the reserve), or the very young age of the candidate father (6 years old).Additional relationships in the pedigree were Broadly, the groups Eastern (EA) and Central (CE) represent an eastern cluster (blue shading), and groups South-Eastern (SE), South-Northern (SN) and South-Western (SW) represent a southern cluster (gold shading).All maternal lineages bar two EEP animals have an eastern origin.Of these two southern haplotypes in the zoo population, one is an individual from the discontinued D. b minor breeding programme, and one is an admixed male with maternal southern ancestry added from observed assignments or previous assignments from molecular data which were not genotyped in this study (Kretzschmar, pers comm.).None of these additional assignments caused the pedigree to 'break', but these relationships should be considered candidates rather than confirmed (SI Table 1).For a visualization of the pedigree, see Fig. 2.
Free-ranging female black rhino may produce their first offspring aged five or six (Okita-Ouma et al. 2021), the earliest age of reproduction documented in the SAx population was six years old (N = 3).Most females will have produced an offspring by age nine (Harvey Sky et al. 2022).Of the 25 sampled candidate mothers aged over six at the birth of the youngest animal, 18 accounted for the 74 births for which maternity could be assigned.Four animals in the candidate parent group had not reproduced by the age of 9.The average inter-calving interval (ICI) for females ranged from 2 to 3.33 years, and four females had an average intercalving interval < 2.5 years indicating a highly productive population.
Male black rhino do not become reproductively competitive until the age of nine or ten (Adcock 1994).For genotyped males aged > 10 at the approximate conception time of the most recent cohort, this span included 67 births spread over 22 of the 26 candidate fathers.Unfortunately, samples were not available for all the population, therefore genotypes for several animals are missing.One ungenotyped male 'Bwana' is suggested to have contributed significantly to the population based on observations.As a result, one of his offspring was assigned a high proportion of the paternal relationships by FRANz.This was subsequently corrected based on the offspring's young age (< 6 years old) and Bwana was listed as the likely parent.

Founder contributions
Under an ideal scenario, where drift is inconsequential, each founder lineage would have equal representation in the population.For the zoo population, founder contributions among females varied from 17.3% (Jimmi, studbook number 175, Tsavo National Park, Kenya) to 0.6% (June, studbook number 17, unknown origin; Jarca, studbook number 178, Tsavo, Tanzania, Fig. 4A).For male founders of the European zoo population, contributions also varied, ranging from 13.1% (Gareth Edwards, studbook number 534, from SAx (Addo) South Africa) to 0.7% (Conni, studbook number 36, recorded as from 'Kenya Game Ward', Fig. 4A).Although Gareth Edwards originated in the SAx population, based on his assigned relationships he is predicted to be unadmixed.Conversely, male Madiba who also originates in SAx is estimated to have 25% southern ancestry by the paternal relationship assigned from the pedigree reconstruction.Likewise, an examination of records compiled by Rookmaaker et al. (1998), revealed that male 'Sammy', imported from Japan has 50% southern ancestry with his maternal grandparents both imported from Namibia (Rookmaaker et al. 1998).Consequently, while the genetic makeup of the European zoo population is predominantly eastern black rhino, admixture is also present (3.5% of the current population).
In the SAx population, founder contributions among males show strong asymmetry (Fig. 4B) with one of the original founders 'Darkie' contributing only 9.6% of total male ancestry (Table 2).The other original male founder from Kenya, 'J.A.' has contributed a higher proportion of ancestry (22.76%) while more recent founders 'Bwana' translocated from Port Lympne Safari Park in the UK in 1995, and 'Richard' from Tanzania contributed intermediate proportions of ancestry (12.57% and 15.72%, respectively).However, the highest proportion of ancestry is assigned to the southern black rhino male (SBR, 39.35%, Table 2).
Actual values may deviate slightly due to recombination and the random segregation of gametes.The southern black rhino male has contributed the greatest proportion of male ancestry to the SAx population.Subspecies assignments are inferred based on capture location.

Discussion
Following a decline of over 96% in the wild (Harley et al. 2005), the two ex situ populations examined here represent an important and substantial proportion of the global black rhino meta-population.In this study we have reconstructed the pedigree of the freely mixing SAx population with a high degree of certainty that enabled a reliable assessment of reproductive skew and the contribution of southern ancestry.This has enabled a real-world comparison of two contrasting management strategies commonly employed for endangered species.Elucidating the genetic effects of population management in black rhino, in contrast to animals with a shorter generation time (e.g. the Tasmanian devil, Sarcophilus harrisii, Farquharson et al. 2021;Wright et al. 2020), will have A detailed understanding of diversity, vulnerability to drift, and individual proportions of admixture in these two key populations offers exceptional scope for detailed management strategies that will maximise their conservation value and begin to bridge the gap between ex situ and in situ conservation.
The SAx population has reasonable levels of diversity in spite of its small number of founders.However, the European zoo population retains a higher diversity of microsatellite alleles and mtDNA haplotypes.Indeed, the European zoos contain high genetic diversity based on microsatellite genotypes, comparable to the diversity observed in Kenya (H O = 0.73 ± 0.03, Muya et al. 2011), and substantially higher than estimates for the southern black rhino (H O = 0.436, Harley et al. 2005) and the south-western black rhino (H O = 0.523, Harley et al. 2005).
Maintaining this diversity over time requires the effects of genetic drift to be minimised so that rare variants are more likely to persist over generations.In this study we show that the mean-kinship breeding strategy in the zoo effectively reduces reproductive skew in comparison to a freely mixing population with low intervention.Both the reduction in reproductive skew and the larger population of founders are reflected in the larger effective population size of the European zoos.Conversely, the SAx population has a small effective population size which, with sustained high reproductive skew, means this population will be subjected to stronger effects of genetic drift and will likely lose diversity more quickly.Increasing the effective population size must be a priority given the precarious future predicted by our findings if this population were to remain genetically isolated.
Acting to counter drift by maximizing the effective population size is most important when the population is small.Selection becomes more effective and drift less extreme as the population becomes larger.In comparison to the effective population size of large and connected historic populations, both of these ex situ populations would be considered small, but there is also an argument to be made to maintain the semi-natural breeding system and associated level of reproductive skew in the context of the SAx population.Implementing a mean-kinship breeding strategy circumvents reproductive selection, which is maintained in the SAx population, as to an extent are many 'natural' ecological processes.Altered selection can result in inefficient purifying selection against mutations that may be deleterious in situ, and an increase in the frequency of these mutations may compromise fitness (Wilder et al. 2020) and the success of reintroductions.
Adaptation to captivity, is a major concern for ex situ populations (Frankham 2008), and a particular concern for populations in zoos.Equalizing contributions to the subsequent generation by mean-kinship breeding, rather than selecting for fecundity, should function to slow the rate of adaptation to captivity (Williams & Hoffman 2009), but purifying selection may still be relaxed.Although the mean-kinship breeding programme is effective at reducing the current degree of reproductive skew in the zoo population, founder contributions in this population are highly unequal.Moreover, Edwards et al. (2015a, b) recorded a potential founder population size of 135 wildcaught individuals, of which only 41 have contributed to the extant zoo population, representing a loss of 69% of potential lineages.Nonetheless, we found no evidence of reproductive skew for females in subsequent generations.
Fig. 3 Multinomial index of reproductive skew (M) for males and females, accounting for age.No skew was present for females in either populations as distributions were not different from zero.Skew is found among males in both populations, but is significantly higher for the free-ranging SAx population than for the European zoo population which is managed by mean-kinship breeding (median of difference 0.72, 95% CI 0.12 to 1.72) It is possible that this later lack of skew reflects improved husbandry, or alternatively early variation in reproductive success may have selected for lineages that do well in captivity.Interestingly, skew was still present among males in zoos.This may reflect the fact that males are often favoured over females for translocations intended to inject new diversity into a breeding programme.Some males will therefore reflect more unique lineages and will be favoured by the mean-kinship breeding strategy due to their low relatedness to the population.Unfortunately, two of these imported males 'Sammy' and 'Madiba' are admixed and this has resulted in the spread of a low proportion of southern ancestry in the European zoo population, although currently this is confined to certain lineages.
Indeed, the diverse founder composition of the zoo population, primarily originating in East Africa, has likely resulted in the interbreeding of several eastern evolutionary units (Sánchez-Barreiro et al. 2023).As a rule, the translocation of admixed animals from ex situ populations into existing native populations in situ is a high-risk strategy and is actively discouraged (Bertola et al. 2022).However, this may not preclude translocations of zoo animals from mixed eastern ESUs into the native range, given that a similar Fig. 4 A Contributions (%) of each founder to the current European zoo population based on the studbook.For founders with lower contributions, these lineages are likely to be lost from the population.B contributions (%) of each male founder to the current SAx population.The southern male founder has contributed the greatest population of ancestry to the SAx population mixing of animals has likely also taken place in Kenya (Muya et al. 2011).The zoo population may therefore serve as a useful reservoir of diversity for future translocations, even where native eastern populations persist, but it is key that any translocations be genomically informed to prevent any further erosion of rare native diversity.
The conservation value of eastern black rhino with southern ancestry remains an open question.Admixture between subspecies can result in increased fitness of hybrid offspring through heterosis due to the masking of deleterious alleles (Onley et al. 2022;Wei & Zhang 2018).However, outbreeding depression may result from the breakdown of co-adapted gene-complexes (Edmands & Timmerman 2003;Frankham et al. 2011).This gamble is not merited by the current status of the black rhino, particularly considering the deep 641Ka divergence between eastern and southern clades (Moodley et al. 2020).This divergence far exceeds the recommendations regarding admixture which are intended to minimize the risk of outbreeding depression (Frankham et al. 2011).
In the zoo population, southern admixture is recent and limited, therefore pairing decisions could be made to prevent further spread of southern alleles.However, the situation in the SAx population is more complex.Admixture from the subspecies D. b. bicornis is extensive in this population, but not reflected among mitochondrial haplotypes due to its paternal origin.Indeed, the southern male, has contributed a greater proportion of ancestry than any other male founder.This southern male, added to the population on Addo, in 1977, to replace the dominant eastern bull after his death (Hall-Martin 1984), was the only breeding male over several years at a time when the population was small (ca.15 individuals).This proliferation of southern ancestry is an accident of timing and historic management rather than sexual selection but has led to the extensive introgression of southern ancestry throughout the SAx population.Management strategies including the translocation of eastern males into the population are currently being implemented to reduce the proportion of admixture and increase N e , but the translocation of admixed animals continues to pose a dilemma.
The original intent of the SAx population was to return animals to where the eastern founders were captured, but the growth of the population in the Tsavo region of Kenya (Chyulu hills, Tsavo East, Tsavo West, Ngulia and the Intensive Protection Zone (IPZ)) to around 363 black rhino (Kenya Wildlife Service 2022), means translocations of admixed SAx animals would be better prioritized elsewhere.Contractually the SAx animals may not be translocated elsewhere in South Africa, being regarded as an alien subspecies (Emslie 2020).Translocations to Tanzania have already been undertaken (Kretzschmar, pers.comm.), however a better strategy may be to introduce these animals to regions where admixture may have historically occurred at the boundary of the D. b. michaeli and D. b. bicornis range in Malawi, Zambia and Southern DRC.Alternatively, these animals would be a potential source for the founding of new populations where black rhino are regionally extinct, as is the case for recent translocations to Rwanda and Chad (Bantlin, pers. comm.)This study contributes to the relatively small number of detailed studies on genetic management (e.g.Tasmanian devil (Sarcophilus harrisii), Farquharson et al. 2021;Wright et al. 2020; Kākāpo (Strigops habroptilus), Guhlin et al. 2023).Our findings therefore have wider implications for how we manage meta-populations of other Critically Endangered species.Firstly, generating reliable data on genetic diversity and reproductive skew under different management strategies can inform decisions on the optimal means to ensure long-term maintenance of genetic diversity in small human-modified populations.Secondly, as for the black rhino, admixture is a frequent outcome of unknown founder diversity in many ex situ populations (Hvilsom et al. 2013;Węcek et al. 2016).This study demonstrates the importance of tracing and identifying ancestry to inform appropriate translocations, integrate ex situ populations into meta-population level strategy and realise their conservation potential.By incorporating genetic information into translocation strategy, risk can be minimized, translocation success improved and the long-term benefits from these costly exercises maximized.
In compliance with the Nagoya agreement, samples from South Africa were imported with prior informed consent and an access and benefit sharing agreement with permit number: AWB No. 7220347514.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Fig. 1
Fig.1Median-joining haplotype network of mitochondrial d-loop sequences.Haplotypes from the two ex situ populations are presented in the context of the data and evolutionary groups from theMoodley et al. (2017) meta-population study.Broadly, the groups Eastern (EA) and Central (CE) represent an eastern cluster (blue shading), and groups South-Eastern (SE), South-Northern (SN) and South-Western (SW) represent a southern cluster (gold shading).All maternal lineages bar two EEP animals have an eastern origin.Of these two southern haplotypes in the zoo population, one is an individual from the discontinued D. b minor breeding programme, and one is an admixed male with maternal southern ancestry

Fig. 2
Fig. 2Complex pedigree of the SAx population.Females are represented by light blue, and males by light green circles.Maternal relationships are indicated by red lines and paternal ones by blue lines.Relationships above the dotted line are derived from observational records from the SAx population while it was at Addo Elephant National Park (Hall-Martin 1986) and some assignments are based on molecular data generated by Harper (unpublished data).Individuals below the black line are assignments from this study.Where samples weren't available for genotyping, for some individuals, paternity was inferred by excluding sampled candidate fathers

Table 1
Metrics from mitochondrial and microsatellite genotypes for the zoo (N = 85) and SAx (N = 110) populations