Sex-biased inbreeding effects on reproductive success and home range size of the critically endangered black rhinoceros

Linking diet switching to reproductive performance across populations of two critically endangered mammalian herbivores

Optimal foraging theory predicts that animals maximise energy intake by consuming the most valuable foods available. When resources are limited, they may include lower-quality fallback foods in their diets. As seasonal herbivore diet switching is understudied, we evaluate its extent and effects across three Kenyan reserves each for Critically Endangered eastern black rhino (Diceros bicornis michaeli) and Grevy's zebra (Equus grevyi), and its associations with habitat quality, microbiome variation, and reproductive performance. Black rhino diet breadth increases with vegetation productivity (NDVI), whereas zebra diet breadth peaks at intermediate NDVI. Black rhino diets associated with higher vegetation productivity have less acacia (Fabaceae: Vachellia and Senegalia spp.) and more grass suggesting that acacia are fallback foods, upending conventional assumptions. Larger dietary shifts are associated with longer calving intervals. Grevy's zebra diets in high rainfall areas are consistently grass-dominated, whereas in arid areas they primarily consume legumes during low vegetation productivity periods. Whilst microbiome composition between individuals is affected by the environment, and diet composition in black rhino, seasonal dietary shifts do not drive commensurate microbiome shifts. Documenting diet shifts across ecological gradients can increase the effectiveness of conservation by informing habitat suitability models and improving understanding of responses to resource limitation.

Effects of dehorning on population productivity in four Namibia sub-populations of black rhinoceros (Diceros bicornis bicornis)

The black rhinoceros (rhino) (Diceros bicornis) is critically endangered, with poaching being one of several threats to the species’ survival. Many reserves across several countries, including Namibia, South Africa, and Zimbabwe, now dehorn their rhinos in an attempt to reduce poaching. Historical data collected by the Namibian Ministry of Environment, Forestry, and Tourism and Save the Rhino Trust were used to investigate whether dehorning has an effect on age at first reproduction (AFR), inter-calving interval (ICI), birth sex ratios, calf survival, cause of death, and lifespan. These were assessed in four Namibian sub-populations (hereafter referred to as ‘populations’) of black rhino (denoted A, B, C, and D) which have undergone varying levels of dehorning. No significant difference was found in any of the variables between dehorned and horned individuals. Population was a significant predictor of AFR (LRT = 7.433, p = 0.024) and ICI (LRT = 13.281, p = 0.004), although pairwise comparisons only found populations A and B to be significantly different (AFR: z = −2.736, p = 0.017, ICI: z = −3.705, p = 0.001). Additionally, a significantly higher number of males than females were born in population D (p = 0.021, CI = 0.544, 0.960). The main cause of death across all individuals was poaching, although there was no significant difference in the proportion of rhinos poached between dehorned and horned individuals (X² = 0.638, p = 0.424, n = 265). No evidence was found to suggest that dehorning has any effect on AFR, ICI, birth sex ratios, calf survival, or lifespan in the black rhino, which is reassuring in the continued use of dehorning as an anti-poaching technique in this species.

Female reproductive skew exacerbates the extinction risk from poaching in the eastern black rhino

Variation in individual demographic rates can have large consequences for populations. Female reproductive skew is an example of structured demographic heterogeneity where females have intrinsic qualities that make them more or less likely to breed. The consequences of reproductive skew for population dynamics are poorly understood in non-cooperatively breeding mammals, especially when coupled with other drivers such as poaching. We address this knowledge gap with population viability analyses using an age-specific, female-only, individual-based, stochastic population model built with long-term data for three Kenyan populations of the Critically Endangered eastern black rhino (Diceros bicornis michaeli). There was substantial reproductive skew, with a high proportion of females not breeding or doing so at very low rates. This had a large impact on the projected population growth rate for the smaller population on Ol Jogi. Moreover, including female reproductive skew exacerbates the effects of poaching, increasing the probability of extinction by approximately 70% under a simulated poaching pressure of 5% offtake per year. Tackling the effects of reproductive skew depends on whether it is mediated by habitat or social factors, with potential strategies including habitat and biological management respectively. Investigating and tackling reproductive skew in other species requires long-term, individual-level data collection.

Female Southern White Rhinoceros Can Select Mates to Avoid Inbreeding

Current management models for many endangered species focus primarily on demographic recovery, often ignoring their intrinsic ecological requirements. Across the protected area network of southern Africa, most southern white rhinoceros are managed in populations of less than 50 individuals, experiencing restricted dispersal opportunities, and limited breeding male numbers due to their exclusive home range requirements. In the absence of information on the breeding structure of these populations, poor management decisions may require females to either forego a breeding opportunity or select to inbreed with close relatives. Here, we use a combination of social pedigree data together with genetic analyses to reconstruct the parentage of all 28 offspring produced in a 5-year period in a managed free-ranging southern white rhinoceros population. During this period, all breeding females (founders and first-generation daughters) had access to both a founder male (father to most of the daughters) and two recently introduced inexperienced males. We report that while founder females were more likely to breed with the founder male, their daughters, in contrast, were more likely to breed with the introduced males, thus avoiding inbreeding. However, we also found evidence of father-daughter inbreeding in this population, and contend that in the absence of choice, rather than forego a breeding opportunity, female white rhinoceros will inbreed with their fathers. We argue that to effectively conserve the southern white rhinoceros, managers need to understand the breeding structure of these small populations, particularly in terms of parentage and kinship.

Using Reciprocal Transplants to Assess Local Adaptation, Genetic Rescue, and Sexual Selection in Newly Established Populations

Small populations establishing on colonization fronts have to adapt to novel environments with limited genetic variation. The pace at which they can adapt, and the influence of genetic variation on their success, are key questions for understanding intraspecific diversity. To investigate these topics, we performed a reciprocal transplant experiment between two recently founded populations of brown trout in the sub-Antarctic Kerguelen Islands. Using individual tagging and genetic assignment methods, we tracked the fitness of local and foreign individuals, as well as the fitness of their offspring over two generations. In both populations, although not to the same extent, gene flow occurred between local and foreign gene pools. In both cases, however, we failed to detect obvious footprints of local adaptation (which should limit gene flow) and only weak support for genetic rescue (which should enhance gene flow). In the population where gene flow from foreign individuals was low, no clear differences were observed between the fitness of local, foreign, and F1 hybrid individuals. In the population where gene flow was high, foreign individuals were successful due to high mating success rather than high survival, and F1 hybrids had the same fitness as pure local offspring. These results suggest the importance of considering sexual selection, rather than just local adaptation and genetic rescue, when evaluating the determinants of success in small and recently founded populations.

Sexual selection on male body size, genital length and heterozygosity: Consistency across habitats and social settings

Spatial and temporal variation in environmental factors and the social setting can help to maintain genetic variation in sexually selected traits if it affects the strength of directional selection. A key social parameter which affects the intensity of, and sometimes predicts the response to, mating competition is the operational sex ratio (OSR; ratio of receptive males to females). How the OSR affects selection for specific male traits is poorly understood. It is also unclear how sexual selection is affected by interactions between the OSR and environmental factors, such as habitat complexity, that alter key male-female interactions such as mate encounter rates. Here, we experimentally manipulated the OSR and habitat complexity and quantified sexual selection on male mosquitofish (Gambusia holbrooki) by directly measuring male reproductive success (i.e. paternity). We show that despite a more equitable sharing of paternity (i.e. higher levels of multiple paternity) under a male-biased OSR, selection on focal male traits was unaffected by the OSR or habitat complexity. Instead, sexual selection consistently, and significantly, favoured smaller bodied males, males with higher genome wide heterozygosity (based on >3,000 SNP markers) and males with a relatively long gonopodium (intromittent organ). Our results show that sexual selection on male body size, relative genital size and heterozygosity in this system is consistent across environments that vary in ecological parameters that are expected to influence mate encounter rates.

Population demography and heterozygosity-fitness correlations in natural guppy populations: An examination using sexually selected fitness traits

Heterozygosity-fitness correlations (HFCs) have been examined in a wide diversity of contexts, and the results are often used to infer the role of inbreeding in natural populations. Although population demography, reflected in population-level genetic parameters such as allelic diversity or identity disequilibrium, is expected to play a role in the emergence and detectability of HFCs, direct comparisons of variation in HFCs across many populations of the same species, with different genetic histories, are rare. Here, we examined the relationship between individual microsatellite heterozygosity and a range of sexually selected traits in 660 male guppies from 22 natural populations in Trinidad. Similar to previous studies, observed HFCs were weak overall. However, variation in HFCs among populations was high for some traits (although these variances were not statistically different from zero). Population-level genetic parameters, specifically genetic diversity levels (number of alleles, observed/expected heterozygosity) and measures of identity disequilibrium (g2 and heterozygosity-heterozygosity correlations), were not associated with variation in population-level HFCs. This latter result indicates that these metrics do not necessarily provide a reliable predictor of HFC effect sizes across populations. Importantly, diversity and identity disequilibrium statistics were not correlated, providing empirical evidence that these metrics capture different essential characteristics of populations. A complex genetic architecture likely underpins multiple fitness traits, including those associated with male fitness, which may have reduced our ability to detect HFCs in guppy populations. Further advances in this field would benefit from additional research to determine the demographic contexts in which HFCs are most likely to occur.

Reproductive senescence: new perspectives in the wild

According to recent empirical studies, reproductive senescence, the decline in reproductive success with increasing age, seems to be nearly ubiquitous in the wild. However, a clear understanding of the evolutionary causes and consequences of reproductive senescence is still lacking and requires new and integrative approaches. After identifying the sequential and complex nature of female reproductive senescence, we show that the relative contributions of physiological decline and alterations in the efficiency of parental care to reproductive senescence remain unknown and need to be assessed in the light of current evolutionary theories of ageing. We demonstrate that, although reproductive senescence is generally studied only from the female viewpoint, age-specific female reproductive success strongly depends on male-female interactions. Thus, a reduction in male fertilization efficiency with increasing age has detrimental consequences for female fitness. Lastly, we call for investigations of the role of environmental conditions on reproductive senescence, which could provide salient insights into the underlying sex-specific mechanisms of reproductive success. We suggest that embracing such directions should allow building new bridges between reproductive senescence and the study of sperm competition, parental care, mate choice and environmental conditions.

Adult survival selection in relation to multilocus heterozygosity and body size in a tropical bird species, the Zenaida dove, Zenaida aurita

Both phenotypic and genetic quality can influence the survival of individuals through time, although their relative influences are rarely addressed simultaneously. Here we used capture-mark-recapture modelling to assess the influence of both multilocus heterozygosity (MLH) and body size on apparent adult survival in a tropical bird species, the Zenaida dove, Zenaida aurita, using a sample of 391 individuals genotyped at 11 microsatellites, while controlling for the effects of sex. No effect of body size on either adult survival or capture rate was found. In the best model, survival was a logit linear function of MLH, whereas detection probability was a sex-dependent logit linear function of the logarithm of field effort, increasing with time and affected by a random individual effect. Using a Bayesian approach, we found that MLH explained 1.14% of the total deviance, as expected from theory and previous studies of heterozygosity-fitness correlations, with no evidence for local effects. However, results from capture-mark-recapture modelling indicated that expected longevity varied from 4.8 years in the least heterozygous individuals (MLH = 0.37) to 10.6 years in the most heterozygous ones (MLH = 1), thus suggesting that MLH had potentially a substantial effect on survival. We discuss our results in relation to current hypotheses about the origin of heterozygosity-fitness correlations.

Male reproductive success is correlated with testosterone in the eastern black rhinoceros (Diceros bicornis michaeli)

Among natural populations of polygynous species, males often vary in their lifetime reproductive success. However, in managed populations of endangered species, either in situ or as part of captive breeding programmes, it is important to understand why differences in reproductive success occur. The European captive population of the critically endangered eastern black rhinoceros is currently under-performing relative to their wild counterparts, with low reproductive output and high reproductive skew limiting growth and genetic diversity. To investigate why over 40% of captive males fail to breed, faecal samples were collected weekly from 23 males at 12 institutions across Europe for 4-32 months. Testosterone metabolite concentration was compared between proven and non-proven males and a number of intrinsic and extrinsic factors that could influence reproductive success were also investigated. Males that sired within the last 3½ years had significantly higher androgen concentrations than non-proven males, and average testosterone was positively correlated with the number of offspring sired per year spent in the reproductive age class. Proven and non-proven males did not differ in their body condition, or in average faecal glucocorticoid concentration. Differences in individual temperament were associated with adrenal activity, but did not correlate with reproductive category. Highest testosterone concentrations were observed in proven males that were housed with females during oestrus, and lowest concentrations in non-proven females not housed with females at all during the study period. Further work is necessary to determine whether proven males had higher testosterone due to underlying differences associated with quality, or whether external stimuli such as access to females could influence testosterone concentration and increase a male's chances of becoming a successful breeder.