Pedigrees and microsatellites among endangered ungulates: what do they tell us?

Pedigree analysis in the mhorr gazelle (Nanger dama mhorr): Genetic variability evolution of the captive population

Breeding programs have an essential role in the recovery of threatened populations through optimal genetic management and mating strategies. The dama gazelle (Nanger dama) is a North African ungulate listed as critically endangered. The mhorr subspecies is extinct in the wild and currently survives thanks to the creation in 1971 of an ex situ breeding program. The aim of the present study was to assess the evolution of genetic variability in this mhorr gazelle captive population, as well as the mating strategy used in two reference populations studied (Almeria and Europe). The entire pedigree, with 2739 animals, was analyzed to measure demographic characters, pedigree completeness level, probability of gene origin, level of relatedness and genetic structure of the population. The population size has been progressively increasing, with up to 264 individuals alive in Europe at the time of the study. The average number of equivalent complete generations was 5.55. The effective number of founders and ancestors was both 3, and the founder genome equivalent was 1.99. The genetic contributions of the four main ancestors were unbalanced. The average values of inbreeding and average relatedness for the whole pedigree were, respectively, 28.34% and 50.14%. The effective population size was 8.7 by individual increase in inbreeding and 9.8 by individual increase in coancestry. F-statistics evidenced a very small level of population subdivision (F ST = 0.033370). The mating strategy used, based on the minimum coancestry of the individuals, has minimized the losses of genetic variability and helped to balance the genetic contributions between ancestors. The strategy also avoided large subdivisions within the population and the appearance of new bottlenecks. This study shows how pedigree analysis can both be used to determine the genetic variability of the population and to assess the influence of the mating strategy used in the breeding program on such variability.

Genetic assessment of captive breeding program of Indian Pangolin: implications for conservation and management

BACKGROUND

Captive breeding programs play a vital role in conservation of threatened species, necessitating an understanding of genetic diversity among captive individuals to ensure long-term genetic viability, appropriate mate selection, and successful reintroduction to native habitats.

METHODS AND RESULTS

We did not observe any recent genetic bottleneck, and population showed moderate genetic diversity. The estimated effective population size, representing individuals capable of contributing genetically to future generations, was estimated as 18.6 individuals (11.4-35.1 at 95% CI). Based on the genetic make-up and allelic diversity, we found seventeen pangolins (11 females and 6 males) were genetically unrelated and relatively more potent than others.

CONCLUSION

In this study, we evaluated the captive breeding program of the Indian pangolin population at the Pangolin Conservation Breeding Centre in Nandankanan Zoological Park, Bhubaneswar, Odisha. We highlight the significance of genetic monitoring within the captive population of Indian pangolin for preserving genetic diversity and ensuring the long-term survival of the species. We established the genetic profiles of all 29 pangolins and identified 17 pangolins to be prioritized for enhanced breeding and future zoo exchange programs. We appreciate the zoo authorities for promoting genetic assessment of pangolin for better and more effective monitoring of the captive breeding of the endangered Indian pangolin.

Microsatellite characterization and development of unified STR panel for big cats in captivity: a case study from a Seoul Grand Park Zoo, Republic of Korea

The zoos manage small populations of endangered big cat species like tiger, lion, and leopard for display, research, and conservation breeding. Genetic management of these populations is essential to ensure long term survival and conservation utility. Here we propose a simple and cost effective microsatellite based protocol for the genetic management of captive big cats. We sampled 36 big cat individuals from Seoul Grand Park Zoo (Republic of Korea) and amplified 33 published microsatellite loci. Overall, allelic richness and gene diversity was found highest for leopards, followed by lions and tigers. Twelve of the thirty-three markers showed a high degree of polymorphism across all target species. These microsatellites provide a high degree of discrimination for tiger (1.45 × 10^-8), lion (1.54 × 10^-10), and leopard (1.88 × 10^-12) and thus can be adopted for the genetic characterization of big cats in accredited zoos globally. During captive breeding, zoo authorities rely on pedigree records maintained in studbooks to ensure mating of genetically fit unrelated individuals. Several studies have reported errors in studbook records of big cat species. Microsatellites are simple and cost effective tool for DNA fingerprinting, estimation of genetic diversity, and paternity assessment. Our unified microsatellite panel (12-plex) for big cats is efficient and can easily be adopted by zoo authorities for regular population management.

Genetic assessment of a conservation breeding program of the houbara bustard (Chlamydotis undulata undulata) in Morocco, based on pedigree and molecular analyses

Protection and restoration of species in the wild may require conservation breeding programs under genetic management to minimize deleterious effects of genetic changes that occur in captivity, while preserving populations' genetic diversity and evolutionary resilience. Here, through interannual pedigree analyses, we first assessed the efficiency of a 21-year genetic management, including minimization of mean kinship, inbreeding avoidance, and regular addition of founders, of a conservation breeding program targeting on Houbara bustard (Chlamydotis undulata undulata) in Morocco. Secondly, we compared pedigree analyses, the classical way of assessing and managing genetic diversity in captivity, to molecular analyses based on seven microsatellites. Pedigree-based results indicated an efficient maintenance of the genetic diversity (99% of the initial genetic diversity retained) while molecular-based results indicated an increase in allelic richness and an increase in unbiased expected heterozygosity across time. The pedigree-based average inbreeding coefficient F remained low (between 0.0004 and 0.003 in 2017) while the proportion of highly inbred individuals (F > .1) decreased over time and reached 0.2% in 2017. Furthermore, pedigree-based F and molecular-based individual multilocus heterozygosity were weakly negatively correlated, (Pearson's r = -.061 when considering all genotyped individuals), suggesting that they cannot be considered as alternatives, but rather as complementary sources of information. These findings suggest that a strict genetic monitoring and management, based on both pedigree and molecular tools can help mitigate genetic changes and allow to preserve genetic diversity and evolutionary resilience in conservation breeding programs.

Last exit before the brink: Conservation genomics of the Cambodian population of the critically endangered southern river terrapin

The southern river terrapin, Batagur affinis is one of the world's 25 most endangered freshwater turtle species. The major portion of the global population is currently found in peninsular Malaysia, with the only remnant Indochinese population in southern Cambodia. For more than a decade, wild nests in this remnant Cambodian population have been fenced and hatchlings reared in captivity. Here we amplified 10 microsatellite markers from all 136 captive individuals, obtained 2,658 presumably unlinked and neutral single nucleotide polymorphisms from 72 samples with ddRAD-seq, and amplified 784 bp of mtDNA from 50 samples. Our results reveal that the last Indochinese population comprised only four kinship groups as of 2012, with all offspring sired from <10 individuals in the wild. We demonstrate an obvious decrease in genetic contributions of breeders in the wild from 2006-2012 and identify high-value breeders instrumental for ex-situ management of the contemporary genetic stock of the species.

Effective population sizes and adaptive genetic variation in a captive bird population

Captive populations are considered a key component of ex situ conservation programs. Research on multiple taxa has shown the differential success of maintaining demographic versus genetic stability and viability in captive populations. In typical captive populations, usually founded by few or related individuals, genetic diversity can be lost and inbreeding can accumulate rapidly, calling into question their ultimate utility for release into the wild. Furthermore, domestication selection for survival in captive conditions is another concern. Therefore, it is crucial to understand the dynamics of population sizes, particularly the effective population size, and genetic diversity at non-neutral and adaptive loci in captive populations. In this study, we assessed effective population sizes and genetic variation at both neutral microsatellite markers, as well as SNP variants from the MHC-B locus of a captive Red Junglefowl population. This population represents a rare instance of a population with a well-documented history in captivity, following a realistic scenario of chain-of-custody, unlike many captive lab populations. Our analyses, which included 27 individuals comprising the entirety of one captive population show very low neutral and adaptive genetic variation, as well as low effective sizes, which correspond with the known demographic history. Finally, our study also shows the divergent impacts of small effective size and inbreeding in captive populations on microsatellite versus adaptive genetic variation in the MHC-B locus. Our study provides insights into the difficulties of maintaining adaptive genetic variation in small captive populations.

Low incidence of inbreeding in a long-lived primate population isolated for 75 years

ABSTRACT

When close relatives mate, offspring are expected to suffer fitness consequences due to inbreeding depression. Inbreeding has previously been quantified in two ways: using a sufficiently large panel of markers or deep and complete pedigrees over several generations. However, the application of both approaches is still limited by the challenge of compiling such data for species with long generation times, such as primates. Here, we assess inbreeding in rhesus macaques living on Cayo Santiago (Puerto Rico), a population genetically isolated since 1938, but descendant of a large set of presumably unrelated founders. Using comprehensive genetic data, we calculated inbreeding coefficients (F) for 2669 individuals with complete three generation pedigrees and 609 individuals with complete four generation pedigrees. We found that 0.79 and 7.39% of individuals had an F > 0 when using data from three and four generation pedigrees, respectively. No evidence of an increase in inbreeding over the study period (up to 23 years) was found. Furthermore, the observed mean relatedness of breeding pairs differed significantly from the distribution of parental relatedness expected as simulated based on previous reproductive data, suggesting that kin generally avoid breeding with each other. Finally, inbreeding was not a predictor of early mortality measured as survival until weaning and sexual maturation, respectively. Our results remain consistent with three estimators of inbreeding (standardized heterozygosity, internal relatedness, and homozygosity by loci) using up to 42 highly polymorphic microsatellites for the same set of individuals. Together, our results demonstrate that close inbreeding may not be prevalent even in populations isolated over long periods when mechanisms of inbreeding avoidance can operate.

SIGNIFICANCE STATEMENT

When close relatives mate, offspring may suffer from such inbreeding, e.g., via lower survival and/or fertility. Using (i) a large panel of genetic markers and (ii) complete three or four generation pedigrees, respectively, we show that incidences of inbreeding in a long-lived primate population are rare, even after genetic isolation for 75 years. Moreover, our simulations suggest that kin in our population generally avoid breeding with each other. Finally, the few inbred individuals detected in our large sample did not suffer from lower survival. Given that many animal species face dramatic habitat loss combined with critical population declines, our study provides important implications for conservation biology in general and for population management in particular.

Quantifying realized inbreeding in wild and captive animal populations

Most molecular measures of inbreeding do not measure inbreeding at the scale that is most relevant for understanding inbreeding depression-namely the proportion of the genome that is identical-by-descent (IBD). The inbreeding coefficient FPed obtained from pedigrees is a valuable estimator of IBD, but pedigrees are not always available, and cannot capture inbreeding loops that reach back in time further than the pedigree. We here propose a molecular approach to quantify the realized proportion of the genome that is IBD (propIBD), and we apply this method to a wild and a captive population of zebra finches (Taeniopygia guttata). In each of 948 wild and 1057 captive individuals we analyzed available single-nucleotide polymorphism (SNP) data (260 SNPs) spread over four different genomic regions in each population. This allowed us to determine whether any of these four regions was completely homozygous within an individual, which indicates IBD with high confidence. In the highly nomadic wild population, we did not find a single case of IBD, implying that inbreeding must be extremely rare (propIBD=0-0.00094, 95% CI). In the captive population, a five-generation pedigree strongly underestimated the average amount of realized inbreeding (FPed=0.013<propIBD=0.064), as expected given that pedigree founders were already related. We suggest that this SNP-based technique is generally useful for quantifying inbreeding at the individual or population level, and we show analytically that it can capture inbreeding loops that reach back up to a few hundred generations.

Estimates of direct and indirect effects for early juvenile survival in captive populations maintained for conservation purposes: the case of Cuvier's gazelle

Together with the avoidance of any negative impact of inbreeding, preservation of genetic variability for life-history traits that could undergo future selective pressure is a major issue in endangered species management programmes. However, most of these programmes ignore that, apart from the direct action of genes on such traits, parents, as contributors of offspring environment, can influence offspring performance through indirect parental effects (when parental genotype and phenotype exerts environmental influences on offspring phenotype independently of additive genetic effects). Using quantitative genetic models, we estimated the additive genetic variance for juvenile survival in a population of the endangered Cuvier's gazelle kept in captivity since 1975. The dataset analyzed included performance recording for 700 calves and a total pedigree of 740 individuals. Results indicated that in this population juvenile survival harbors significant additive genetic variance. The estimates of heritability obtained were in general moderate (0.115-0.457) and not affected by the inclusion of inbreeding in the models. Maternal genetic contribution to juvenile survival seems to be of major importance in this gazelle's population as well. Indirect genetic and indirect environmental effects assigned to mothers (i.e., maternal genetic and maternal permanent environmental effects) roughly explain a quarter of the total variance estimated for the trait analyzed. These findings have major evolutionary consequences for the species as show that offspring phenotypes can evolve strictly through changes in the environment provided by mothers. They are also relevant for the captive breeding programme of the species. To take into account, the contribution that mothers have on offspring phenotype through indirect genetic effects when designing pairing strategies might serve to identify those females with better ability to recruit, and, additionally, to predict reliable responses to selection in the captive population.

Recent breeding history of dog breeds in Sweden: modest rates of inbreeding, extensive loss of genetic diversity and lack of correlation between inbreeding and health

One problem in modern dogs is a high occurrence of physical diseases, defects and disorders. Many breeds exhibit physical problems that affect individual dogs throughout life. A potential cause of these problems is inbreeding that is known to reduce the viability of individuals. We investigated the possible correlation between recent inbreeding and health problems in dogs and used studbook data from 26 breeds provided by the Swedish Kennel Club for this purpose. The pedigrees date back to the mid-20th century and comprise 5-10 generations and 1 000-50 000 individuals per pedigree over our study period of 1980-2010. We compared levels of inbreeding and loss of genetic variation measured in relation to the number of founding animals during this period in the investigated dog breeds that we classified as 'healthy' (11 breeds) or 'unhealthy' (15) based on statistics on the extent of veterinary care obtained from Sweden's four largest insurance companies for pets. We found extensive loss of genetic variation and moderate levels of recent inbreeding in all breeds examined, but no strong indication of a difference in these parameters between healthy versus unhealthy breeds over this period. Thus, recent breeding history with respect to rate of inbreeding does not appear to be a main cause of poor health in the investigated dog breeds in Sweden. We identified both strengths and weaknesses of the dog pedigree data important to consider in future work of monitoring and conserving genetic diversity of dog breeds.

Mitochondrial and pedigree analysis in Przewalski's horse populations: implications for genetic management and reintroductions

BACKGROUND AND AIMS

Przewalski's horses have been imported from the western zoos to China since 1985. Yet the genetic diversity in China's populations has not been studied, thus lacking of such knowledge inevitably affects this population's management. The aim of this study was to assess genetic diversity in Chinese population of Przewalski's horses via mitochondrial DNA (mtDNA) control region and pedigree analysis.

MATERIALS AND METHODS

Two captive and one reintroduced populations were examined based on mitochondrial DNA control region variation via fecal sampling from 2010 to 2012, together with pedigree analysis.

RESULTS

Amplification success rates of fecal mtDNA were as high as 96.2% (93.8%-100%), and were higher for sample in winter than in summer and autumn. Two haplotypes were identified and shared among three populations, but the proportion of individuals with each haplotype varied among the three populations (F(ST) = 0.10874, p = 0.00978). Haplotype diversity in the released population (0.153) was much lower than that in the two captive populations (0.4011 and 0.4966), in accordance with the direction of increase in probability of identity at the dam lines.

CONCLUSION

Future concerns in Przewalski's horse population management should emphasize on strict reproduction control to minimize inbreeding in captivity, followed by long-term genetic diversity guidelines and non-invasive monitoring in the reintroduction programmes.

Heterozygosity-fitness correlations and inbreeding depression in two critically endangered mammals

The relation among inbreeding, heterozygosity, and fitness has been studied primarily among outbred populations, and little is known about these phenomena in endangered populations. Most researchers conclude that the relation between coefficient of inbreeding estimated from pedigrees and fitness traits (inbreeding-fitness correlations) better reflects inbreeding depression than the relation between marker heterozygosity and fitness traits (heterozygosity-fitness correlations). However, it has been suggested recently that heterozygosity-fitness correlations should only be expected when inbreeding generates extensive identity disequilibrium (correlations in heterozygosity and homozygosity across loci throughout the genome). We tested this hypothesis in Mohor gazelle (Gazella dama mhorr) and Iberian lynx (Lynx pardinus). For Mohor gazelle, we calculated the inbreeding coefficient and measured heterozygosity at 17 microsatellite loci. For Iberian lynx, we measured heterozygosity at 36 microsatellite loci. In both species we estimated semen quality, a phenotypic trait directly related to fitness that is controlled by many loci and is affected by inbreeding depression. Both species showed evidence of extensive identity disequilibrium, and in both species heterozygosity was associated with semen quality. In the Iberian lynx the low proportion of normal sperm associated with low levels of heterozygosity was so extreme that it is likely to limit the fertility of males. In Mohor gazelle, although heterozygosity was associated with semen quality, inbreeding coefficient was not. This result suggests that when coefficient of inbreeding is calculated on the basis of a genealogy that begins after a long history of inbreeding, the coefficient of inbreeding fails to capture previous demographic information because it is a poor estimator of accumulated individual inbreeding. We conclude that among highly endangered species with extensive identity disequilibrium, examination of heterozygosity-fitness correlations may be an effective way to detect inbreeding depression, whereas inbreeding-fitness correlations may be poor indicators of inbreeding depression if the pedigree does not accurately reflect the history of inbreeding.

The imprecision of heterozygosity-fitness correlations hinders the detection of inbreeding and inbreeding depression in a threatened species

In nonpedigreed wild populations, inbreeding depression is often quantified through the use of heterozygosity-fitness correlations (HFCs), based on molecular estimates of relatedness. Although such correlations are typically interpreted as evidence of inbreeding depression, by assuming that the marker heterozygosity is a proxy for genome-wide heterozygosity, theory predicts that these relationships should be difficult to detect. Until now, the vast majority of empirical research in this area has been performed on generally outbred, nonbottlenecked populations, but differences in population genetic processes may limit extrapolation of results to threatened populations. Here, we present an analysis of HFCs, and their implications for the interpretation of inbreeding, in a free-ranging pedigreed population of a bottlenecked species: the endangered takahe (Porphyrio hochstetteri). Pedigree-based inbreeding depression has already been detected in this species. Using 23 microsatellite loci, we observed only weak evidence of the expected relationship between multilocus heterozygosity and fitness at individual life-history stages (such as survival to hatching and fledging), and parameter estimates were imprecise (had high error). Furthermore, our molecular data set could not accurately predict the inbreeding status of individuals (as 'inbred' or 'outbred', determined from pedigrees), nor could we show that the observed HFCs were the result of genome-wide identity disequilibrium. These results may be attributed to high variance in heterozygosity within inbreeding classes. This study is an empirical example from a free-ranging endangered species, suggesting that even relatively large numbers (>20) of microsatellites may give poor precision for estimating individual genome-wide heterozygosity. We argue that pedigree methods remain the most effective method of quantifying inbreeding in wild populations, particularly those that have gone through severe bottlenecks.

Paternal levels of DNA damage in spermatozoa and maternal parity influence offspring mortality in an endangered ungulate

Understanding which factors influence offspring mortality rates is a major challenge since it influences population dynamics and may constrain the chances of recovery among endangered species. Most studies have focused on the effects of maternal and environmental factors, but little is known about paternal factors. Among most polygynous mammals, males only contribute the haploid genome to their offspring, but the possibility that sperm DNA integrity may influence offspring survival has not been explored. We examined several maternal, paternal and individual factors that may influence offspring survival in an endangered species (Gazella cuvieri). Levels of sperm DNA damage had the largest impact upon offspring mortality rates, followed by maternal parity. In addition, there was a significant interaction between these two variables, so that offspring born to primiparous mothers were more likely to die if their father had high levels of sperm DNA damage, but this was not the case among multiparous mothers. Thus, multiparous mothers seem to protect their offspring from the deleterious effects of sperm DNA damage. Since levels of sperm DNA damage seem to be higher among endangered species, more attention should be paid to the impact of this largely ignored factor among the viability of endangered species.

High levels of DNA fragmentation in spermatozoa are associated with inbreeding and poor sperm quality in endangered ungulates

Inbreeding is known to cause deleterious effects upon reproduction and survival, but its effects upon sperm DNA integrity have not been examined. In the present study, we analyzed this relationship among three endangered ungulates: Gazella cuvieri, Gazella dama mhorr, and Gazella dorcas neglecta. In addition, we examined whether levels of sperm DNA fragmentation are associated with semen quality. The magnitude of sperm DNA damage in the two species with high levels of inbreeding (G. cuvieri and G. dama mhorr) was extremely high when compared to the species with low levels of inbreeding (G. dorcas neglecta) and to values previously reported for outbred populations. Levels of sperm DNA fragmentation significantly increased with inbreeding and age. Increased DNA damage in sperm was associated with increased sperm head abnormalities, lower percentage of sperm with an intact acrosome, and poor motility. Our findings suggest that the link between inbreeding and semen quality is mediated by the effects of inbreeding upon sperm DNA damage. The deleterious effects of inbreeding upon the paternal genome likely decrease male fertility and may cause genetic damage to future generations. Because inbreeding is common among endangered species, high levels of sperm DNA damage may have considerable impact upon the viability of their populations.

Low heterozygosity at microsatellite markers in Iberian red deer with small antlers

Deer antlers are costly structures subjected to directional sexual selection that may be sensitive to heterozygosity. However, a relationship between heterozygosity and antler development has only been found for select protein-coding loci and MHC genes in one deer species (the white-tailed deer Odocoileus virginianus). Here, we study the relationship between multilocus heterozygosity at 11 microsatellite markers and antler size (AS) in a sample of 367 Iberian red deer males (Cervus elaphus hispanicus) from two study areas with different ecological and genetic conditions. We found that males with very small antlers (10% of the sampled individuals with the lowest values of AS) had lower levels of heterozygosity than those with bigger antlers (significant effect in an analysis of variance, P = 0.011). This relationship was noticeable mainly in situations of low genetic diversity, where the differences in heterozygosity between groups of males were greater. Finally, we conducted analyses to address the hypotheses proposed by the heterozygosity-fitness correlation, and we found the local effect as the most likely hypothesis. Our findings reveal an expected but not previously detected association between low heterozygosity and reduced AS, with implications for red deer evolution and management.

Heterozygosity-fitness correlations: a time for reappraisal

Owing to the remarkable progress of molecular techniques, heterozygosity-fitness correlations (HFCs) have become a popular tool to study the impact of inbreeding in natural populations. However, their underlying mechanisms are often hotly debated. Here we argue that these "debates" rely on verbal arguments with no basis in existing theory and inappropriate statistical testing, and that it is time to reconcile HFC with its historical and theoretical fundaments. We show that available data are quantitatively and qualitatively consistent with inbreeding-based theory. HFC can be used to estimate the impact of inbreeding in populations, although such estimates are bound to be imprecise, especially when inbreeding is weak. Contrary to common belief, linkage disequilibrium is not an alternative to inbreeding, but rather comes with some forms of inbreeding, and is not restricted to closely linked loci. Finally, the contribution of local chromosomal effects to HFC, while predicted by inbreeding theory, is expected to be small, and has rarely if ever proven statistically significant using adequate tests. We provide guidelines to safely interpret and quantify HFCs, and present how HFCs can be used to quantify inbreeding load and unravel the structure of natural populations.

Parents' genetic dissimilarity and offspring sex in a polygynous mammal

Offspring quality may benefit from genetic dissimilarity between parents. However, genetic dissimilarity may trade-off with additive genetic benefits. We hypothesized that when sexual selection produces sex-specific selective scenarios, the relative benefits of additive genetic vs. dissimilarity may differ for sons and daughters. Here we study a sample of 666 red deer (Cervus elaphus) microsatellite genotypes, including males, females and their foetuses, from 20 wild populations in Spain (the main analyses are based on 241 different foetuses and 190 mother-foetus pairs). We found that parental lineages were more dissimilar in daughters than in sons. On average, every mother was less related to her mate than to the sample of fathers in the population when producing daughters not sons. Male foetuses conceived early in the rutting season were much more inbred than any other foetuses. These differences maintained through gestation length, ruling out intrauterine mortality as a cause for the results, and indicating that the potential mechanism producing the association between parents' dissimilarity and offspring sex should operate close to mating or conception time. Our findings highlight the relevance of considering the sex of offspring when studying genetic similarity between parents.