DNA-based methods for pedigree reconstruction and kinship analysis in natural populations

MalKinID: A Likelihood-Based Model for Identifying Malaria Parasite Genealogical Relationships Using Identity-by-Descent

Pathogen genomics is a powerful tool for tracking infectious disease transmission. In malaria, identity-by-descent (IBD) is used to assess the genetic relatedness between parasites and has been used to study transmission and importation. In theory, IBD can be used to distinguish genealogical relationships to reconstruct transmission history or identify parasites for genotype-to-phenotype quantitative-trait-locus experiments. MalKinID (Malaria Kinship Identifier) is a new likelihood-based classification model designed to identify genealogical relationships among malaria parasites based on genome-wide IBD proportions and IBD segment distributions. MalKinID was calibrated to the genomic data from three laboratory-based genetic crosses (yielding 440 parent-child and 9060 full-sibling comparisons). MalKinID identified lab generated F1 progeny with >80% sensitivity and showed that 0.39 (95% CI 0.28, 0.49) of the second-generation progeny of a NF54 and NHP4026 cross were F1s and 0.56 (0.45, 0.67) were backcrosses of an F1 with the parental NF54 strain. In simulated outcrossed importations, MalKinID accurately reconstructs genealogy history with high precision and sensitivity, with F1-scores exceeding 0.84. However, when importation involves inbreeding, such as during serial co-transmission, the precision and sensitivity of MalKinID declined, with F1-scores of 0.76 (0.56, 0.92) and 0.23 (0.0, 0.4) for PC and FS and <0.05 for second-degree and third-degree relatives. Genealogical inference is most powered 1) when outcrossing is the norm or 2) when multi-sample comparisons based on a predefined pedigree are used. MalKinID lays the foundations for using IBD to track parasite transmission history and for separating progeny for quantitative-trait-locus experiments.

Distinguishing recrudescence from reinfection in lymphatic filariasis

BACKGROUND

The Global Program to Eliminate Lymphatic Filariasis (GPELF) is the largest public health program based on mass drug administration (MDA). Despite decades of MDA, ongoing transmission in some countries remains a challenge. To optimise interventions, it is critical to differentiate between recrudescence and new infections. Since adult filariae are inaccessible in humans, deriving a method that relies on the offspring microfilariae (mf) is necessary.

METHODS

We developed a genome amplification and kinship analysis-based approach using Brugia malayi samples from gerbils, and applied it to analyse Wuchereria bancrofti mf from humans in Côte d'Ivoire. We examined the pre-treatment genetic diversity in 269 mf collected from 18 participants, and further analysed 1-year post-treatment samples of 74 mf from 4 participants. Hemizygosity of the male X-chromosome allowed for direct inference of haplotypes, facilitating robust maternal parentage inference. To enrich parasite DNA from samples contaminated with host DNA, a whole-exome capture panel was created for W. bancrofti.

FINDINGS

By reconstructing and temporally tracking sibling relationships across pre- and post-treatment samples, we differentiated between new and established maternal families, suggesting reinfection in one participant and recrudescence in three participants. The estimated number of reproductively active adult females ranged between 3 and 11 in the studied participants. Population structure analysis revealed genetically distinct parasites in Côte d'Ivoire compared to samples from other countries. Exome capture identified protein-coding variants with ∼95% genotype concordance rate.

INTERPRETATION

We have generated resources to facilitate the development of molecular genetic tools that can estimate adult worm burdens and monitor parasite populations, thus providing essential information for the successful implementation of GPELF.

FUNDING

This work was financially supported by the Bill and Melinda Gates Foundation (https://www.gatesfoundation.org) under grant OPP1201530 (Co-PIs PUF & Gary J. Weil). B. malayi parasite material was generated with support of the Foundation for Barnes Jewish Hospital (PUF). In addition, the development of computational methods was supported by the National Institutes of Health under grants AI144161 (MM) and AI146353 (MM). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Past volcanic activity predisposes an endemic threatened seabird to negative anthropogenic impacts

Humans are regularly cited as the main driver of current biodiversity extinction, but the impact of historic volcanic activity is often overlooked. Pre-human evidence of wildlife abundance and diversity are essential for disentangling anthropogenic impacts from natural events. Réunion Island, with its intense and well-documented volcanic activity, endemic biodiversity, long history of isolation and recent human colonization, provides an opportunity to disentangle these processes. We track past demographic changes of a critically endangered seabird, the Mascarene petrel Pseudobulweria aterrima, using genome-wide SNPs. Coalescent modeling suggested that a large ancestral population underwent a substantial population decline in two distinct phases, ca. 125,000 and 37,000 years ago, coinciding with periods of major eruptions of Piton des Neiges. Subsequently, the ancestral population was fragmented into the two known colonies, ca. 1500 years ago, following eruptions of Piton de la Fournaise. In the last century, both colonies declined significantly due to anthropogenic activities, and although the species was initially considered extinct, it was rediscovered in the 1970s. Our findings suggest that the current conservation status of wildlife on volcanic islands should be firstly assessed as a legacy of historic volcanic activity, and thereafter by the increasing anthropogenic impacts, which may ultimately drive species towards extinction.

Development of genetic markers for reproductive management of toucans

Retention of genetic diversity in successive generations is key to successful ex situ programs and will become increasingly important to restore wild populations of threatened animals. When animal genealogy is partly unknown or gaps exist in studbook records, the application of molecular resources facilitates informed breeding. Here, we apply molecular resources to an ex situ breeding population of toucans (Ramphastidae), a bird family zoos commonly maintain. Toucans face population declines from illegal poaching and habitat degradation. We developed novel microsatellite markers using blood samples from 15 Keel-billed Toucans (Ramphastos sulfuratus Lesson 1830). Parentage of two individuals was known a priori, but possible sibship among 13 putative founders-including the parents-was unknown. We compared available avian heterologous and novel microsatellite markers to recover known relationships and reconstruct sibship. Eight of 61 heterologous markers amplified consistently and were polymorphic, but less so than the 18 novel markers. Known sibship (and three sibling pairs whose relatedness was unknown a priori) and paternity-though not maternity except in one case-were well-recovered using both likelihood and pairwise relatedness methods, when incorporating novel but not heterologous markers. Zoo researchers seeking microsatellite primer sets for their breeding toucan populations will likely benefit from our heterologous markers, which can be leveraged both to assess relatedness and select breeding pairs. We recommend that zoo biologists rely on species-specific primers and not optimize heterologous primers for toucan species without molecular resources. We conclude with a brief discussion of modern genotyping methods of interest to zoo researchers.

Advances in wildlife abundance estimation using pedigree reconstruction

The conservation and management of wildlife populations, particularly for threatened and endangered species are greatly aided with abundance, growth rate, and density measures. Traditional methods of estimating abundance and related metrics represent trade-offs in effort and precision of estimates. Pedigree reconstruction is an emerging, attractive alternate approach because its use of one-time, noninvasive sampling of individuals to infer the existence of unsampled individuals. However, advances in pedigree reconstruction could improve its utility, including forming a measure of precision for the method, establishing required spatial sampling effort for accurate estimates, ascertaining the spatial extent of abundance estimates derived from pedigree reconstruction, and assessing how population density affects the estimator's performance. Using established relationships for a stochastic, spatially explicit simulated moose (Alces americanus) population, pedigree reconstruction provided accurate estimates of the adult moose population size and trend. Novel bootstrapped confidence intervals performed as expected with intensive sampling but underperformed with moderate sampling efforts that could produce abundance estimates with low bias. Adult population estimates more closely reflected the total number of adults in the extant population, rather than number of adults inhabiting the area where sampling occurred. Increasing sampling effort, measured as the proportion of individuals sampled and as the proportion of a hypothetical study area, yielded similar asymptotic patterns over time. Simulations indicated a positive relationship between animal density and sampling effort required for unbiased estimates. These results indicate that pedigree reconstruction can produce accurate abundance estimates and may be particularly valuable for surveying smaller areas and low-density populations.

Human genetic history on the Tibetan Plateau in the past 5100 years

Using genome-wide data of 89 ancient individuals dated to 5100 to 100 years before the present (B.P.) from 29 sites across the Tibetan Plateau, we found plateau-specific ancestry across plateau populations, with substantial genetic structure indicating high differentiation before 2500 B.P. Northeastern plateau populations rapidly showed admixture associated with millet farmers by 4700 B.P. in the Gonghe Basin. High genetic similarity on the southern and southwestern plateau showed population expansion along the Yarlung Tsangpo River since 3400 years ago. Central and southeastern plateau populations revealed extensive genetic admixture within the plateau historically, with substantial ancestry related to that found in southern and southwestern plateau populations. Over the past ~700 years, substantial gene flow from lowland East Asia further shaped the genetic landscape of present-day plateau populations. The high-altitude adaptive EPAS1 allele was found in plateau populations as early as in a 5100-year-old individual and showed a sharp increase over the past 2800 years.

Genetic diversity and structure in wild Robusta coffee (Coffea canephora A. Froehner) populations in Yangambi (DR Congo) and their relation to forest disturbance

Degradation and regeneration of tropical forests can strongly affect gene flow in understorey species, resulting in genetic erosion and changes in genetic structure. Yet, these processes remain poorly studied in tropical Africa. Coffea canephora is an economically important species, found in the understorey of tropical rainforests of Central and West Africa, and the genetic diversity harboured in its wild populations is vital for sustainable coffee production worldwide. Here, we aimed to quantify genetic diversity, genetic structure, and pedigree relations in wild C. canephora populations, and we investigated associations between these descriptors and forest disturbance and regeneration. Therefore, we sampled 256 C. canephora individuals within 24 plots across three forest categories in Yangambi (DR Congo), and used genotyping-by-sequencing to identify 18,894 SNPs. Overall, we found high genetic diversity, and no evidence of genetic erosion in C. canephora in disturbed old-growth forest, as compared to undisturbed old-growth forest. In addition, an overall heterozygosity excess was found in all populations, which was expected for a self-incompatible species. Genetic structure was mainly a result of isolation-by-distance, reflecting geographical location, with low to moderate relatedness at finer scales. Populations in regrowth forest had lower allelic richness than populations in old-growth forest and were characterised by a lower inter-individual relatedness and a lack of isolation-by-distance, suggesting that they originated from different neighbouring populations and were subject to founder effects. Wild Robusta coffee populations in the study area still harbour high levels of genetic diversity, yet careful monitoring of their response to ongoing forest degradation remains required.

A simulation-based evaluation of methods for estimating census population size of terrestrial game species from genetically-identified parent-offspring pairs

Estimates of wildlife population size are critical for conservation and management, but accurate estimates are difficult to obtain for many species. Several methods have recently been developed that estimate abundance using kinship relationships observed in genetic samples, particularly parent-offspring pairs. While these methods are similar to traditional Capture-Mark-Recapture, they do not need physical recapture, as individuals are considered recaptured if a sample contains one or more close relatives. This makes methods based on genetically-identified parent-offspring pairs particularly interesting for species for which releasing marked animals back into the population is not desirable or not possible (e.g., harvested fish or game species). However, while these methods have successfully been applied in commercially important fish species, in the absence of life-history data, they are making several assumptions unlikely to be met for harvested terrestrial species. They assume that a sample contains only one generation of parents and one generation of juveniles of the year, while more than two generations can coexist in the hunting bags of long-lived species, or that the sampling probability is the same for each individual, an assumption that is violated when fecundity and/or survival depend on sex or other individual traits. In order to assess the usefulness of kin-based methods to estimate population sizes of terrestrial game species, we simulated population pedigrees of two different species with contrasting demographic strategies (wild boar and red deer), applied four different methods and compared the accuracy and precision of their estimates. We also performed a sensitivity analysis, simulating population pedigrees with varying fecundity characteristics and various levels of harvesting to identify optimal conditions of applicability of each method. We showed that all these methods reached the required levels of accuracy and precision to be effective in wildlife management under simulated circumstances (i.e., for species within a given range of fecundity and for a given range of sampling intensity), while being robust to fecundity variation. Despite the potential usefulness of the methods for terrestrial game species, care is needed as several biases linked to hunting practices still need to be investigated (e.g., when hunting bags are biased toward a particular group of individuals).

A new GTSeq resource to facilitate multijurisdictional research and management of walleye Sander vitreus

Conservation and management professionals often work across jurisdictional boundaries to identify broad ecological patterns. These collaborations help to protect populations whose distributions span political borders. One common limitation to multijurisdictional collaboration is consistency in data recording and reporting. This limitation can impact genetic research, which relies on data about specific markers in an organism's genome. Incomplete overlap of markers between separate studies can prevent direct comparisons of results. Standardized marker panels can reduce the impact of this issue and provide a common starting place for new research. Genotyping-in-thousands (GTSeq) is one approach used to create standardized marker panels for nonmodel organisms. Here, we describe the development, optimization, and early assessments of a new GTSeq panel for use with walleye (Sander vitreus) from the Great Lakes region of North America. High genome-coverage sequencing conducted using RAD capture provided genotypes for thousands of single nucleotide polymorphisms (SNPs). From these markers, SNP and microhaplotype markers were chosen, which were informative for genetic stock identification (GSI) and kinship analysis. The final GTSeq panel contained 500 markers, including 197 microhaplotypes and 303 SNPs. Leave-one-out GSI simulations indicated that GSI accuracy should be greater than 80% in most jurisdictions. The false-positive rates of parent-offspring and full-sibling kinship identification were found to be low. Finally, genotypes could be consistently scored among separate sequencing runs >94% of the time. Results indicate that the GTSeq panel that we developed should perform well for multijurisdictional walleye research throughout the Great Lakes region.

Preliminary insight into parental contributions to Japanese eel (Anguilla japonica) preleptocephali spawned on different nights

Parentage sibship-inference analyses were conducted using mtDNA sequencing and six microsatellite genotypes of 182 Japanese eel preleptocephali that were collected from one net-tow near the West Mariana Ridge in May 2014. At least 328 parents were involved in producing the 182 preleptocephali, and several parents may have spawned a few times during 3 days of a spawning period. Half-sibs suggested that a few parents mated with 1-3 partners, indicating that the Japanese eel can form spawning aggregations in which several parents mate with each other in the ocean.

Parentage Assignment Using Microsatellites Reveals Multiple Mating in Aedes aegypti (Diptera: Culicidae): Implications for Mating Dynamics

The mosquito Aedes aegypti is the primary vector of the dengue, yellow fever, and chikungunya viruses. Evidence shows that Ae. aegypti males are polyandrous whereas Ae. aegypti females are monandrous in mating. However, the degree to which Ae. aegypti males and females can mate with different partners has not been rigorously tested. Therefore, this study examined the rates of polyandry via parentage assignment in three sets of competitive mating experiments using wild-type male and female Ae. aegypti. Parentage assignment was monitored using nine microsatellite DNA markers. All Ae. aegypti offspring were successfully assigned to parents with 80% or 95% confidence using CERVUS software. The results showed that both male and female Ae. aegypti mated with up to 3-4 different partners. Adults contributed differentially to the emergent offspring, with reproductive outputs ranging from 1 to 25 viable progeny. This study demonstrates a new perspective on the capabilities of male and female Ae. aegypti in mating. These findings are significant because successful deployment of reproductive control methods using genetic modification or sterile Ae. aegypti must consider the following criteria regarding their mating fitness: 1) choosing Ae. aegypti males that can mate with many different females; 2) testing how transformed Ae. aegypti male perform with polyandrous females; and 3) prioritizing the selection of polyandrous males and/or females Ae. aegypti that have the most offspring.

Pedigree reconstruction and genetic analysis of major ornamental characters of ornamental crabapple (Malus spp.) based on paternity analysis

Ornamental crabapple is an important woody ornamental plant in the Northern Hemisphere. Its flowers, fruits, leaves and tree habit are all important ornamental characters. As there has been no research on the selection of superior parents and phenotypic variation, new varieties of ornamental crabapple are mainly selected from open-pollination progeny. In order to explore the transmission rule of ornamental traits between parents and offspring of crabapple, and to provide a basis for the selection of hybrid parents for directional breeding, 14 pairs of SSR markers were used in this study for paternity analysis of 384 offspring from 4 female parents crossed with 91 candidate male parents. And 273 offspring (71.1%) were matched with only the father at a 95% strict confidence level. We reconstructed 7 full-sib families (number of progeny ≥ 10) on the basis of the paternity analysis results. Genetic analysis of characters in the full-sib families revealed that green leaves and white flowers were dominant traits. All the hybrid offspring from the white flower (♀) × non-white flower (♂) cross produced white flowers, while 7.04% produced non-white flowers when both parents had white flowers. The results showed that white flowers might be a dominant qualitative trait in crabapple, while the depth of red was a quantitative trait. The genetic characteristics of green and non-green leaves and the depth of red of the peel were similar to flower color. Compared with the upright and spreading traits, the weeping trait was recessive. Some progeny showed an earlier blooming period, indicating the possibility of breeding for blooming period. Our findings are important for parent screening and improving the breeding efficiency of new varieties in ornamental crabapple hybridization.

Evaluating sources of bias in pedigree-based estimates of breeding population size

Applications of genetic-based estimates of population size are expanding, especially for species for which traditional demographic estimation methods are intractable due to the rarity of adult encounters. Estimates of breeding population size (N_S ) are particularly amenable to genetic-based approaches as the parameter can be estimated using pedigrees reconstructed from genetic data gathered from discrete juvenile cohorts, therefore eliminating the need to sample adults in the population. However, a critical evaluation of how genotyping and sampling effort influence bias in pedigree reconstruction, and how these biases subsequently influence estimates of N_S , is needed to evaluate the efficacy of the approach under a range of scenarios. We simulated a model system to understand the interactive effects of genotyping and sampling effort on error in genetic pedigrees reconstructed from the program COLONY. We then evaluated how errors in pedigree reconstruction influenced bias and precision in estimates of N_S using three different rarefaction estimators. Results indicated that pedigree error can be minimal when adequate genetic data are available, such as when juvenile sample sizes are large and/or individuals are genotyped at many informative loci. However, even in cases for which data are limited, using results of the simulation analysis to understand the magnitude and sources of bias in reconstructed pedigrees can still be informative when estimating N_S . We applied results of the simulation analysis to evaluate N ̂ $$ \hat{N} $$ _S for a population of federally endangered Atlantic sturgeon (Acipenser oxyrinchus oxyrinchus) in the Delaware River, USA. Our results indicated that N_S is likely to be three orders of magnitude lower compared with historic breeding population sizes, which is a considerable advancement in our understanding of the population status of Atlantic sturgeon in the Delaware River. Our analyses are broadly applicable in the design and interpretation of studies seeking to estimate N_S and can help to guide conservation decisions when ecological uncertainty is high. The utility of these results is expected to grow as rapid advances in genetic technologies increase the popularity of genetic population monitoring and estimation.

A Noninvasive Genetic Insight into the Spatial and Social Organization of an Endangered Population of the Eurasian Otter (Lutra lutra, Mustelidae, Carnivora)

The Eurasian otter is endangered in Italy, only surviving in southern river basins. The spatial and social structure of a population living at the border of the current range was explored through a noninvasive genetic study along 174 km of the Sangro river. Sampling was conducted in 2011 and 2012, collecting spraints and anal jellies at 62 marking sites. Samples were successfully genotyped at 13 nuclear microsatellites and the ZFX/ZFY locus for molecular sexing, resulting in 14 distinct genotypes (4 females, 2 possible females, 8 males), from 35 marking sites. Mean captures/recaptures rate was 3.8 captures/individual, with males being recaptured more frequently than females. Spatial overlap among individuals was analyzed through a linear regression model fitted against sibship categories and sex pairing. Nine out of the fourteen genotyped individuals belonged to three full-sib clusters, while five individuals had no full-sibs in the population. Full-sibs overlapped more than half-sibs, while male–male pairs showed significantly higher spatial overlap than both male–female and female–female pairs. Estimated mean density was 0.152 otters/km and 2.4 individuals/10 × 10 km grid cell. Accordingly, the 3440 grid cells of otter occurrence in Italy could likely host about 8000 otters, suggesting the current population has become larger than the minimum viable population size.

Kinship Analysis and Pedigree Reconstruction of a Natural Regenerated Cork Oak (Quercus suber) Population

Cork oak (Quercus suber L.) is a valuable forest species in the western Mediterranean Basin due to its ecological value and the production of cork (a renewable natural material). Cork quality depends on the genetic background and cork oak environment, which has long been recognized. As no cork oak genetic trials with pedigree information were available, the inference of the genetic relatedness between individuals from molecular markers can potentially be applied to natural populations. This work aimed to investigate the potential of performing kinship prediction and pedigree reconstruction by SNP genotyping a natural cork oak population. A total of 494 trees located in Portugal were genotyped with 8K SNPs. The raw SNP set was filtered differently, producing four SNP sets that were further filtered by missing data, genotype frequency, and minor allele frequency. For each set, an identity by descent (IBD) matrix was generated to perform the relationship prediction, revealing from 22,114 to 23,859 relationships. Familial categories from the first to the third degree were able to be assigned. The feasibility of SNP genotyping for future studies on the kinship analysis and pedigree reconstruction of cork oak populations was demonstrated. The information produced may be used in further breeding and conservation programs for cork oak.

Phenotypic Diversity and Association Mapping of Ascorbic Acid Content in Spinach

Ascorbic acid (AsA), or vitamin C, is an essential nutrient for humans. In plants, AsA functions as an antioxidant during normal metabolism or in response to stress. Spinach is a highly nutritious green leafy vegetable that is consumed fresh, cooked or as a part of other dishes. One current goal in spinach breeding programs is to enhance quality and nutritional content. However, little is known about the diversity of nutritional content present in spinach germplasm, especially for AsA content. In this study, a worldwide panel of 352 accessions was screened for AsA content showing that variability in spinach germplasm is high and could be utilized for cultivar improvement. In addition, a genome-wide association study for marker-trait association was performed using three models, and associated markers were searched in the genome for functional annotation analysis. The generalized linear model (GLM), the compressed mixed linear model (CMLM) based on population parameters previously determined (P3D) and the perMarker model together identified a total of 490 significant markers distributed across all six spinach chromosomes indicating the complex inheritance of the trait. The different association models identified unique and overlapping marker sets, where 27 markers were identified by all three models. Identified high AsA content accessions can be used as parental lines for trait introgression and to create segregating populations for further genetic analysis. Bioinformatic analysis indicated that identified markers can differentiate between high and low AsA content accessions and that, upon validation, these markers should be useful for breeding programs.

Microsatellite Polymorphism and the Population Structure of Dugongs (Dugong dugon) in Thailand

Simple Summary

For this study, skin samples were analyzed from 77 individual stranded dugongs collected in Thai waters from 1994–2019 using six microsatellite markers to assess the genetic diversity and population structure. Dugongs in the Andaman Sea had higher genetic variation than those in the Gulf of Thailand. Populations in Trang, Satun and some areas of Krabi had highest diversity compared to other regions of Thailand. The analysis of Bayesian genetic clustering showed that dugongs in Thailand consist of five genetic groups. Furthermore, dugongs in the middle and lower Andaman Sea presented the greatest gene flow compared to other regions. Based on calculation of inbreeding coefficients, dugong populations in the Sea of Thailand are experiencing some levels of inbreeding, and so may warrant special protections. Results of this study provide important information on genetic diversity and genetic population structuring of dugongs in Thailand and for understanding the genetic status of dugongs that can lead to improved management and conservation of this endangered species.

Abstract

The dugong (Dugong dugon) is an endangered species of marine mammals, so knowledge of genetic diversity of these populations is important for conservation planning within different habitats. In this study, six microsatellite markers were used to assess the genetic diversity and population structure of 77 dugongs from skin samples of stranded animals collected from 1994–2019 (69 from Andaman Sea and 8 from the Gulf of Thailand). Our results found that dugongs in the Andaman Sea had higher genetic variation than those in the Gulf of Thailand. Populations in Trang, Satun, and some areas of Krabi had highest diversity compared to other regions of Thailand. Bayesian genetic clustering analysis revealed that dugongs in Thailand consist of five genetic groups. Moreover, dugongs in the middle and lower Andaman Sea presented the greatest gene flow compared to other regions. However, based on calculation of inbreeding coefficients (Fis value = 0.239), dugong populations in the Sea of Thailand are experiencing some levels of inbreeding, and so may warrant special protections. These results provide important information for understanding the genetic status of dugongs that can lead to improved management and conservation of this endangered species.

Impact of model assumptions on demographic inferences: the case study of two sympatric mouse lemurs in northwestern Madagascar

BACKGROUND

Quaternary climate fluctuations have been acknowledged as major drivers of the geographical distribution of the extraordinary biodiversity observed in tropical biomes, including Madagascar. The main existing framework for Pleistocene Malagasy diversification assumes that forest cover was strongly shaped by warmer Interglacials (leading to forest expansion) and by cooler and arid glacials (leading to forest contraction), but predictions derived from this scenario for forest-dwelling animals have rarely been tested with genomic datasets.

RESULTS

We generated genomic data and applied three complementary demographic approaches (Stairway Plot, PSMC and IICR-simulations) to infer population size and connectivity changes for two forest-dependent primate species (Microcebus murinus and M. ravelobensis) in northwestern Madagascar. The analyses suggested major demographic changes in both species that could be interpreted in two ways, depending on underlying model assumptions (i.e., panmixia or population structure). Under panmixia, the two species exhibited larger population sizes across the Last Glacial Maximum (LGM) and towards the African Humid Period (AHP). This peak was followed by a population decline in M. ravelobensis until the present, while M. murinus may have experienced a second population expansion that was followed by a sharp decline starting 3000 years ago. In contrast, simulations under population structure suggested decreasing population connectivity between the Last Interglacial and the LGM for both species, but increased connectivity during the AHP exclusively for M. murinus.

CONCLUSION

Our study shows that closely related species may differ in their responses to climatic events. Assuming that Pleistocene climatic conditions in the lowlands were similar to those in the Malagasy highlands, some demographic dynamics would be better explained by changes in population connectivity than in population size. However, changes in connectivity alone cannot be easily reconciled with a founder effect that was shown for M. murinus during its colonization of the northwestern Madagascar in the late Pleistocene. To decide between the two alternative models, more knowledge about historic forest dynamics in lowland habitats is necessary. Altogether, our study stresses that demographic inferences strongly depend on the underlying model assumptions. Final conclusions should therefore be based on a comparative evaluation of multiple approaches.

Are immigrants outbred and unrelated? Testing standard assumptions in a wild metapopulation

Immigration into small recipient populations is expected to alleviate inbreeding and increase genetic variation, and hence facilitate population persistence through genetic and/or evolutionary rescue. Such expectations depend on three standard assumptions: that immigrants are outbred, unrelated to existing natives at arrival, and unrelated to each other. These assumptions are rarely explicitly verified, including in key field systems in evolutionary ecology. Yet, they could be violated due to non-random or repeated immigration from adjacent small populations. We combined molecular genetic marker data for 150-160 microsatellite loci with comprehensive pedigree data to test the three assumptions for a song sparrow (Melospiza melodia) population that is a model system for quantifying effects of inbreeding and immigration in the wild. Immigrants were less homozygous than existing natives on average, with mean homozygosity that closely resembled outbred natives. Immigrants can therefore be considered outbred on the focal population scale. Comparisons of homozygosity of real or hypothetical offspring of immigrant-native, native-native and immigrant-immigrant pairings implied that immigrants were typically unrelated to existing natives and to each other. Indeed, immigrants' offspring would be even less homozygous than outbred individuals on the focal population scale. The three standard assumptions of population genetic and evolutionary theory were consequently largely validated. Yet, our analyses revealed some deviations that should be accounted for in future analyses of heterosis and inbreeding depression, implying that the three assumptions should be verified in other systems to probe patterns of non-random or repeated dispersal and facilitate precise and unbiased estimation of key evolutionary parameters.

Meta-analytic evidence that animals rarely avoid inbreeding

Animals are usually expected to avoid mating with relatives (kin avoidance) as incestuous mating can lead to the expression of inbreeding depression. Yet, theoretical models predict that unbiased mating with regards to kinship should be common, and that under some conditions, the inclusive fitness benefits associated with inbreeding can even lead to a preference for mating with kin. This mismatch between empirical and theoretical expectations generates uncertainty as to the prevalence of inbreeding avoidance in animals. Here, we synthesized 677 effect sizes from 139 experimental studies of mate choice for kin versus non-kin in diploid animals, representing 40 years of research, using a meta-analytical approach. Our meta-analysis revealed little support for the widely held view that animals avoid mating with kin, despite clear evidence of publication bias. Instead, unbiased mating with regards to kinship appears widespread across animals and experimental conditions. The significance of a variety of moderators was explored using meta-regressions, revealing that the degree of relatedness and prior experience with kin explained some variation in the effect sizes. Yet, we found no difference in kin avoidance between males and females, choice and no-choice experiments, mated and virgin animals or between humans and animals. Our findings highlight the need to rethink the widely held view that inbreeding avoidance is a given in experimental studies.

Identity-by-state-based haplotyping expands the application of comprehensive preimplantation genetic testing

STUDY QUESTION

Is it possible to haplotype parents using parental siblings to leverage preimplantation genetic testing (PGT) for monogenic diseases and aneuploidy (comprehensive PGT) by genome-wide haplotyping?

SUMMARY ANSWER

We imputed identity-by-state (IBS) sharing of parental siblings to phase parental genotypes.

WHAT IS KNOWN ALREADY

Genome-wide haplotyping of preimplantation embryos is being implemented as a generic approach for genetic diagnosis of inherited single-gene disorders. To enable the phasing of genotypes into haplotypes, genotyping the direct family members of the prospective parent carrying the mutation is required. Current approaches require genotypes of either (i) both or one of the parents of the affected prospective parent or (ii) an affected or an unaffected child of the couple. However, this approach cannot be used when parents or children are not attainable, prompting an investigation into alternative phasing options.

STUDY DESIGN, SIZE, DURATION

This is a retrospective validation study, which applied IBS-based phasing of parental haplotypes in 56 embryos derived from 12 PGT families. Genome-wide haplotypes and copy number profiles generated for each embryo using the new phasing approach were compared with the reference PGT method to evaluate the diagnostic concordance.

PARTICIPANTS/MATERIALS, SETTING, METHODS

This study included 12 couples with a known hereditary genetic disorder, participating in the comprehensive PGT program and with at least one parental sibling available (e.g. brother and/or sister). Genotyping data from both prospective parents and the parental sibling(s) were used to perform IBS-based phasing and to trace the disease-associated alleles. The outcome of the IBS-based PGT was compared with the results of the clinically implemented reference haplotyping-based PGT method.

MAIN RESULTS AND THE ROLE OF CHANCE

IBS-based haplotyping was performed for 12 PGT families. In accordance with the theoretical prediction of allele sharing between sibling pairs, 6 out of 12 (50%) couples or 23 out of 56 embryos could be phased using parental siblings. In families where phasing was possible, haplotype calling in the locus of interest was 100% concordant between the reference PGT method and IBS-based approach using parental siblings.

LARGE SCALE DATA

N/A.

LIMITATIONS, REASONS FOR CAUTION

Phasing of parental haplotypes will only be possible when the disease locus lies in an informative region (categorized as IBS1). Phasing prospective parents using relatives with reduced genetic relatedness as a reference (e.g. siblings) decreases the size and the occurrence of informative IBS1 regions, necessary for haplotype calling. By including more than one extended family member, the chance of obtaining IBS1 coverage in the interrogated locus can be increased. A pre-PGT work-up can define whether the carrier couple could benefit from this approach.

WIDER IMPLICATIONS OF THE FINDINGS

Phasing by relatives extends the potential of comprehensive PGT, since it allows the inclusion of couples who do not have access to the standard phasing references, such as parents or offspring.

STUDY FUNDING/COMPETING INTEREST(S)

The study was funded by the KU Leuven grant (C14/18/092), Research Foundation Flanders (FWO; GA09311N), Horizon 2020 innovation programme (WIDENLIFE, 692065) and Agilent Technologies. J.R.V., T.V. and M.Z.E. are co-inventors of a patent ZL910050-PCT/EP2011/060211-WO/2011/157846 'Methods for haplotyping single-cells' and ZL913096-PCT/EP2014/068315-WO/2015/028576 'Haplotyping and copy number typing using polymorphic variant allelic frequencies' licensed to Agilent Technologies. The other authors have no conflict of interest to declare.

Kin relationships in cultural species of the marine realm: case study of a matrilineal social group of sperm whales off Mauritius island, Indian Ocean

Understanding the organization and dynamics of social groups of marine mammals through the study of kin relationships is particularly challenging. Here, we studied a stable social group of sperm whales off Mauritius, using underwater observations, individual-specific identification, non-invasive sampling and genetic analyses based on mitochondrial sequencing and microsatellite profiling. Twenty-four sperm whales were sampled between 2017 and 2019. All individuals except one adult female shared the same mitochondrial DNA (mtDNA) haplotype-one that is rare in the western Indian Ocean-thus confirming with near certainty the matrilineality of the group. All probable first- and second-degree kin relationships were depicted in the sperm whale social group: 13 first-degree and 27 second-degree relationships were identified. Notably, we highlight the likely case of an unrelated female having been integrated into a social unit, in that she presented a distinct mtDNA haplotype and no close relationships with any members of the group. Investigating the possible matrilineality of sperm whale cultural units (i.e. vocal clans) is the next step in our research programme to elucidate and better apprehend the complex organization of sperm whale social groups.

Performing parentage analysis for polysomic inheritances based on allelic phenotypes

Polyploidy poses several problems for parentage analysis. We present a new polysomic inheritance model for parentage analysis based on genotypes or allelic phenotypes to solve these problems. The effects of five factors are simultaneously accommodated in this model: (1) double-reduction, (2) null alleles, (3) negative amplification, (4) genotyping errors and (5) self-fertilization. To solve genotyping ambiguity (unknown allele dosage), we developed a new method to establish the likelihood formulas for allelic phenotype data and to simultaneously include the effects of our five chosen factors. We then evaluated and compared the performance of our new method with three established methods by using both simulated data and empirical data from the cultivated blueberry (Vaccinium corymbosum). We also developed and compared the performance of two additional estimators to estimate the genotyping error rate and the sample rate. We make our new methods freely available in the software package polygene, at http://github.com/huangkang1987/polygene.

Dispersal patterns of a recovering brown bear (Ursus arctos) population in a human-dominated landscape

Despite increasing habitat fragmentation, large carnivore populations in parts of Europe have been recovering and expanding into human-dominated areas. Knowledge of animal dispersal patterns in such areas is important for their conservation, management, and coexistence with humans. We used genetic data based on 15 microsatellite markers from 312 individuals (98 females, 214 males) to assess kinship and dispersal patterns during the recovery and spatial expansion of a wild brown bear (Ursus arctos) population (2003–2010) in the human-dominated landscape of Greece. We hypothesized that bear dispersal in Greece was sex-biased, with females being more philopatric and males dispersing more frequently and over greater distances. Dispersal indeed was sex-biased, with males dispersing more frequently and farther than females. Overall, females were found to be philopatric; males also appeared to be philopatric, but to a lesser degree. However, a high proportion of females displayed dispersal behavior, which may be indicative of a pre-saturation stage of the population in that part of the country. Our results indicate that dispersal may be due to evading competition and avoiding inbreeding. We also documented long-distance dispersal of bears, which is considered to be indicative of a spatially expanding population. Our results highlight the value of using noninvasive genetic monitoring data to assess kinship among individuals and study dispersal patterns in human-dominated landscapes. Brown bears remain threatened in Greece; we therefore recommend systematic genetic monitoring of the species in combination with careful habitat management to protect suitable habitat (i.e., dispersal corridors) and ultimately ensure co-existence with humans and survival of brown bears in the country.

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.

Mode and Tempo of Microsatellite Evolution across 300 Million Years of Insect Evolution

Microsatellites are short, repetitive DNA sequences that can rapidly expand and contract due to slippage during DNA replication. Despite their impacts on transcription, genome structure, and disease, relatively little is known about the evolutionary dynamics of these short sequences across long evolutionary periods. To address this gap in our knowledge, we performed comparative analyses of 304 available insect genomes. We investigated the impact of sequence assembly methods and assembly quality on the inference of microsatellite content, and we explored the influence of chromosome type and number on the tempo and mode of microsatellite evolution across one of the most speciose clades on the planet. Diploid chromosome number had no impact on the rate of microsatellite evolution or the amount of microsatellite content in genomes. We found that centromere type (holocentric or monocentric) is not associated with a difference in the amount of microsatellite content; however, in those species with monocentric chromosomes, microsatellite content tends to evolve faster than in species with holocentric chromosomes.

Large-scale connectivity, cryptic population structure, and relatedness in Eastern Pacific Olive ridley sea turtles (Lepidochelys olivacea)

Endangered species are grouped into genetically discrete populations to direct conservation efforts. Mitochondrial control region (mtCR) haplotypes are used to elucidate deep divergences between populations, as compared to nuclear microsatellites that can detect recent structuring. When prior populations are unknown, it is useful to subject microsatellite data to clustering and/or ordination population inference. Olive ridley sea turtles (Lepidochelys olivacea) are the most abundant sea turtle, yet few studies have characterized olive ridley population structure. Recently, clustering results of olive ridleys in the Eastern Tropical Pacific Ocean suggested weak structuring (F ST = 0.02) between Mexico and Central America. We analyzed mtCR haplotypes, new microsatellite genotypes from Costa Rica, and preexisting microsatellite genotypes from olive ridleys across the Eastern Tropical Pacific, to further explore population structuring in this region. We subjected inferred populations to multiple analyses to explore the mechanisms behind their structuring. We found 10 mtCR haplotypes from 60 turtles nesting at three sites in Costa Rica, but did not detect divergence between Costa Rican sites, or between Central America and Mexico. In Costa Rica, clustering suggested one population with no structuring, but ordination suggested four cryptic clusters with moderate structuring (F ST = 0.08, p < .001). Across the Eastern Tropical Pacific, ordination suggested nine cryptic clusters with moderate structuring (F ST = 0.103, p < .001) that largely corresponded to Mexican and Central American populations. All ordination clusters displayed significant internal relatedness relative to global relatedness (p < .001) and contained numerous sibling pairs. This suggests that broadly dispersed family lineages have proliferated in Eastern Tropical Pacific olive ridleys and corroborates previous work showing basin-wide connectivity and shallow population structure in this region. The existence of broadly dispersed kin in Eastern Tropical Pacific olive ridleys has implications for management of olive ridleys in this region, and adds to our understanding of sea turtle ecology and life history, particularly in light of the natal-homing paradigm.

Assessing the genetic diversity in Argopecten nucleus (Bivalvia: Pectinidae), a functional hermaphrodite species with extremely low population density and self-fertilization: Effect of null alleles

Argopecten nucleus is a functional hermaphroditic pectinid species that exhibits self-fertilization, whose natural populations have usually very low densities. In the present study, the genetic diversity of a wild population from Neguanje Bay, Santa Marta (Colombia), was estimated using microsatellite markers, and the effect of the presence of null alleles on this estimation was assessed. A total of 8 microsatellite markers were developed, the first described for this species, and their amplification conditions were standardized. They were used to determine the genotype of 48 wild individuals from Naguanje Bay, and 1,010 individuals derived from the offspring of 38 directed crosses. For each locus, the frequencies of the identified alleles, including null alleles, were estimated using the statistical package Micro-Checker, and the parental genotypes were confirmed using segregation analysis. Three to 8 alleles per locus with frequencies from 0.001 to 0.632 were detected. The frequencies of null alleles ranged from 0.10 to 0.45, with Ho from 0.0 to 0.79, and He from 0.53 to 0.80. All loci were in H-W disequilibrium. The null allele frequencies values were high, with lower estimations using segregation analysis than estimated using Micro-Checker. The present results show high levels of population genetic diversity and indicate that null alleles were not the only cause of deviation from H-W equilibrium in all loci, suggesting that the wild population under study presents signs of inbreeding and Wahlund effect.

Multi-targeted management of upland game birds at the agroecosystem interface in midwestern North America

Despite its imperative, biodiversity conservation is chronically underfunded, a deficiency that often forces management agencies to prioritize. Single-species recovery thus becomes a focus (often with socio-political implications), whereas a more economical approach would be the transition to multi-targeted management (= MTM). This challenge is best represented in Midwestern North America where biodiversity has been impacted by 300+ years of chronic anthropogenic disturbance such that native tall-grass prairie is now supplanted by an agroecosystem. Here, we develop an MTM with a population genetic metric to collaboratively manage three Illinois upland gamebirds: common pheasant (Phasianus colchicus; pheasant), northern bobwhite quail (Colinus virginianus; quail), and threatened-endangered (T&E) greater prairie chicken (Tympanuchus cupido pinnatus; prairie chicken). We first genotyped our study pheasant at 19 microsatellite DNA loci and identified three captive breeding stocks (N = 143; IL Department of Natural Resources) as being significantly bottlenecked, with relatedness >1st-cousin (μR = 0.158). 'Wild' (non-stocked) pheasant [N = 543; 14 Pheasant-Habitat-Areas (PHAs)] were also bottlenecked, significantly interrelated (μR = 0.150) and differentiated (μFST = 0.047), yet distinct from propagation stock. PHAs that encompassed significantly with larger areas also reflected greater effective population sizes (μNE = 43; P<0.007). We juxtaposed these data against previously published results for prairie chicken and quail, and found population genetic structure driven by drift, habitat/climate impacts, and gender-biased selection via hunter-harvest. Each species (hunter-harvested or T&E) is independently managed, yet their composite population genetic baseline provides the quantitative criteria needed for an upland game bird MTM. Its implementation would require agricultural plots to be rehabilitated/reclaimed using a land-sharing/sparing portfolio that differs markedly from the Conservation Reserve Program (CRP), where sequestered land decreases as agricultural prices escalate. Cost-savings for an MTM would accrue by synchronizing single-species management with a dwindling hunter-harvest program, and by eliminating propagation/stocking programs. This would sustain not only native grasslands and their resident species, but also accelerate conservation at the wildlife-agroecosystem interface.

Unsustainable anthropogenic mortality disrupts natal dispersal and promotes inbreeding in leopards

Anthropogenic mortality of wildlife is typically inferred from measures of the absolute decline in population numbers. However, increasing evidence suggests that indirect demographic effects including changes to the age, sex, and social structure of populations, as well as the behavior of survivors, can profoundly impact population health and viability. Specifically, anthropogenic mortality of wildlife (especially when unsustainable) and fragmentation of the spatial distribution of individuals (home-ranges) could disrupt natal dispersal mechanisms, with long-term consequences to genetic structure, by compromising outbreeding behavior and gene flow. We investigate this threat in African leopards (Panthera pardus pardus), a polygynous felid with male-biased natal dispersal. Using a combination of spatial (home-range) and genetic (21 polymorphic microsatellites) data from 142 adult leopards, we contrast the structure of two South African populations with markedly different histories of anthropogenically linked mortality. Home-range overlap, parentage assignment, and spatio-genetic autocorrelation together show that historical exploitation of leopards in a recovering protected area has disrupted and reduced subadult male dispersal, thereby facilitating opportunistic male natal philopatry, with sons establishing territories closer to their mothers and sisters. The resultant kin-clustering in males of this historically exploited population is comparable to that of females in a well-protected reserve and has ultimately led to localized inbreeding. Our findings demonstrate novel evidence directly linking unsustainable anthropogenic mortality to inbreeding through disrupted dispersal in a large, solitary felid and expose the genetic consequences underlying this behavioral change. We therefore emphasize the importance of managing and mitigating the effects of unsustainable exploitation on local populations and increasing habitat fragmentation between contiguous protected areas by promoting in situ recovery and providing corridors of suitable habitat that maintain genetic connectivity.

Spatial networks differ when food supply changes: Foraging strategy of Egyptian fruit bats

Animals are faced with a range of ecological constraints that shape their behavioural decisions. Habitat features that affect resource abundance will also have an impact, especially as regards spatial distribution, which will in turn affect associations between the animals. Here we utilised a network approach, using spatial and genetic data, to describe patterns in use of space (foraging sites) by free-ranging Egyptian fruit bats (Rousettus aegyptiacus) at the Dakhla Oasis in Egypt. We observed a decrease in home range size during spring, when food availability was lowest, which was reflected by differences in space sharing networks. Our data showed that when food was abundant, space sharing networks were less connected and more related individuals shared more foraging sites. In comparison, when food was scarce the bats had few possibilities to decide where and with whom to forage. Overall, both networks had high mean degree, suggesting communal knowledge of predictable food distribution.

RAPTURE (RAD capture) panel facilitates analyses characterizing sea lamprey reproductive ecology and movement dynamics

Genomic tools are lacking for invasive and native populations of sea lamprey (Petromyzon marinus). Our objective was to discover single nucleotide polymorphism (SNP) loci to conduct pedigree analyses to quantify reproductive contributions of adult sea lampreys and dispersion of sibling larval sea lampreys of different ages in Great Lakes tributaries. Additional applications of data were explored using additional geographically expansive samples. We used restriction site-associated DNA sequencing (RAD-Seq) to discover genetic variation in Duffins Creek (DC), Ontario, Canada, and the St. Clair River (SCR), Michigan, USA. We subsequently developed RAD capture baits to genotype 3,446 RAD loci that contained 11,970 SNPs. Based on RAD capture assays, estimates of variance in SNP allele frequency among five Great Lakes tributary populations (mean F ST 0.008; range 0.00-0.018) were concordant with previous microsatellite-based studies; however, outlier loci were identified that contributed substantially to spatial population genetic structure. At finer scales within streams, simulations indicated that accuracy in genetic pedigree reconstruction was high when 200 or 500 independent loci were used, even in situations of high spawner abundance (e.g., 1,000 adults). Based on empirical collections of larval sea lamprey genotypes, we found that age-1 and age-2 families of full and half-siblings were widely but nonrandomly distributed within stream reaches sampled. Using the genomic scale set of SNP loci developed in this study, biologists can rapidly genotype sea lamprey in non-native and native ranges to investigate questions pertaining to population structuring and reproductive ecology at previously unattainable scales.

Genomic and phenotypic effects of inbreeding across two different hatchery management regimes in Chinook salmon

Genomic approaches permit direct estimation of inbreeding and its effect on fitness. We used genomic-based estimates of inbreeding to investigate their relationship with eight adult traits in a captive-reared Pacific salmonid that is released into the wild. Estimates were also used to determine whether alternative broodstock management approaches reduced risks of inbreeding. Specifically, 1,100 unlinked restriction-site associated (RAD) loci were used to compare pairwise relatedness, derived from a relationship matrix, and individual inbreeding, estimated by comparing observed and expected homozygosity, across four generations in two hatchery lines of Chinook salmon that were derived from the same source. The lines are managed as "integrated" with the founding wild stock, with ongoing gene flow, and as "segregated" with no gene flow. While relatedness and inbreeding increased in the first generation of both lines, possibly due to population subdivision caused by hatchery initiation, the integrated line had significantly lower levels in some subsequent generations (relatedness: F₂ -F₄ ; inbreeding F₂ ). Generally, inbreeding was similar between the lines despite large differences in effective numbers of breeders. Inbreeding did not affect fecundity, reproductive effort, return timing, fork length, weight, condition factor, and daily growth coefficient. However, it delayed spawn timing by 1.75 days per one standard deviation increase in F (~0.16). The results indicate that integrated management may reduce inbreeding but also suggest that it is relatively low in a small, segregated hatchery population that maximized number of breeders. Our findings demonstrate the utility of genomics to monitor inbreeding under alternative management strategies in captive breeding programs.

Comparisons of between-group differentiation in male kinship between bonobos and chimpanzees

Patterns of kinship among individuals in different groups have been rarely examined in animals. Two closest living relatives of humans, bonobos and chimpanzees share many characteristics of social systems including male philopatry, whereas one major difference between the two species is the nature of intergroup relationship. Intergroup relationship is basically antagonistic and males sometimes kill individuals of other groups in chimpanzees, whereas it is much more moderate in bonobos and copulations between individuals of different groups are often observed during intergroup encounters. Such behavioural differences may facilitate more frequent between-group male gene flow and greater between-group differentiation in male kinship in bonobos than in chimpanzees. Here we compared differences between average relatedness among males within groups and that among males of neighbouring groups, and between-group male genetic distance between bonobos and chimpanzees. Contrary to expectation, the differences between average relatedness among males within groups and that among males of neighbouring groups were significantly greater in bonobos than in chimpanzees. There were no significant differences in autosomal and Y-chromosomal between-group male genetic distance between the two species. Our results showed that intergroup male kinship is similarly or more differentiated in bonobos than in chimpanzees.

A population genetic study of feral cats on Christmas Island

Feral and stray cats are a major threat for endemic species on Christmas Island and have been contributing to their decline. Cats were introduced to Christmas Island in 1888 and are now distributed across the whole island. We analysed the genetic population structure and diversity of feral and stray cats on Christmas Island to evaluate connectivity across the island and the possibility of discernible populations that could be targeted separately. Results indicate no differentiated population structure across the island, with cats facing no habitat obstacles to reduce their dispersal abilities across the island. We found high kin structure, suggesting individuals breeding successfully on the whole island. With the management of domestic and feral/stray cats since 2010, removal efforts targeting the whole island have successfully reduced the effective population size of feral/stray cats in the last five years. We suggest the use of various management techniques to facilitate future removal efforts, especially in areas on the island that are difficult to access.

Minor alleles are associated with white rust (Albugo occidentalis) susceptibility in spinach (Spinacia oleracea)

Minor alleles (MA) have been associated with disease incidence in human studies, enabling the identification of diagnostic risk factors for various diseases. However, allelic mapping has rarely been performed in plant systems. The goal of this study was to determine whether a difference in MA prevalence is a strong enough risk factor to indicate a likely significant difference in disease resistance against white rust (WR; Albugo occidentalis) in spinach (Spinacia oleracea). We used WR disease severity ratings (WR-DSRs) in a diversity panel of 267 spinach accessions to define resistant- and susceptibility-associated groups within the distribution scores and then tested the single-nucleotide polymorphism (SNP) variants to interrogate the MA prevalence in the most susceptible (MS) vs. most resistant (MR) individuals using permutation-based allelic association tests. A total of 448 minor alleles associated with WR severity were identified in the comparison between the 25% MS and the 25% MR accessions, while the MA were generally similar between the two halves of the interquartile range. The minor alleles in the MS group were distributed across all six chromosomes and made up ~71% of the markers that were also strongly associated with WR in parallel performed genome-wide association study. These results indicate that susceptibility may be highly determined by the disproportionate overrepresentation of minor alleles, which could be used to select for resistant plants. Furthermore, by focusing on the distribution tails, allelic mapping could be used to identify plant markers associated with quantitative traits on the most informative segments of the phenotypic distribution.

How relatedness between mates influences reproductive success: An experimental analysis of self-fertilization and biparental inbreeding in a marine bryozoan

Kin associations increase the potential for inbreeding. The potential for inbreeding does not, however, make inbreeding inevitable. Numerous factors influence whether inbreeding preference, avoidance, or tolerance evolves, and, in hermaphrodites where both self-fertilization and biparental inbreeding are possible, it remains particularly difficult to predict how selection acts on the overall inbreeding strategy, and to distinguish the type of inbreeding when making inferences from genetic markers. Therefore, we undertook an empirical analysis on an understudied type of mating system (spermcast mating in the marine bryozoan, Bugula neritina) that provides numerous opportunities for inbreeding preference, avoidance, and tolerance. We created experimental crosses, containing three generations from two populations to estimate how parental reproductive success varies across parental relatedness, ranging from self, siblings, and nonsiblings from within the same population. We found that the production of viable selfed offspring was extremely rare (only one colony produced three selfed offspring) and biparental inbreeding more common. Paternity analysis using 16 microsatellite markers confirmed outcrossing. The production of juveniles was lower for sib mating compared with nonsib mating. We found little evidence for consistent inbreeding, in terms of nonrandom mating, in adult samples collected from three populations, using multiple population genetic inferences. Our results suggest several testable hypotheses that potentially explain the overall mating and dispersal strategy in this species, including early inbreeding depression, inbreeding avoidance through cryptic mate choice, and differential dispersal distances of sperm and larvae.

A core set of microsatellite loci for yellow-throated marten, Martes flavigula: a case of inferences of family relationships

BACKGROUND

Microsatellite markers are an ideal molecular marker for population genetic studies such as population structure, pedigree, and kinship. The yellow-throated marten (Martes flavigula) is widely distributed in coniferous and deciduous forests of eastern Asia and plays the role of an indicator and umbrella species in South Korea, given the absence of top predators such as tiger and leopard.

OBJECTIVE

The aim of our study was to establish a core set of microsatellite markers that could be used for a population genetics study on M. flavigula.

METHODS

We characterized 21 di-motif microsatellites for M. flavigula by Illumina next-generation sequencing. We evaluated them for a population genetics study against five established criteria together with 33 previously developed microsatellites. We calculated relatedness values between individual yellow-throated martens in two groups that were suspected to be siblings using the selected core set of markers to confirm applicability.

RESULTS

Twenty-three loci were determined as the core set of microsatellite markers. The probability of identity P_(ID) and probability of identity between siblings P_(ID)sib of the core set was estimated as 2^-15 and 2.2^-7, respectively. Relatedness values between individuals in the two groups of M. flavigula revealed that one of the pairs was sisters, confirming that the core set can be applied to kinship studies.

CONCLUSION

The developed microsatellite core set in this study is expected to contribute to studies on molecular ecology and population structure of M. flavigula.

Joint Estimation of Pedigrees and Effective Population Size Using Markov Chain Monte Carlo

Pedigrees provide the genealogical relationships among individuals at a fine resolution and serve an important function in many areas of genetic studies. One such use of pedigree information is in the estimation of the short-term effective population size [Formula: see text], which is of great relevance in fields such as conservation genetics. Despite the usefulness of pedigrees, however, they are often an unknown parameter and must be inferred from genetic data. In this study, we present a Bayesian method to jointly estimate pedigrees and [Formula: see text] from genetic markers using Markov Chain Monte Carlo. Our method supports analysis of a large number of markers and individuals within a single generation with the use of a composite likelihood, which significantly increases computational efficiency. We show, on simulated data, that our method is able to jointly estimate relationships up to first cousins and [Formula: see text] with high accuracy. We also apply the method on a real dataset of house sparrows to reconstruct their previously unreported pedigree.

From microsatellites to single nucleotide polymorphisms for the genetic monitoring of a critically endangered sturgeon

The use of genetic information is crucial in conservation programs for the establishment of breeding plans and for the evaluation of restocking success. Short tandem repeats (STRs) have been the most widely used molecular markers in such programs, but next-generation sequencing approaches have prompted the transition to genome-wide markers such as single nucleotide polymorphisms (SNPs). Until now, most sturgeon species have been monitored using STRs. The low diversity found in the critically endangered European sturgeon (Acipenser sturio), however, makes its future genetic monitoring challenging, and the current resolution needs to be increased. Here, we describe the discovery of a highly informative set of 79 SNPs using double-digest restriction-associated DNA (ddRAD) sequencing and its validation by genotyping using the MassARRAY system. Comparing with STRs, the SNP panel proved to be highly efficient and reproducible, allowing for more accurate parentage and kinship assignments' on 192 juveniles of known pedigree and 40 wild-born adults. We explore the effectiveness of both markers to estimated relatedness and inbreeding, using simulated and empirical datasets. Interestingly, we found significant correlations between STRs and SNPs at individual heterozygosity and inbreeding that give support to a reasonable representation of whole genome diversity for both markers. These results are useful for the conservation program of A. sturio in building a comprehensive studbook, which will optimize conservation strategies. This approach also proves suitable for other case studies in which highly discriminatory genetic markers are needed to assess parentage and kinship.