Relative accuracy of three common methods of parentage analysis in natural populations

Parentage Verification and Segregation Distortion Patterns of Microsatellite Markers in Olive Flounder (Paralichthys olivaceus) Full-Sib Families

Microsatellite markers are widely used in aquaculture for genetic analysis and breeding programs, but challenges such as segregation distortion and allelic instability can impact their effectiveness in parentage verification and inheritance studies. This study evaluated 15 microsatellite loci in seven experimental olive flounder (Paralichthys olivaceus) families bred through 1:1 full-sibling crosses, assessing their utility for accurate parentage and inheritance stability. Parentage assignments were conducted within an expanded pool of 647 candidate parents (including the actual 14 parents), encompassing both closely related and moderately distant individuals. Despite increased genetic diversity, assignments maintained a high accuracy rate (99.6%), demonstrating marker robustness. Differences in delta values highlighted the influence of genetic backgrounds on assignment clarity, with some groups showing reduced distinctiveness in the expanded pool. Segregation distortion was observed at certain loci, deviating from Mendelian inheritance, likely due to meiotic drive and post-zygotic selection. These findings underscore the need for empirical validation of microsatellite loci for stable inheritance and reliable parentage in diverse breeding programs, especially with genetically similar spawners.

Genomic-Guided Conservation Actions to Restore the Most Endangered Conifer in the Mediterranean Basin

Species with extremely small population sizes are critically endangered because of reduced genetic diversity, increased inbreeding and hybridisation threats. Genomic tools significantly advance conservation by revealing genetic insights into endangered species, notably in monitoring frameworks. Sicilian fir (Abies nebrodensis) is the most endangered conifer in Europe with only 30 adult trees in an 84-ha area. Using 20,824 SNPs from RAD-seq, employing genome assembly and a custom 120 SNP-array, we evaluated genetic diversity, mating patterns, and effective population size in adult trees, 118 natural seedlings, and 2064 nursery seedlings from past conservation actions. We assessed introgression from neighbouring non-native fir plantations (~6%) and established an intra-population assisted gene flow (AGF) program selecting the most genetically dissimilar individuals and investigating the outcome through simulations. Genomic analysis unveiled significant genetic diversity among adult Sicilian firs, comparable to non-endangered Mediterranean firs with larger populations. However, the genetic diversity of the forthcoming generation declined due to high self-fertilisation, leading to marked inbreeding (FIS = 0.38) and an alarmingly low effective population size (Ne = 6). Nursery seedling monitoring revealed similar selfing rates and introgression (~2%) from non-native firs. Although intra-population AGF could help to mitigate genetic loss, it may not alleviate the species vulnerability to imminent environmental challenges, perpetuating the risk of an extinction vortex. Hence, investigating the impact of Sicilian fir population decline and selfing on inbreeding depression, along with exploring the potential of hybrids for genetic load alleviation and future adaptation, is crucial for effective conservation strategies.

Genetic diversities in wild and cultivated populations of the two closely-related medical plants species, Tripterygium Wilfordii and T. Hypoglaucum (Celastraceae)

BACKGROUND

The sustainable supply of medicinal plants is important, and cultivating and domesticating them has been suggested as an optimal strategy. However, this can lead to a loss of genetic diversity. Tripterygium wilfordii Hook. f. is a medicinal plant commonly used in traditional Chinese medicine, but its wild populations are dwindling due to excessive harvesting. To protect the species and meet the increasing demand, it is urgent to cultivate it on a large scale. However, distinguishing between T. wilfordii and T. hypoglaucum, two similar species with different medicinal properties, is challenging. Therefore, it is crucial to understand the genetic diversity and population structure of these species for their sustainable utilization.

RESULTS

In this study, we investigated the genetic diversity and population structure of the two traditional medicinal semiwoody vines plant species, Tripterygium wilfordii and T. hypoglaucum, including wild and cultivated populations using chloroplast DNA (cpDNA) sequences and microsatellite loci. Our results indicated that the two species maintain a high level of genetic divergence, indicating possible genetic bases for the different contents of bioactive compounds of the two species. T. wilfordii showed lower genetic diversity and less subdivided population structures of both markers than T. hypoglaucum. The potential factors in shaping these interesting differences might be differentiated pollen-to-seed migration rates, interbreeding, and history of population divergence. Analyses of cpDNA and microsatellite loci supported that the two species are genetically distinct entities. In addition, a significant reduction of genetic diversity was observed for cultivated populations of the two species, which mainly resulted from the small initial population size and propagated vegetative practice during their cultivation.

CONCLUSION

Our findings indicate significant genetic divergence between T. wilfordii and T. hypoglaucum. The genetic diversity and population structure analyses provide important insights into the sustainable cultivation and utilization of these medicinal plants. Accurate identification and conservation efforts are necessary for both species to ensure the safety and effectiveness of crude drug use. Our study also highlighted the importance of combined analyses of different DNA markers in addressing population genetics of medicinal plants because of the contrasts of inheritance and rates of gene flow. Large-scale cultivation programs should consider preserving genetic diversity to enhance the long-term sustainability of T. wilfordii and T. hypoglaucum. Our study proposed that some populations showed higher genetic diversity and distinctness, which can be considered with priority for conservation and as the sources for future breeding and genetic improvement.

Relatedness dynamics and sex-biased dispersal in a seasonal cycle of corn mice from intensively managed agroecosystems

The sex-biased dispersal and kinship dynamics are important factors shaping the spatial distribution of individuals and are key parameters affecting a variety of ecological and evolutionary processes. Here, we studied the spatial distribution of related individuals within a population of corn mice Calomys musculinus in a seasonal cycle to infer dispersal patterns. The sampling was carried out from spring 2005 to winter 2006 in field borders of intensively managed agroecosystems. Genotyping data from 346 individuals with 9 microsatellites showed spatial genetic structure was weak for males, but not for females. The results indicate a complex spatial kinship dynamic of related females across all seasons. Which, contrary to our expectations, dispersal distances decrease with the increase of the population abundance. Meanwhile, male dispersal distances were greater when population abundance increased and thus the availability of active females. Males disperse greater distances to mate and sire offspring with distant females as a possible inbreeding avoidance mechanism. This study shows that C. musculinus is capable of much greater scattering distances than previously reported and that dispersal occurs fluidly and without barriers across the agroecosystem. The indirect benefit of dispersal on individual fitness could be related to relaxing the competition in the natal area and increasing the mating rate. Our study highlights the value of combining genetic relatedness, fieldwork observations, and behavioral data to estimate dispersal at a fine geographical scale.

Development of eighteen novel microsatellite loci for Masked Bobwhite, Colinus virginianus ridgwayi

BACKGROUND

Masked Bobwhite (Colinus virginianus ridgwayi) is a critically-endangered New World quail species endemic to Sonoran Desert grasslands of North America. It suffered severe population declines during the nineteenth and twentieth centuries, with its persistence now reliant upon a captive breeding program that requires careful genetic management to maintain extant genetic diversity. Although nuclear microsatellite DNA markers existed for the closely related Northern Bobwhite (C. virginianus), none were available for Masked Bobwhite to inform necessary management decisions.

METHODS AND RESULTS

Paired-end Illumina© sequencing was conducted to screen the Masked Bobwhite genome for microsatellite loci. We identified 18 loci exhibiting high polymorphism and limited deviations from genetic equilibrium expectations. These loci were amplified in 78 individuals. Familial relationships were reconstructed via sibship methods and compared to manually-curated pedigree data. Thirteen of fifteen full-sibling groups in the pedigree were exactly reconstructed (86.6%). Three other full-sibling groups partially matched pedigree relationships with high statistical confidence, and likely represented pedigree inaccuracies. Four additional full-sibling pairs were identified with low statistical confidence and likely resulted from analytical artifacts.

CONCLUSIONS

The novel microsatellite loci accurately reconstructed parent-offspring and sibling relationships. These loci will be useful for guiding genetic management decisions and identifying pedigree inaccuracies in the captive breeding program.

Propensity for selfing varies within a population of hermaphroditic snails: coexistence of selfers, outcrossers and mixed-mating individuals

To understand mating-system evolution in self-compatible hermaphrodites, variation in selfing rates is highly relevant. Empirical studies are rarely designed to capture variation between individuals, instead often comparing species and populations. Yet, evolution primarily occurs within populations, rendering among-individual variation essential. Observed individual selfing rates depend on the environment (e.g. differences in mate availability) and individuals' propensity for selfing. We quantified individual variation in selfing propensity in the snail Radix balthica by conducting laboratory mating trials that manipulated mate availability (low versus moderate) and estimating selfing rates from progeny arrays. We also measured female lifetime fitness. We found substantial among-individual variation in selfing propensity, including pure selfers (32%), pure outcrossers (31%) and mixed-mating individuals that selfed and outcrossed (37%). Experimental levels of mate availability did not significantly affect selfing rates. Selfers had reduced female liftetime fitness. Our results show that the propensity for selfing can differ considerably among individuals, with similar proportions of selfers, outcrossers and mixed maters. As mate availability did not affect selfing, our 'moderate' experimental level of mate availability might still have been too low to prompt selfers to outcross. This and the observed fitness differences also cautiously suggest that investigating the heritability of selfing propensities might be worthwhile in this population.

Individual-based analyses reveal effects of behavioral and demographic variables associated with multi-annual reproductive success of male and female lake sturgeon

Quantifying effects of individual attributes and population demographic characteristics that affect inter- and intrasexual interactions and adult reproductive success, and the spatial and temporal contexts in which they are expressed is important to effective species management. Multi-year individual-based analyses using genetically determined parentage allowed the examination of variables associated with the reproductive success of male and female lake sturgeon (Acipenser fulvescens) in the well-studied population in Black Lake, Michigan, USA. Spawning lake sturgeon (a total of 599 individuals where many were captured more than once based on 1024 total captures) and larvae (N = 3436) were genotyped during each of seven consecutive years (2001-2007). Factors associated with individual reproductive success differed between sexes and varied among spawning groups within a year and among years depending on spawning date (higher reproductive success earlier in the season for females) and spawning locations (higher reproductive success in upstream spawning zones for females). Female reproductive success increased nonlinearly with increasing body size. Male reproductive success increased with increasing residence time in spawning areas and, to a modest degree, with increasing body size in a nonlinear fashion. Fixed effects of repeatability in spawn timing and location across years led to consistently higher or lower reproductive success for females. Results identified factors, including time spent at spawning areas by males and intersexual encounters and mate number, that contributed to higher interindividual variance in reproductive success and affected population levels of recruitment, the degree of subpopulation genetic structure (lack of isolation by time), and effective population size.

Characterization of 25 new microsatellite markers for the green turtle (Chelonia mydas) and cross-species amplification in other marine turtle species

BACKGROUND

The green sea turtle, Chelonia mydas, is a migratory species with a strong natal homing behavior leading to a complex population structure worldwide. The species has suffered severe declines in local populations; it is therefore crucial to understand its population dynamics and genetic structure to adopt appropriate management policies. Here, we describe the development of 25 new microsatellite markers specific to C. mydas and suitable for these analyses.

METHODS AND RESULTS

They were tested on 107 specimens from French Polynesia. An average allelic diversity of 8 alleles per locus was reported and observed heterozygosity ranged from 0.187 to 0.860. Ten loci were significantly deviant from the Hardy-Weinberg equilibrium, and 16 loci showed a moderate to high level of linkage disequilibrium (4-22%). The overall F_is was positive (0.034, p-value < 0.001), and sibship analysis revealed 12 half- or full-sibling dyads, suggesting possible inbreeding in this population. Cross-amplification tests were performed on two other marine turtle species, Caretta caretta and Eretmochelys imbricata. All loci successfully amplified on these two species, though 1 to 5 loci were monomorphic.

CONCLUSION

These new markers will not only be relevant for further analyses on the population structure of the green turtle and the two other species, but they will also be invaluable for parentage studies, for which a high number of polymorphic loci are necessary. This can provide important insight into male reproductive behavior and migration, an aspect of sea turtle biology that is of critical importance for the conservation of the species.

Pollinator foraging tactics have divergent consequences for the mating system of a tropical plant

Resolving the consequences of pollinator foraging behaviour for plant mating systems is a fundamental challenge in evolutionary ecology. Pollinators may adopt particular foraging tactics: complete trapline foraging (repeated movements along a fixed route), sample-and-shift trapline foraging (a variable route that incorporates information from previous experiences) and territorial foraging (stochastic movements within a restricted area). Studies that integrate these pollinator foraging tactics with plant mating systems are generally lacking. We investigate the consequences of particular pollinator foraging tactics for Heliconia tortuosa. We combine parentage and sibship inference analysis with simulation modelling to: estimate mating system parameters; infer the foraging tactic adopted by the pollinators; and quantify the impact of pollinator foraging tactics on mating system parameters. We found high outcrossing rates, ubiquitous multiple paternity and a pronounced departure from near-neighbour mating. We also found that plants repeatedly receive pollen from a series of particular donors. We infer that the pollinators primarily adopt complete trapline foraging and occasionally engage in sample-and-shift trapline foraging. This enhances multiple paternity without a substantial increase in near-neighbour mating. The particular pollinator foraging tactics have divergent consequences for multiple paternity and near-neighbour mating. Thus, pollinator foraging behaviour is an important driver of the ecology and evolution of plant mating systems.

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.

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.

Reduced relative fitness in hatchery-origin Pink Salmon in two streams in Prince William Sound, Alaska

Previous studies generally report that hatchery-origin Pacific Salmon (Oncorhynchus spp.) have lower relative reproductive success (RRS) than their natural-origin counterparts. We estimated the RRS of Pink Salmon (O. gorbuscha) in Prince William Sound (PWS), Alaska using incomplete pedigrees. In contrast to other RRS studies, Pink Salmon have a short freshwater life history, freshwater habitats in PWS are largely unaltered by development, and sampling was conducted without the aid of dams or weirs resulting in incomplete sampling of spawning individuals. Pink Salmon released from large-scale hatchery programs in PWS have interacted with wild populations for more than 15 generations. Hatchery populations were established from PWS populations but have subsequently been managed as separate broodstocks. Gene flow is primarily directional, from hatchery strays to wild populations. We used genetic-based parentage analysis to estimate the RRS of a single generation of stray hatchery-origin Pink Salmon in two streams, and across the odd- and even-year lineages. Despite incomplete sampling, we assigned 1745 offspring to at least one parent. Reproductive success (RS), measured as sampled adult offspring that returned to their natal stream, was significantly lower for hatchery- vs. natural-origin parents in both lineages, with RRS ranging from 0.03 to 0.47 for females and 0.05 to 0.86 for males. Generalized linear modeling for the even-year lineage indicated that RRS was lower for hatchery-origin fish, ranging from 0.42 to 0.60, after accounting for sample date (run timing), sample location within the stream, and fish length. Our results strongly suggest that hatchery-origin strays have lower fitness in the wild. The consequences of reduced RRS on wild productivity depend on whether the mechanisms underlying reduced RRS are environmentally driven, and likely ephemeral, or genetically driven, and likely persistent across generations.

An evaluation of the potential factors affecting lifetime reproductive success in salmonids

Lifetime reproductive success (LRS), the number of offspring produced over an organism's lifetime, is a fundamental component of Darwinian fitness. For taxa such as salmonids with multiple species of conservation concern, understanding the factors affecting LRS is critical for the development and implementation of successful conservation management practices. Here, we reviewed the published literature to synthesize factors affecting LRS in salmonids including significant effects of hatchery rearing, life history, and phenotypic variation, and behavioral and spawning interactions. Additionally, we found that LRS is affected by competitive behavior on the spawning grounds, genetic compatibility, local adaptation, and hybridization. Our review of existing literature revealed limitations of LRS studies, and we emphasize the following areas that warrant further attention in future research: (1) expanding the range of studies assessing LRS across different life-history strategies, specifically accounting for distinct reproductive and migratory phenotypes; (2) broadening the variety of species represented in salmonid fitness studies; (3) constructing multigenerational pedigrees to track long-term fitness effects; (4) conducting LRS studies that investigate the effects of aquatic stressors, such as anthropogenic effects, pathogens, environmental factors in both freshwater and marine environments, and assessing overall body condition, and (5) utilizing appropriate statistical approaches to determine the factors that explain the greatest variation in fitness and providing information regarding biological significance, power limitations, and potential sources of error in salmonid parentage studies. Overall, this review emphasizes that studies of LRS have profoundly advanced scientific understanding of salmonid fitness, but substantial challenges need to be overcome to assist with long-term recovery of these keystone species in aquatic ecosystems.

Pedigree simulations reveal that maternity assignment is reliable in populations with conspecific brood parasitism, incomplete parental sampling and kin structure

Modern genetic parentage methods reveal that alternative reproductive strategies are common in both males and females. Under ideal conditions, genetic methods accurately connect the parents to offspring produced by extra-pair matings or conspecific brood parasitism. However, some breeding systems and sampling scenarios present significant complications for accurate parentage assignment. We used simulated genetic pedigrees to assess the reliability of parentage assignment for a series of challenging sampling regimes that reflect realistic conditions for many brood-parasitic birds: absence of genetic samples from sires, absence of samples from brood parasites and female kin-structured populations. Using 18 microsatellite markers and empirical allele frequencies from two populations of a conspecific brood parasite, the wood duck (Aix sponsa), we simulated brood parasitism and determined maternity using two widely used programs, cervus and colony. Errors in assignment were generally modest for most sampling scenarios but differed by program: cervus suffered from false assignment of parasitic offspring, whereas colony sometimes failed to assign offspring to their known mothers. Notably, colony was able to accurately infer unsampled parents. Reducing the number of markers (nine loci rather than 18) caused the assignment error to slightly worsen with colony but balloon with cervus. One potential error with important biological implications was rare in all cases-few nesting females were incorrectly excluded as the mother of their own offspring, an error that could falsely indicate brood parasitism. We consider the implications of our findings for both a retrospective assessment of previous studies and suggestions for best practices for future studies.

Genetic connectivity between Atlantic bluefin tuna larvae spawned in the Gulf of Mexico and in the Mediterranean Sea

The highly migratory Atlantic bluefin tuna (ABFT) is currently managed as two distinct stocks, in accordance with natal homing behavior and population structuring despite the absence of barriers to gene flow. Larval fish are valuable biological material for tuna molecular ecology. However, they have hardly been used to decipher the ABFT population structure, although providing the genetic signal from successful breeders. For the first time, cooperative field collection of tuna larvae during 2014 in the main spawning area for each stock, the Gulf of Mexico (GOM) and the Mediterranean Sea (MED), enabled us to assess the ABFT genetic structure in a precise temporal and spatial frame exclusively through larvae. Partitioning of genetic diversity at nuclear microsatellite loci and in the mitochondrial control region in larvae spawned contemporarily resulted in low significant fixation indices supporting connectivity between spawners in the main reproduction area for each population. No structuring was detected within the GOM after segregating nuclear diversity in larvae spawned in two hydrographically distinct regions, the eastern GOM (eGOM) and the western GOM (wGOM), with the larvae from eGOM being more similar to those collected in the MED than the larvae from wGOM. We performed clustering of genetically characterized ABFT larvae through Bayesian analysis and by Discriminant Analysis of Principal Components (DAPC) supporting the existence of favorable areas for mixing of ABFT spawners from Western and Eastern stocks, leading to gene flow and apparent connectivity between weakly structured populations. Our findings suggest that the eastern GOM is more prone for the mixing of breeders from the two ABFT populations. Conservation of this valuable resource exploited for centuries calls for intensification of tuna ichthyoplankton research and standardization of genetic tools for monitoring population dynamics.

Minimum size limits and the reproductive value of numerous, young, mature female fish

Fisheries management relies on various catch and effort controls to preserve spawning stock biomass and maximize sustainable yields while limiting fishery impacts on marine ecosystems. These include species-specific minimum or maximum size limits to protect either small non-reproductive subadults, a portion of reproductively mature adults, or large highly fecund individuals. Protecting size classes of mature fish is expected to yield a viable source of larvae for replenishing populations and reduce the risk of recruitment overfishing, yet size-specific recruitment contributions have rarely been assessed empirically. Here, we apply genetic parentage analysis to measure the reproductive success of a size-structured population of a commercially important species of coral grouper (Plectropomus maculatus-Serranidae) in no-take marine reserves (NTMRs) in the Great Barrier Reef Marine Park, Australia. Although the per capita reproductive success of individual fish increases rapidly with body length, the numerous young mature female fish, below the minimum size limit (MSL) (38 cm total length), were responsible for generating disproportionately large contributions (36%) towards larval replenishment of both fished and reserve reefs. Our findings indicate that MSLs are an effective harvest control measure to safeguard a portion of the spawning stock biomass for coral grouper and supplement recruitment subsidies assured from NTMRs.

The smell of success: Reproductive success related to rub behavior in brown bears

Several species of bears are known to rub deliberately against trees and other objects, but little is known about why bears rub. Patterns in rubbing behavior of male and female brown bears (Ursus arctos) suggest that scent marking via rubbing functions to communicate among potential mates or competitors. Using DNA from bear hairs collected from rub objects in southwestern Alberta from 2011–2014 and existing DNA datasets from Montana and southeastern British Columbia, we determined sex and individual identity of each bear detected. Using these data, we completed a parentage analysis. From the parentage analysis and detection data, we determined the number of offspring, mates, unique rub objects where an individual was detected, and sampling occasions during which an individual was detected for each brown bear identified through our sampling methods. Using a Poisson regression, we found a positive relationship between bear rubbing behavior and reproductive success; both male and female bears with a greater number of mates and a greater number of offspring were detected at more rub objects and during more occasions. Our results suggest a fitness component to bear rubbing, indicate that rubbing is adaptive, and provide insight into a poorly understood behaviour.

Great Genetic Diversity but High Selfing Rates and Short-Distance Gene Flow Characterize Populations of a Tree (Foetidia; Lecythidaceae) in the Fragmented Tropical Dry Forest of the Mascarene Islands

Following the global trend of deforestation and degradation, tropical dry forests in the Mascarenes archipelago on Reunion has undergone harsh reduction and fragmentation within 3 centuries of human occupation. We investigated the genetic diversity, mating system, and gene flow in fragmented populations of the native tree Foetidia mauritiana (Lecythidaceae) on Reunion, using microsatellite genotyping of adults (in- and ex situ) and seed progenies (in situ only). To test genetic isolation between the Mascarene islands, we also genotyped conspecific adults on Mauritius, and trees of Foetidia rodriguesiana on Rodrigues. We found a high genetic diversity among the trees on Reunion, but no population structure (G'ST: 0.039-0.090), and an increase of the fixation index (FIS) from adults to progenies. A subsequent analysis of mating systems from progeny arrays revealed selfing rates >50% in fragmented populations and close to 100% in lone trees. A paternity analysis revealed pollen flow ranging from 15.6 to 296.1 m within fragments. At broader scale, the populations of F. mauritiana on Reunion and Mauritius are genetically differentiated. The morphologically allied taxa F. rodriguesiana and F. mauritiana are clearly isolated. Therefore, this case study shows that genetic diversity may persist after deforestation, especially in long-lived tree species, but the reproductive features may be deeply altered during this process. This would explain the low seed production and the absence of recruitment in F. mauritiana. Restoration programs should take into account these features, as well as the importance that trees ex situ represent in restoring and conserving diversity.

Fast sequence-based microsatellite genotyping development workflow

Application of high-throughput sequencing technologies to microsatellite genotyping (SSRseq) has been shown to remove many of the limitations of electrophoresis-based methods and to refine inference of population genetic diversity and structure. We present here a streamlined SSRseq development workflow that includes microsatellite development, multiplexed marker amplification and sequencing, and automated bioinformatics data analysis. We illustrate its application to five groups of species across phyla (fungi, plant, insect and fish) with different levels of genomic resource availability. We found that relying on previously developed microsatellite assay is not optimal and leads to a resulting low number of reliable locus being genotyped. In contrast, de novo ad hoc primer designs gives highly multiplexed microsatellite assays that can be sequenced to produce high quality genotypes for 20-40 loci. We highlight critical upfront development factors to consider for effective SSRseq setup in a wide range of situations. Sequence analysis accounting for all linked polymorphisms along the sequence quickly generates a powerful multi-allelic haplotype-based genotypic dataset, calling to new theoretical and analytical frameworks to extract more information from multi-nucleotide polymorphism marker systems.

Mixed mating system and variable mating patterns in tropical woody bamboos

BACKGROUND

So far, little is known in detail about mating systems of woody bamboos. Paternity analysis of offspring improved our understanding of these systems, and contributed to their germplasm conservation and genetic improvement.

RESULTS

In this study, a paternity analysis of offspring from two consecutive mass or sporadically flowering events of Dendrocalamus membranaceus and D. sinicus were conducted to determine their mating system and pollen dispersal using the program COLONY based on simple sequence repeat (SSR) markers. Two sporadically flowering populations of D. sinicus (C1, C2) obtained relatively high paternity assignments rates (69.0-71.4%). Meanwhile, among three populations of D. membranaceus, the sporadically flowering population A also had much higher paternity assignments rates (56.4%) than mass flowering populations B1(28.6%) and B2 (42.5%). Both D. membranaceus and D. sinicus had mixed mating systems while their mating patterns were variable depending on pollination conditions. The maximum pollen dispersal distances were 90 m and 4378 m for D. membranaceus and D. sinicus populations, respectively, and the mating distances of these two species focused on ranges of ca. 0-50 m and 0-1500 m, respectively.

CONCLUSIONS

These results revealed for the first time variable mating patterns in woody bamboos. This suggests half-sib seeds from the same bamboo clump may have different male parents and it is crucial to clarify genetic origin in woody bamboos' breeding programs. The results also indicate the importance of pollinators in the mating systems of tropical woody bamboos.

Greater opportunities for sexual selection in male than in female obligate brood parasitic birds

Females are expected to have evolved to be more discriminatory in mate choice than males as a result of greater reproductive investment into larger gametes (eggs vs. sperm). In turn, males are predicted to be more promiscuous than females, showing both a larger variance in the number of mates and a greater increase in reproductive success with more mates, yielding more intense sexual selection on males vs. females (Bateman's Paradigm). However, sex differences in costly parental care strategies can either reinforce or counteract the initial asymmetry in reproductive investment, which may be one cause for some studies failing to conform with predictions of Bateman's Paradigm. For example, in many bird species with small female-biased initial investment but extensive biparental care, both sexes should be subject to similar strengths of sexual selection because males and females are similarly restricted in their ability to pursue additional mates. Unlike 99% of avian species, however, obligate brood parasitic birds lack any parental care in either sex, predicting a conformation to Bateman's Paradigm. Here we use microsatellite genotyping to demonstrate that in brood parasitic brown-headed cowbirds (Molothrus ater), per capita annual reproductive success increases with the number of mates in males, but not in females. Furthermore, also as predicted, the variance of the number of mates and offspring is greater in males than in females. Thus, contrary to previous findings in this species, our results conform to predictions of the Bateman's Paradigm for taxa without parental care.

Extra-pair mating in a socially monogamous and paternal mouth-brooding cardinalfish

Many vertebrates form monogamous pairs to mate and care for their offspring. However, genetic tools have increasingly shown that offspring often arise from matings outside of the monogamous pair bond. Social monogamy is relatively common in coral reef fishes, but there have been few studies that have confirmed monogamy or extra-pair reproduction, either for males or for females. Here, long-term observations and genetic tools were applied to examine the parentage of embryos in a paternally mouth-brooding cardinalfish, Sphaeramia nematoptera. Paternal care in fishes, such as mouth-brooding, is thought to be associated with a high degree of confidence in paternity. Two years of observations confirmed that S. nematoptera form long-term pair bonds within larger groups. However, genetic parentage revealed extra-pair mating by both sexes. Of 105 broods analysed from 64 males, 30.1% were mothered by a female that was not the partner and 11.5% of broods included eggs from two females. Despite the high paternal investment associated with mouth-brooding, 7.6% of broods were fertilized by two males. Extra-pair matings appeared to be opportunistic encounters with individuals from outside the immediate group. We argue that while pair formation contributes to group cohesion, both males and females can maximize lifetime reproductive success by taking advantage of extra-pair mating opportunities.

Tracing individuals and populations of the tambaqui, Colossoma macropomum (Cuvier, 1818), from Brazilian hatcheries using microsatellite markers

BACKGROUND

In recent years, tracing of alimentary produce of animal origin has become increasingly important, for economic, food safety and ecological reasons. The tambaqui, Colossoma macropomum, is the native fish most farmed in Brazil. The reliable identification of the origin of tambaquis (wild or farmed) offered for sale to the general public has become necessary to satisfy regulatory norms and uphold consumer confidence. Molecular methods based on the analysis of DNA sequences have often been used to evaluate the potential for tracing farmed fish, given their reliability and precision.

RESULTS

Full likelihood and Bayesian approaches proved to be the most efficient for the identification, respectively, of individuals and populations for most of the fish sampled from seven hatcheries and one wild stock. The exclusion method and genetic distances were the least effective approaches for the identification of individuals and populations. The Bayesian method identified correctly more than 99% of the fry from most stocks, except those of the Santarém hatchery and River Amazon wild stock, which presented the best results for individual identification.

CONCLUSIONS

The identification of populations was effective for most hatcheries, although the identification of individuals from most stocks was hampered by the reduced genetic variability. © 2018 Society of Chemical Industry.

'apparent': a simple and flexible R package for accurate SNP-based parentage analysis in the absence of guiding information

Background

The accurate determination of parent-progeny relationships within both in situ natural populations and ex situ genetic resource collections can greatly enhance plant breeding/domestication efforts and support plant genetic resource conservation strategies. Although a range of parentage analysis tools are available, none are designed to infer such relationships using genome-wide single nucleotide polymorphism (SNP) data in the complete absence of guiding information, such as generational groups, partial pedigrees, or genders. The R package (‘apparent’) developed and presented here addresses this gap.

Results

‘apparent’ adopts a novel strategy of parentage analysis based on a test of genetic identity between a theoretically expected progeny (EP_ij), whose genotypic state can be inferred at all homozygous loci for a pair of putative parents (i and j), and all potential offspring (PO_k), represented by the k individuals of a given germplasm collection. Using the Gower Dissimilarity metric (GD), genetic identity between EP_ij and PO_k is taken as evidence that individuals i and j are the true parents of offspring k. Significance of a given triad (parental pair_ij + offspring_k) is evaluated relative to the distribution of all GD_ij|k values for the population. With no guiding information provided, ‘apparent’ correctly identified the parental pairs of 15 lines of known pedigree within a test population of 77 accessions of Actinidia arguta, a performance unmatched by five other commonly used parentage analysis tools. In the case of an inconclusive triad analysis due to the absence of one parent from the test population, ‘apparent’ can perform a subsequent dyad analysis to identify a likely single parent for a given offspring. Average dyad analysis accuracy was 73.3% in the complete absence of pedigree information but increased to 100% when minimal generational information (adults vs. progeny) was provided.

Conclusions

The ‘apparent’ R package is a fast and accurate parentage analysis tool that uses genome-wide SNP data to identify parent-progeny relationships within populations for which no a priori knowledge of family structure exists.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-2662-3) contains supplementary material, which is available to authorized users.

A simulation framework for evaluating multi-stage sampling designs in populations with spatially structured traits

Selecting an appropriate and efficient sampling strategy in biological surveys is a major concern in ecological research, particularly when the population abundance and individual traits of the sampled population are highly structured over space. Multi-stage sampling designs typically present sampling sites as primary units. However, to collect trait data, such as age or maturity, only a sub-sample of individuals collected in the sampling site is retained. Therefore, not only the sampling design, but also the sub-sampling strategy can have a major impact on important population estimates, commonly used as reference points for management and conservation. We developed a simulation framework to evaluate sub-sampling strategies from multi-stage biological surveys. Specifically, we compare quantitatively precision and bias of the population estimates obtained using two common but contrasting sub-sampling strategies: the random and the stratified designs. The sub-sampling strategy evaluation was applied to age data collection of a virtual fish population that has the same statistical and biological characteristics of the Eastern Bering Sea population of Pacific cod. The simulation scheme allowed us to incorporate contributions of several sources of error and to analyze the sensitivity of the different strategies in the population estimates. We found that, on average across all scenarios tested, the main differences between sub-sampling designs arise from the inability of the stratified design to reproduce spatial patterns of the individual traits. However, differences between the sub-sampling strategies in other population estimates may be small, particularly when large sub-sample sizes are used. On isolated scenarios (representative of specific environmental or demographic conditions), the random sub-sampling provided better precision in all population estimates analyzed. The sensitivity analysis revealed the important contribution of spatial autocorrelation in the error of population trait estimates, regardless of the sub-sampling design. This framework will be a useful tool for monitoring and assessment of natural populations with spatially structured traits in multi-stage sampling designs.

Variance in Reproductive Success is Driven by Environmental Factors, not Mating System, in Bonytails

Studying the reproductive ecology of aggregate broadcast spawning fishes is difficult because it generally is not feasible to sample all potential parents and unambiguously assign their offspring. We used molecular-based parentage analysis to gain insights into the reproductive ecology of the endangered Bonytail, and to evaluate whether protected off-channel habitats could be used as an alternative to hatchery production. By genotyping adults and offspring stocked (n = 4130) into two experimental backwaters across three years, we determined that most adults (82-97%) contributed to progeny production across years and backwaters, with one exception. Both sexes mated multiply and the number of mates and family size were positively correlated. There was also a positive correlation between adult size and metrics of reproductive success. There were strong interactions between sample years and backwaters suggesting that environmental factors are the primary driver of variance in reproductive success. Knowledge of mating systems and sources of variance in reproductive success is important for management of endangered fish because high variance in reproductive success leads to substantial losses of genetic variation when few individuals reproduce successfully. Although variance in reproductive success was observed, most adults contributed to genetically diverse progeny in experimental backwaters. These results support the use of predator-free, but otherwise natural, backwaters as an effective conservation tool for reintroducing Bonytail to its native habitat.

Molecular Genetics Unveiled Unknown Family Relationships and Hybrids in an Ex-Situ Colony of African Penguins (Spheniscus demersus)

Genealogical relationships among colony members, inbreeding status, and presence of hybrids are crucial data that can assist zoo curators in captive colony management and decision-making on relocation for reproduction. This study employed molecular markers to study a large colony (n = 56) of African Penguin hosted in an Italian biopark. A panel of 15 STRs (single tandem repeats) was selected, and genotype data were analyzed using COLONY software to determine parentage relationships and compare the existing studbook information to a pedigree built from genetic analyses. The existence of extra-pair mating and the presence of hybrids were investigated: discrepancies in kinship relationships emerged following molecular parentage analysis and 10 unknown genetic relationships were revealed. Infidelity of one member of the pair was observed in 6 cases and extra-pair copulation was assessed by genetic analysis in 2 episodes. One member of the colony was found to be a hybrid (S. demersus × S. humboldti); his progeny, derived by extra-pair copulation, was traced. Three other hidden hybrids were discovered and assessed using the identified candidate private alleles. Overall, our results demonstrate that molecular methods to confirm parentage and analyze relatedness among colony members are a valuable tool to complement studbook-based genetic management of African penguin captive populations. Because a variety of behavioral dynamics (e.g., extra-pair mating) can make observations ineffective in some species and because molecular markers outperform studbook in identifying the presence of hybrids, reliance on studbook information alone is not recommended.

Phenological match drives pollen-mediated gene flow in a temporally dimorphic tree

Variation in flowering phenology is common in natural populations, and is expected to be, together with inter-mate distance, an important driver of effective pollen dispersal. In populations composed of plants with temporally separated sexual phases (i.e. dichogamous or heterodichogamous populations), pollen-mediated gene flow is assumed to reflect phenological overlap between complementary sexual phases. In this study, we conducted paternity analyses to test this hypothesis in the temporally dimorphic tree Acer opalus. We performed spatially explicit analyses based on categorical and fractional paternity assignment, and included tree size, pair-wise genetic relatedness and morph type as additional predictors. Because differences between morphs in flowering phenology may also influence pollination distances, we modelled separate pollen dispersal kernels for the two morphs. Extended phenological overlap between male and female phases (mainly associated with inter-morph crosses) resulted in higher siring success after accounting for the effects of genetic relatedness, morph type and tree size, while reduced phenological overlap (mainly associated with intra-morph crosses) resulted in longer pollination distances achieved. Siring success also increased in larger trees. Mating patterns could not be predicted by phenology alone. However, as heterogeneity in flowering phenology was the single morph-specific predictor of siring success, it is expected to be key in maintaining the temporal dimorphism in A. opalus, by promoting not only a prevalent pattern of inter-morph mating, but also long-distance pollination resulting from intra-morph mating events.

Effects of Sample Size and Full Sibs on Genetic Diversity Characterization: A Case Study of Three Syntopic Iberian Pond-Breeding Amphibians

Accurate characterization of genetic diversity is essential for understanding population demography, predicting future trends and implementing efficient conservation policies. For that purpose, molecular markers are routinely developed for nonmodel species, but key questions regarding sampling design, such as calculation of minimum sample sizes or the effect of relatives in the sample, are often neglected. We used accumulation curves and sibship analyses to explore how these 2 factors affect marker performance in the characterization of genetic diversity. We illustrate this approach with the analysis of an empirical dataset including newly optimized microsatellite sets for 3 Iberian amphibian species: Hyla molleri, Epidalea calamita, and Pelophylax perezi. We studied 17-21 populations per species (total n = 547, 652, and 516 individuals, respectively), including a reference locality in which the effect of sample size was explored using larger samples (77-96 individuals). As expected, FIS and tests for Hardy-Weinberg equilibrium and linkage disequilibrium were affected by the presence of full sibs, and most initially inferred disequilibria were no longer statistically significant when full siblings were removed from the sample. We estimated that to obtain reliable estimates, the minimum sample size (potentially including full sibs) was close to 20 for expected heterozygosity, and between 50 and 80 for allelic richness. Our pilot study based on a reference population provided a rigorous assessment of marker properties and the effects of sample size and presence of full sibs in the sample. These examples illustrate the advantages of this approach to produce robust and reliable results for downstream analyses.

Demographic mechanisms underpinning genetic assimilation of remnant groups of a large carnivore

Current range expansions of large terrestrial carnivores are occurring following human-induced range contraction. Contractions are often incomplete, leaving small remnant groups in refugia throughout the former range. Little is known about the underlying ecological and evolutionary processes that influence how remnant groups are affected during range expansion. We used data from a spatially explicit, long-term genetic sampling effort of grizzly bears (Ursus arctos) in the Northern Continental Divide Ecosystem (NCDE), USA, to identify the demographic processes underlying spatial and temporal patterns of genetic diversity. We conducted parentage analysis to evaluate how reproductive success and dispersal contribute to spatio-temporal patterns of genetic diversity in remnant groups of grizzly bears existing in the southwestern (SW), southeastern (SE) and east-central (EC) regions of the NCDE. A few reproductively dominant individuals and local inbreeding caused low genetic diversity in peripheral regions that may have persisted for multiple generations before eroding rapidly (approx. one generation) during population expansion. Our results highlight that individual-level genetic and reproductive dynamics play critical roles during genetic assimilation, and show that spatial patterns of genetic diversity on the leading edge of an expansion may result from historical demographic patterns that are highly ephemeral.

Marine protected area restricts demographic connectivity: Dissimilarity in a marine environment can function as a biological barrier

The establishment of marine protected areas (MPAs) can often lead to environmental differences between MPAs and fishing zones. To determine the effects on marine dispersal of environmental dissimilarity between an MPA and fishing zone, we examined the abundance and recruitment patterns of two anemonefishes (Amphiprion frenatus and A. perideraion) that inhabit sea anemones in different management zones (i.e., an MPA and two fishing zones) by performing a field survey and a genetic parentage analysis. We found lower levels of abundance per anemone in the MPA compared to the fishing zones for both species (n = 1,525 anemones, p = .032). The parentage analysis also showed that lower numbers of fishes were recruited from the fishing zones and outside of the study area into each anemone in the MPA than into each anemone in the fishing zones (n = 1,525 anemones, p < .017). However, the number of self-recruit production per female did not differ between the MPA and fishing zones (n = 384 females, p = .516). Because the ocean currents around the study site were unlikely to cause a lower settlement intensity of larvae in the MPA, the ocean circulation was not considered crucial to the observed abundance and recruitment patterns. Instead, stronger top-down control and/or a lower density of host anemones in the MPA were potential factors for such patterns. Our results highlight the importance of dissimilarity in a marine environment as a factor that affects connectivity.

Microsatellite Mutation Rate in Atlantic Sturgeon (Acipenser oxyrinchus)

Understanding mutation rates can greatly extend the utility of population and conservation genetic analyses. Herein, we present an estimate of genome-wide microsatellite mutation rate in Atlantic sturgeon (Acipenser oxyrinchus) based on parent-offspring transmission patterns. We screened 307 individuals for parentage and mutation-rate analysis applying 43 variable markers. Out of 13228 allele transfers, 11 mutations were detected, producing a mutation rate of 8.3 × 10-4 per locus per generation (95% confidence interval: 1.48 × 10-3, 4.15 × 10-4). Single-step mutations predominated and there were trends toward mutations in loci with greater polymorphism and allele length. Two of the detected mutations were most probably cluster mutations, being identified in 12 and 28 sibs, respectively. Finally, we observed evidences of polyploidy based on the sporadic presence of 3 or 4 alleles per locus in the genotyped individuals, supporting previous reports of incomplete diploidization in Atlantic sturgeon.

Increase in multiple paternity across the reproductive lifespan in a sperm-storing, hermaphroditic freshwater snail

Polyandry is a common phenomenon and challenges the traditional view of stronger sexual selection in males than in females. In simultaneous hermaphrodites, the physical proximity of both sex functions was long thought to preclude the operation of sexual selection. Laboratory studies suggest that multiple mating and polyandry in hermaphrodites may actually be common, but data from natural populations are sparse. We therefore estimated the rate of multiple paternity and its seasonal variability in the annual, sperm-storing, simultaneously hermaphroditic freshwater snail Radix balthica for the entire duration of the reproductive lifespan. We also tested whether multiple paternity was associated with clutch size or embryonic development. To obtain these data, we measured and genotyped 60 field-collected egg clutches using nine highly polymorphic microsatellite markers. Overall, 50% of the clutches had multiple fathers, and both the frequency (20-93% of clutches) and magnitude of multiple paternity (mean 1.3-3.8 fathers per clutch) substantially increased over time, probably because of extensive sperm storage. Most multiply sired clutches (83%) had a dominant father, but neither clutch size nor the proportion of developed embryos per clutch was associated with levels of multiple paternity. Both the evident promiscuity and the frequent skew of paternity shares suggest that sexual selection may be an important evolutionary force in the study population.

Nature vs. Nurture: Evidence for Social Learning of Conflict Behaviour in Grizzly Bears

The propensity for a grizzly bear to develop conflict behaviours might be a result of social learning between mothers and cubs, genetic inheritance, or both learning and inheritance. Using non-invasive genetic sampling, we collected grizzly bear hair samples during 2011-2014 across southwestern Alberta, Canada. We targeted private agricultural lands for hair samples at grizzly bear incident sites, defining an incident as an occurrence in which the grizzly bear caused property damage, obtained anthropogenic food, or killed or attempted to kill livestock or pets. We genotyped 213 unique grizzly bears (118 M, 95 F) at 24 microsatellite loci, plus the amelogenin marker for sex. We used the program COLONY to assign parentage. We evaluated 76 mother-offspring relationships and 119 father-offspring relationships. We compared the frequency of problem and non-problem offspring from problem and non-problem parents, excluding dependent offspring from our analysis. Our results support the social learning hypothesis, but not the genetic inheritance hypothesis. Offspring of problem mothers are more likely to be involved in conflict behaviours, while offspring from non-problem mothers are not likely to be involved in incidents or human-bear conflicts themselves (Barnard's test, p = 0.05, 62.5% of offspring from problem mothers were problem bears). There was no evidence that offspring are more likely to be involved in conflict behaviour if their fathers had been problem bears (Barnard's test, p = 0.92, 29.6% of offspring from problem fathers were problem bears). For the mother-offspring relationships evaluated, 30.3% of offspring were identified as problem bears independent of their mother's conflict status. Similarly, 28.6% of offspring were identified as problem bears independent of their father's conflict status. Proactive mitigation to prevent female bears from becoming problem individuals likely will help prevent the perpetuation of conflicts through social learning.