Friends and family: A software program for identification of unrelated individuals from molecular marker data

The identification of related and unrelated individuals from molecular marker data is often difficult, particularly when no pedigree information is available and the data set is large. High levels of relatedness or inbreeding can influence genotype frequencies and thus genetic marker evaluation, as well as the accurate inference of hidden genetic structure. Identification of related and unrelated individuals is also important in breeding programmes, to inform decisions about breeding pairs and translocations. We present Friends and Family, a Windows executable program with a graphical user interface that identifies unrelated individuals from a pairwise relatedness matrix or table generated in programs such as coancestry and genalex . Friends and Family outputs a list of samples that are all unrelated to each other, based on a user‐defined relatedness cut‐off value. This unrelated data set can be used in downstream analyses, such as marker evaluation or inference of genetic structure. The results can be compared to that of the full data set to determine the effect related individuals have on the analyses. We demonstrate one of the applications of the program: how the removal of related individuals altered the Hardy–Weinberg equilibrium test outcome for microsatellite markers in an empirical data set. Friends and Family can be obtained from https://github.com/DeondeJager/Friends-and-Family .     

Estimating the inbreeding level and genetic relatedness in an isolated population of critically endangered Sichuan taimen (Hucho Bleekeri) using genome‐wide SNP markers

Sichuan taimen (Hucho bleekeri) is critically endangered fish listed in The Red List of Threatened Species compiled by the International Union for Conservation of Nature (IUCN). Specific locus amplified fragment sequencing (SLAF‐seq)‐based genotyping was performed for Sichuan taimen with 43 yearling individuals from three locations in Taibai River (a tributary of Yangtze River) that has been sequestered from its access to the ocean for more than 30 years since late 1980s. Applying the inbreeding level and genetic relatedness estimation using 15,396 genome‐wide SNP markers, we found that the inbreeding level of this whole isolated population was at a low level (2.6 × 10^?3 ± 0.079), and the means of coancestry coefficients within and between the three sampling locations were all very low (close to 0), too. Genomic differentiation was negatively correlated with the geographical distances between the sampling locations (p < .001), and the 43 individuals could be considered as genetically independent two groups. The low levels of genomic inbreeding and relatedness indicated a relatively large number of sexually mature individuals were involved in reproduction in Taibai River. This study suggested a genomic‐relatedness‐guided breeding and conservation strategy for wild fish species without pedigree information records.     

When genes move farther than offspring: gene flow by male gamete dispersal in the highly philopatric bat species Thyroptera tricolor

For species characterized by philopatry of both sexes, mate selection represents an important behaviour for inbreeding avoidance, yet the implications for gene flow are rarely quantified. Here, we present evidence of male gamete‐mediated gene flow resulting from extra‐group mating in Spix's disc‐winged bat, Thyroptera tricolor, a species which demonstrates all‐offspring philopatry. We used microsatellite and capture–recapture data to characterize social group structure and the distribution of mated pairs at two sites in southwestern Costa Rica over four breeding seasons. Relatedness and genetic spatial autocorrelation analyses indicated strong kinship within groups and over short distances (<50 m), resulting from matrilineal group structure and small roosting home ranges (~0.2 ha). Despite high relatedness among‐group members, observed inbreeding coefficients were low (F_IS = 0.010 and 0.037). Parentage analysis indicated mothers and offspring belonged to the same social group, while fathers belonged to different groups, separated by large distances (~500 m) when compared to roosting home ranges. Simulated random mating indicated mate choice was not based on intermediate levels of relatedness, and mated pairs were less related than adults within social groups on average. Isolation‐by‐distance (IBD) models of genetic neighbourhood area based on father–offspring distances provided direct estimates of mean gamete dispersal distances () > 10 roosting home range equivalents. Indirect estimates based on genetic distance provided even larger estimates of , indicating direct estimates were biased low. These results suggest extra‐group mating reduces the incidence of inbreeding in T. tricolor, and male gamete dispersal facilitates gene flow in lieu of natal dispersal of young.     

High‐throughput SNPs for all: genotyping‐in‐thousands

Understanding the genetic structure of species is essential for conservation. It is only with this information that managers, academics, user groups and land‐use planners can understand the spatial scale of migration and local adaptation, source‐sink dynamics and effective population size. Such information is essential for a multitude of applications including delineating management units, balancing management priorities, discovering cryptic species and implementing captive breeding programmes. Species can range from locally adapted by hundreds of metres (Pavey et al. 2010 ) to complete species panmixia (Côté et al. 2013 ). Even more remarkable is that this essential information can be obtained without fully sequenced or annotated genomes, but from mere (putatively) nonfunctional variants. First with allozymes, then microsatellites and now SNPs, this neutral genetic variation carries a wealth of information about migration and drift. For many of us, it may be somewhat difficult to remember our understanding of species conservation before the widespread usage of these useful tools. However most species on earth have yet to give us that ‘peek under the curtain’. With the current diversity on earth estimated to be nearly 9 million species (Mora et al. 2011 ), we have a long way to go for a comprehensive meta‐phylogeographic understanding. A method presented in this issue by Campbell and colleagues (Campbell et al. 2015 ) is a tool that will accelerate the pace in this area. Genotyping‐in‐thousands (GT‐seq) leverages recent advancements in sequencing technology to save many hours and dollars over previous methods to generate this important neutral genetic information.     

Fine‐scale spatiotemporal patterns of genetic variation reflect budding dispersal coupled with strong natal philopatry in a cooperatively breeding mammal

The relatedness structure of animal populations is thought to be a critically important factor underlying the evolution of mating systems and social behaviours. While previous work has shown that population structure is shaped by many biological processes, few studies have investigated how these factors vary over time. Consequently, we explored the fine‐scale spatiotemporal genetic structure of an intensively studied population of cooperatively breeding banded mongooses (Mungos mungo) over a 10‐year period. Overall population structure was strong (average F_ST = 0.129) but groups with spatially overlapping territories were not more genetically similar to one another than noncontiguous groups. Instead, genetic differentiation was associated with historical group‐fission (budding) events, with new groups diverging from their parent groups over time. Within groups, relatedness was high within but not between the sexes, although the latter increased over time since group formation due to group founders being replaced by philopatric young. This trend was not mirrored by a decrease in average offspring heterozygosity over time, suggesting that close inbreeding may often be avoided, even when immigration into established groups is virtually absent and opportunities for extra‐group matings are rare. Fine‐scale spatiotemporal population structure could have important implications in social species, where relatedness between interacting individuals is a vital component in the evolution of patterns of inbreeding avoidance, reproductive skew and kin‐selected helping and harming.     

Sampling related individuals within ponds biases estimates of population structure in a pond‐breeding amphibian

Effective conservation and management of pond‐breeding amphibians depends on the accurate estimation of population structure, demographic parameters, and the influence of landscape features on breeding‐site connectivity. Population‐level studies of pond‐breeding amphibians typically sample larval life stages because they are easily captured and can be sampled nondestructively. These studies often identify high levels of relatedness between individuals from the same pond, which can be exacerbated by sampling the larval stage. Yet, the effect of these related individuals on population genetic studies using genomic data is not yet fully understood. Here, we assess the effect of within‐pond relatedness on population and landscape genetic analyses by focusing on the barred tiger salamanders (Ambystoma mavortium) from the Nebraska Sandhills. Utilizing genome‐wide SNPs generated using a double‐digest RADseq approach, we conducted standard population and landscape genetic analyses using datasets with and without siblings. We found that reduced sample sizes influenced parameter estimates more than the inclusion of siblings, but that within‐pond relatedness led to the inference of spurious population structure when analyses depended on allele frequencies. Our landscape genetic analyses also supported different models across datasets depending on the spatial resolution analyzed. We recommend that future studies not only test for relatedness among larval samples but also remove siblings before conducting population or landscape genetic analyses. We also recommend alternative sampling strategies to reduce sampling siblings before sequencing takes place. Biases introduced by unknowingly including siblings can have significant implications for population and landscape genetic analyses, and in turn, for species conservation strategies and outcomes.     

metapop2: Re‐implementation of software for the analysis and management of subdivided populations using gene and allelic diversity

Management programmes often have to make decisions based on the analysis of the genetic properties and diversity of populations. Expected heterozygosity (or gene diversity) and population structure parameters are often used to make recommendations for conservation, such as avoidance of inbreeding or migration across subpopulations. Allelic diversity, however, can also provide complementary and useful information for conservation programmes, as it is highly sensitive to population bottlenecks, and is more related to long‐term selection response than heterozygosity. Here we present a completely revised and updated re‐implementation of the software metapop for the analysis of diversity in subdivided populations, as well as a tool for the management and dynamic estimation of optimal contributions in conservation programmes. This new update includes computation of allelic diversity for population analysis and management, as well as a simulation mode to forecast the consequences of taking different management strategies over time. Furthermore, the new implementation in C++ includes code optimization and improved memory usage, allowing for fast analysis of large data sets including single nucleotide polymorphism markers, as well as enhanced cross‐software and cross‐platform compatibility.     

Ex situ population management in the absence of pedigree information

For captive breeding to play a significant role in conservation, ex situ populations must be scientifically managed to meet objective goals for retaining representative genetic variation. Imperfect genealogical information requires fundamental assumptions to be made that may bias downstream measures of genetic importance, upon which management decisions are based. The impacts of such assumptions are most pronounced within breeding programmes characterized by a high proportion of individuals of unknown ancestry, as exemplified by the large captive population of the St Vincent parrot (Amazona guildingii). The degree to which microsatellite-based estimates of relatedness may improve upon the assumptions of conventional pedigree-based management was investigated using genotypic data collected at eight microsatellite loci and two marker-based relatedness estimators. The measure, rxyLR, was found to explain the highest amount of variation in true relatedness. Integration of pairwise estimates of founder relatedness with studbook data transformed current understanding of the relatedness structure of the A. guildingii population from two subgroups characterized by a high and low degree of relatedness, respectively, to a situation where all 72 individuals are prioritized for breeding according to their estimated mean kinships. Furthermore, the discovery of opposing, directional bias exhibited by rxyLR and rxyQG in assigning dyads to a given relationship category suggests that an approach that utilizes a combination of pairwise relatedness estimators may provide the most genetic information for balancing the dual considerations of maximizing gene diversity and minimizing inbreeding in developing breeding recommendations.     

Pedigree analysis for the genetic management of group‐living species

Captive breeding programs are an important tool for the conservation of endangered species. These programs are commonly managed using pedigrees containing information about the history of each individual's family, such as breeding pairs and parentage. However, there are some species that are kept in groups where it is hard to distinguish between particular individuals within the group, making it very difficult to record any information at an individual level. Currently, software and methods commonly used for registering and analyzing pedigrees to help manage populations at an individual level are not adequate for managing these group‐living species. Therefore, there is a need to further develop these tools and methodologies for pedigree analysis to better manage group‐living species. PMx is a program used for the management of ex situ populations in zoos and aquariums. We adapted the pedigree analysis method implemented in PMx to analyze pedigrees (records of descendant lineages) of group‐living species. In addition, we developed a group pedigree data entry sheet and group2PMx, a converter program that enables group datasets to be imported into PMx. We show how pedigree analysis of a group‐living species can be used for population management using the studbook of the endangered Texas blind cave salamander Eurycea rathbuni. Such analyses of the pedigree of groups can improve the management of group‐living species in ex situ breeding programs. Firstly, it enables better management decisions based on more accurate genetic measures between groups, allowing for greater control of inbreeding. Secondly, it can improve the conditions in which group‐living species are held by adapting husbandry practices to better reflect conditions of these species living in the wild. The use of the spreadsheet and group2PMx extends the application of PMx_, allowing conservation managers and other institutions outside the zoo and aquarium community to easily import and analyze their pedigree data.     

Close‐kin mating,but not inbred parents,reduces hatching rates and offspring quality in a threatened tortoise

Inbreeding depression, the reduction in fitness due to mating of related individuals, is of particular conservation concern in species with small, isolated populations. Although inbreeding depression is widespread in natural populations, long‐lived species may be buffered from its effects during population declines due to long generation times and thus are less likely to have evolved mechanisms of inbreeding avoidance than species with shorter generation times. However, empirical evidence of the consequences of inbreeding in threatened, long‐lived species is limited. In this study, we leverage a well‐studied population of gopher tortoises, Gopherus polyphemus, to examine the role of inbreeding depression and the potential for behavioural inbreeding avoidance in a natural population of a long‐lived species. We tested the hypothesis that increased parental inbreeding leads to reduced hatching rates and offspring quality. Additionally, we tested for evidence of inbreeding avoidance. We found that high parental relatedness results in offspring with lower quality and that high parental relatedness is correlated with reduced hatching success. However, we found that hatching success and offspring quality increase with maternal inbreeding, likely due to highly inbred females mating with more distantly related males. We did not find evidence for inbreeding avoidance in males and outbred females, suggesting sex‐specific evolutionary trade‐offs may have driven the evolution of mating behaviour. Our results demonstrate inbreeding depression in a long‐lived species and that the evolution of inbreeding avoidance is shaped by multiple selective forces.     

Gene dispersal and outbreeding in a philopatric mammal

Extensive mark-recapture data from banner-tailed kangaroo rats, Dipodomys spectabilis, have shown that both males and females are highly philopatric and suggest the possibility of close inbreeding. However, indirect analyses based on genetic structure appear to contradict direct observations, suggesting longer dispersal distances. Using microsatellite genotypes from most members of a banner-tailed kangaroo rat population during five successive breeding seasons, we ask how relatedness is influenced by dispersal and how it in turn influences mating patterns. The data confirm that, because of philopatry, neighbours are often close relatives. However, patterns of parentage also show that the average distance between mates is large relative to natal dispersal distances and larger than the average distance between nearest opposite-sexed neighbours. Females' mates were often not their nearest male neighbour and many were less related than the nearest male neighbour. We detected multiple paternity in some females' litters; both sexes produce offspring with multiple mates within and between breeding seasons. At the population level, heterozygosities were high and estimates of F were low, indicating that levels of inbreeding were low. Using individual inbreeding coefficients of all juveniles to estimate their parents' relatedness, we found that parental relatedness was significantly lower than relatedness between nearest opposite-sexed adult neighbours. Thus in philopatric populations, long breeding forays can cause genes to move further than individuals disperse, and polyandry may serve to reduce relatedness between mates.     

Rhh: an R extension for estimating multilocus heterozygosity and heterozygosity-heterozygosity correlation

Individual multilocus heterozygosity estimates based on a limited number of loci are expected to correlate only weakly with the inbreeding level of an individual. Before using multilocus heterozygosity estimates in studies of inbreeding, their ability to capture information on inbreeding in the given setting should be tested. A convenient method for this is to compute the heterozygosity-heterozygosity correlation, i.e. the mean correlation between multilocus heterozygosity estimates calculated from random samples of loci, which should be positive if multilocus heterozygosity carries a signature of inbreeding. Rhh is an extension package for the statistical software r that estimates this correlation and calculates three measures of individual multilocus heterozygosity: homozygosity by loci, internal relatedness and standardized heterozygosity. The extension package is available through the CRAN (http://cran.r-project.org) and has a homepage at http://www.helsinki.fi/biosci/egru/research/software.     

Evaluating ex situ conservation projects: Genetic structure of the captive population of the Arabian sand cat

Ex situ conservation plays an increasingly important role in the conservation of endangered species. Molecular genetic markers can be helpful to assess the status of captive breeding programmes. We present the first molecular genetic analysis of the captive population of the Arabian sand cat (Felis margarita harrisoni) using microsatellites. Our data indicates that the captive population of F. m. harrisoni comprises three genetic clusters, which are based on different founder lineages. Genetic diversity was relatively high, the effective population size even exceeded the number of founders. This was presumably caused by subsequently integrating unrelated, genetically diverse founders into the captive population and a careful management based on minimizing kinship. However, we detected an error in the studbook records, which might have led to incestuous matings and underlines the usefulness of molecular evaluations in captive breeding programmes for endangered species.     

How many came home? Evaluating ex situ conservation of green turtles in the Cayman Islands

Ex situ management is an important conservation tool that allows the preservation of biological diversity outside natural habitats while supporting survival in the wild. Captive breeding followed by re‐introduction is a possible approach for endangered species conservation and preservation of genetic variability. The Cayman Turtle Centre Ltd was established in 1968 to market green turtle (Chelonia mydas) meat and other products and replenish wild populations, thought to be locally extirpated, through captive breeding. We evaluated the effects of this re‐introduction programmme using molecular markers (13 microsatellites, 800‐bp D‐loop and simple tandem repeat mitochondrial DNA sequences) from captive breeders (N = 257) and wild nesting females (N = 57) (sampling period: 2013–2015). We divided the captive breeders into three groups: founders (from the original stock), and then two subdivisions of F₁ individuals corresponding to two different management strategies, cohort 1995 (“C1995”) and multicohort F₁ (“MCF1”). Loss of genetic variability and increased relatedness was observed in the captive stock over time. We found no significant differences in diversity among captive and wild groups, and similar or higher levels of haplotype variability when compared to other natural populations. Using parentage and sibship assignment, we determined that 90% of the wild individuals were related to the captive stock. Our results suggest a strong impact of the re‐introduction programmme on the present recovery of the wild green turtle population nesting in the Cayman Islands. Moreover, genetic relatedness analyses of captive populations are necessary to improve future management actions to maintain genetic diversity in the long term and avoid inbreeding depression.     

Population structure of a vector‐borne plant parasite

Parasites are among the most diverse groups of life on Earth, yet complex natural histories often preclude studies of their speciation processes. The biology of parasitic plants facilitates in situ collection of data on both genetic structure and the mechanisms responsible for that structure. Here, we studied the role of mating, dispersal and establishment in host race formation of a parasitic plant. We investigated the population genetics of a vector‐borne desert mistletoe (Phoradendron californicum) across two legume host tree species (Senegalia greggii and Prosopis velutina) in the Sonoran desert using microsatellites. Consistent with host race formation, we found strong host‐associated genetic structure in sympatry, little genetic variation due to geographic site and weak isolation by distance. We hypothesize that genetic differentiation results from differences in the timing of mistletoe flowering by host species, as we found initial flowering date of individual mistletoes correlated with genetic ancestry. Hybrids with intermediate ancestry were detected genetically. Individuals likely resulting from recent, successful establishment events following dispersal between the host species were detected at frequencies similar to hybrids between host races. Therefore, barriers to gene flow between the host races may have been stronger at mating than at dispersal. We also found higher inbreeding and within‐host individual relatedness values for mistletoes on the more rare and isolated host species (S. greggii). Our study spanned spatial scales to address how interactions with both vectors and hosts influence parasitic plant structure with implications for parasite virulence evolution and speciation.     

Banded mongooses avoid inbreeding when mating with members of the same natal group

Inbreeding and inbreeding avoidance are key factors in the evolution of animal societies, influencing dispersal and reproductive strategies which can affect relatedness structure and helping behaviours. In cooperative breeding systems, individuals typically avoid inbreeding through reproductive restraint and/or dispersing to breed outside their natal group. However, where groups contain multiple potential mates of varying relatedness, strategies of kin recognition and mate choice may be favoured. Here, we investigate male mate choice and female control of paternity in the banded mongoose (Mungos mungo), a cooperatively breeding mammal where both sexes are often philopatric and mating between relatives is known to occur. We find evidence suggestive of inbreeding depression in banded mongooses, indicating a benefit to avoiding breeding with relatives. Successfully breeding pairs were less related than expected under random mating, which appeared to be driven by both male choice and female control of paternity. Male banded mongooses actively guard females to gain access to mating opportunities, and this guarding behaviour is preferentially directed towards less closely related females. Guard–female relatedness did not affect the guard's probability of gaining reproductive success. However, where mate‐guards are unsuccessful, they lose paternity to males that are less related to the females than themselves. Together, our results suggest that both sexes of banded mongoose use kin discrimination to avoid inbreeding. Although this strategy appears to be rare among cooperative breeders, it may be more prominent in species where relatedness to potential mates is variable, and/or where opportunities for dispersal and mating outside of the group are limited.     

Estimating relatedness and inbreeding using molecular markers and pedigrees: the effect of demographic history

Estimates of inbreeding and relatedness are commonly calculated using molecular markers, although the accuracy of such estimates has been questioned. As a further complication, in many situations, such estimates are required in populations with reduced genetic diversity, which is likely to affect their accuracy. We investigated the correlation between microsatellite‐ and pedigree‐based coefficients of inbreeding and relatedness in laboratory populations of Drosophila melanogaster that had passed through bottlenecks to manipulate their genetic diversity. We also used simulations to predict expected correlations between marker‐ and pedigree‐based estimates and to investigate the influence of linkage between loci and null alleles. Our empirical data showed lower correlations between marker‐ and pedigree‐based estimates in our control (nonbottleneck) population than were predicted by our simulations or those found in similar studies. Correlations were weaker in bottleneck populations, confirming that extreme reductions in diversity can compromise the ability of molecular estimates to detect recent inbreeding events. However, this result was highly dependent on the strength of the bottleneck and we did not observe or predict any reduction in correlations in our population that went through a relatively severe bottleneck of N = 10 for one generation. Our results are therefore encouraging, as molecular estimates appeared robust to quite severe reductions in genetic diversity. It should also be remembered that pedigree‐based estimates may not capture realized identity‐by‐decent and that marker‐based estimates may actually be more useful in certain situations.     

On the use of large marker panels to estimate inbreeding and relatedness: empirical and simulation studies of a pedigreed zebra finch population typed at 771 SNPs

In recent years there has been a dramatic increase in the availability of high density genetic marker data for both model and non‐model organisms. A potential application of these data is to infer relatedness in the absence of a complete pedigree. Using a marker panel of 771 SNPs genotyped in three generations of an extensive zebra finch pedigree, correlations between pedigree relatedness and seven marker‐based estimates of relatedness were examined, as was the relationship between heterozygosity and inbreeding. Although marker‐based and pedigree relatedness were highly correlated, the variance in estimated relatedness was high. Further, the correlation between heterozygosity and inbreeding was weak, even though mean inbreeding coefficient is typical of that seen in wild vertebrate pedigrees; the weak relationship was in part due to the small variance in inbreeding in the pedigree. Our data suggest that using marker information to reconstruct the pedigree, and then calculating relatedness from the pedigree, is likely to give more accurate relatedness estimates than using marker‐based estimators directly.     

Molecular and pedigree measures of relatedness provide similar estimates of inbreeding depression in a bottlenecked population

Individual‐based estimates of the degree of inbreeding or parental relatedness from pedigrees provide a critical starting point for studies of inbreeding depression, but in practice wild pedigrees are difficult to obtain. Because inbreeding increases the proportion of genomewide loci that are identical by descent, inbreeding variation within populations has the potential to generate observable correlations between heterozygosity measured using molecular markers and a variety of fitness related traits. Termed heterozygosity‐fitness correlations (HFCs), these correlations have been observed in a wide variety of taxa. The difficulty of obtaining wild pedigree data, however, means that empirical investigations of how pedigree inbreeding influences HFCs are rare. Here, we assess evidence for inbreeding depression in three life‐history traits (hatching and fledging success and juvenile survival) in an isolated population of Stewart Island robins using both pedigree‐ and molecular‐derived measures of relatedness. We found results from the two measures were highly correlated and supported evidence for significant but weak inbreeding depression. However, standardized effect sizes for inbreeding depression based on the pedigree‐based kin coefficients (k) were greater and had smaller standard errors than those based on molecular genetic measures of relatedness (RI), particularly for hatching and fledging success. Nevertheless, the results presented here support the use of molecular‐based measures of relatedness in bottlenecked populations when information regarding inbreeding depression is desired but pedigree data on relatedness are unavailable.     

Using relatedness networks to infer contemporary dispersal: Application to the endangered mammal Galemys pyrenaicus

Information about the degree of contemporary dispersal is important when trying to understand how populations interchange individuals and identify the specific barriers that prevent these movements. In the case of endangered species, this can represent crucial information when designing appropriate conservation strategies. Here we analyse relatedness between individuals from different localities and use these data to infer whether dispersal occurred in recent generations. We applied this approach to the Pyrenean desman (Galemys pyrenaicus), a semiaquatic and endangered species endemic to the Iberian Peninsula. We studied this species in four primary rivers in the Iberian Range, where two ancient mitochondrial lineages are separated by a strict contact zone, suggesting the existence of complex dispersal patterns. Using next‐generation sequencing, we obtained 912 SNPs from each specimen and estimated relatedness values between them. While relatedness networks were dense within each river, we found surprisingly few relationships between individuals from different rivers despite their close proximity in some cases, indicating much lower dispersal between rivers compared to dispersal within a single river. In agreement with this result, the degree of inbreeding was exceedingly high in most individuals. These data show that relatedness information can be crucial to understand the contemporary dispersal patterns and conservation status of specific populations of endangered species.