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

Fine-scale functional connectivity of two syntopic pond-breeding amphibians with contrasting life-history traits: an integrative assessment of direct and indirect estimates of dispersal

Assessing patterns of functional connectivity among amphibian demes is crucial to unravel their population dynamics and prevent their isolation and eventual extinction. Integrative studies based on direct (capture-mark-recapture) and indirect (genetic) estimates of dispersal provide robust, biologically realistic inferences on population structure and connectivity, with applications for conservation efforts. We focused on two pond-breeding amphibians with contrasting life-history traits: the short-lived, semi-arboreal Hyla molleri and the long-lived, fossorial Pelobates cultripes. We PIT-tagged 2150 individuals of both species in two ponds (Laguna and Gravera, separated by 700 m) and monitored them from 2009 to 2021 to document the frequency and spatial extent of dispersal events. In addition, we genotyped individuals from these and two additional breeding populations at a maximum distance of 5 km with 15–16 microsatellites to characterize fine-scale patterns of genetic structure. We detected dispersal events connecting Laguna and Gravera in both species, albeit at low frequencies (4.8% and 7.7% of recaptured individuals of H. molleri and P. cultripes, respectively). However, both species were capable of covering long distances, with individual accumulated displacements up to 3.5 km (Hyla) and 1.8 km (Pelobates). Breeding populations > 2 km apart were genetically differentiated, indicating lower connectivity at this spatial scale. Estimates of pairwise migration rates differed between species and were asymmetrical, with different ponds representing “source” populations contributing more migrants to other populations in each species. We discuss the role of differences in life history traits and ecological preferences in shaping population dynamics in the two species and highlight management implications of our results.

A low-density single nucleotide polymorphism panel for brown trout (Salmo trutta L.) suitable for exploring genetic diversity at a range of spatial scales

The rivers of southern England and northern France which drain into the English Channel contain several genetically unique groups of trout (Salmo trutta L.) that have suffered dramatic declines in numbers over the past 40 years. Knowledge of levels and patterns of genetic diversity is essential for effective management of these vulnerable populations. Using restriction site-associated DNA sequencing (RADseq) data, we describe the development and characterisation of a panel of 95 single nucleotide polymorphism (SNP) loci for trout from this region and investigate their applicability and variability in both target (i.e., southern English) and non-target trout populations from northern Britain and Ireland. In addition, we present three case studies which demonstrate the utility and resolution of these genetic markers at three levels of spatial separation:(a) between closely related populations in nearby rivers, (b) within a catchment and (c) when determining parentage and familial relationships between fish sampled from a single site, using both empirical and simulated data. The SNP loci will be useful for population genetic and assignment studies on brown trout within the UK and beyond.

Atypical life cycle does not lead to inbreeding or selfing in parasites despite clonemate accumulation in intermediate hosts

Many parasites utilize asexual and sexual reproduction and multiple hosts to complete their life cycles. How these taxa avoid inbreeding is an essential question for understanding parasite evolution and ecology. Aquatic trematodes that require multiple host species may benefit from diverse genetic parasite assemblages accumulating within second intermediate hosts prior to sexual reproduction in definitive hosts. However, Cotylurus species are able to utilize the same snail species as first and second intermediate hosts, potentially resulting in the accumulation of genetically identical clones (clonemates) prior to sexual reproduction. In this study, we developed and analysed novel microsatellite loci to determine if clones are accumulating within snail hosts prior to ingestion by bird hosts and the effects this could have on parasite inbreeding. Contrary to previous studies of aquatic trematodes, significantly large numbers of clonemates were present within snails, but full-sibs were not. Genetic structure was present over a relatively small geographical scale despite the use of vagile definitive hosts. Phylogenetic analysis identified the Cotylurus sp. clones as belonging to a single species. Despite the presence of clones within snails, mating between clones/selfing was not common and heterozygosity is maintained within individuals. Potential issues with clones mating may be mitigated by the presence of snails with numerous clones, the consumption of many snails by bird hosts and parasite clone recognition/avoidance. Use of the same host species for multiple life stages may have advantages when parasites are able to avoid inbreeding and the required hosts are common.

Mitochondrial DNA diversity divulges high levels of haplotype diversity and lack of genetic structure in the Indian camels

Camels represent an important genetic resource of the desert ecosystems of India, with the dromedary and Bactrian camels inhabiting the hot and cold deserts, respectively. This study is the first attempt to investigate mitochondrial DNA based genetic diversity in the Indian camel populations and explores their relationship in the context of global genetic diversity of all the three large camel species (Camelus ferus, Camelus bactrianus and Camelus dromedaries). A mitochondrial DNA fragment encompassing part of cytochrome b gene, tRNA^Thr, tRNA^Pro and the beginning of the control region was amplified and analyzed in 72 dromedary and 8 Bactrian camels of India. Sequence analysis revealed that the haplotype and nucleotide diversity (Hd: 0.937 and π: 0.00431) in the Indian dromedaries was higher than the indices reported so far for the dromedary or Bactrian camels across the globe. The corresponding values in the Indian Bactrian camels were 1.000 and 0.00393, respectively. Signals of population expansion were evident in the dromedaries of India on the basis of mismatch analysis and Fu's Fs values. The analysis of molecular variance attributed most of the genetic variance (92.15%) between the dromedary, wild Bactrian and domestic Bactrian camels indicating separate maternal origins. The existence of three mitochondrial lineages in the old world camels (C. bactrianus: Lineage A; C. ferus: Lineage B and C. dromedarius: Lineage C) was also substantiated by the topology of the Median-Joining network.

Genetic Diversity, Admixture and Analysis of Homozygous-by-Descent (HBD) Segments of Russian Wild Boar

The wild boar is the wild ancestor of the domestic pig and one of the most common species of ungulates. At the beginning of the 20th century, the wild boar was practically exterminated in the European part of Russia. In the period 1935-1988, 7705 boars were caught in various regions of the European part of Russia, the Far East, Ukraine, Belarus, Kyrgyzstan, Kazakhstan, Latvia, Lithuania, Estonia, Tajikistan and resettled in the territory of Russia. Asian and European wild boars dwell the territory of Russia. The aim of our research was to study the genetic diversity and structure of wild boar populations in different regions of Russia using genome-wide genotyping. We have determined the genetic distances, population structure, parameters of genetic diversity and significantly expanded our understanding of the genetic state of the Russian wild boar. For the first time, we calculated autozygosity of the wild boar of the European and Asian subspecies using Homozygous-by-Descent (HBD) Segments analysis, which is important in terms of population recovery. We also found evidence of hybridization between Russian wild boar and domestic pigs. A group of European wild boars showed introgression of the Asian boar into population. The mean level of the inbreeding coefficient in European wild boar was higher than in Asian wild boar, and combined groups of the European boar had higher inbreeding coefficient than Russian wild boars. These results obtained can be used in population management.

Patterns of Genetic Diversity and Mating Systems in a Mass-Reared Black Soldier Fly Colony

The black soldier fly (BSF), Hermetia illucens, is a promising candidate for the emerging insect farming industry with favourable characteristics for both bioremediation and production of animal delivered nutritive and industrial compounds. The genetic management of commercial colonies will become increasingly important for the sustainability of the industry. However, r-selected life history traits of insects pose challenges to conventional animal husbandry and breeding approaches. In this study, the long-term genetic effects of mass-rearing were evaluated as well as mating systems in the species to establish factors that might influence genetic diversity, and by implication fitness and productivity in commercial colonies. Population genetic parameters, based on microsatellite markers, were estimated and compared amongst two temporal wild sampling populations and four generations (F28, F48, F52, and F62) of a mass-reared colony. Furthermore, genetic relationships amongst mate pairs were evaluated and parentage analysis was performed to determine the oc-currence of preferential mate choice and multiple paternity. The mass-reared colony showed a reduction in genetic diversity and evidence for inbreeding with significant successive generational genetic differentiation from the wild progenitor population. Population-level analysis also gave the first tentative evidence of positive assortative mating and genetic polyandry in BSF. The homoge-neity of the mass-reared colony seems to result from a dual action caused by small effective popu-lation size and increased homozygosity due to positive assortative mating. However, the high ge-netic diversity in the wild and a polyandrous mating system might suggest the possible restoration of diversity in mass-reared colonies through augmentation with the wild population.

Evaluating surrogates of genetic diversity for conservation planning

Protected-area systems should conserve intraspecific genetic diversity. Because genetic data require resources to obtain, several approaches have been proposed for generating plans for protected-area systems (prioritizations) when genetic data are not available. Yet such surrogate-based approaches remain poorly tested. We evaluated the effectiveness of potential surrogate-based approaches based on microsatellite genetic data collected across the Iberian Peninsula for 7 amphibian and 3 reptilian species. Long-term environmental suitability did not effectively represent sites containing high genetic diversity (allelic richness). Prioritizations based on long-term environmental suitability had similar performance to random prioritizations. Geographic distances and resistance distances based on contemporary environmental suitability were not always effective surrogates for identification of combinations of sites that contain individuals with different genetic compositions. Our results demonstrate that population genetic data based on commonly used neutral markers can inform prioritizations, and we could not find an adequate substitute. Conservation planners need to weigh the potential benefits of genetic data against their acquisition costs.

Genetic Diversity and Structure of Common Carp (Cyprinus carpio L.) in the Centre of Carpathian Basin: Implications for Conservation

Hungary is one of the largest common carp-production countries in Europe and now, there is a large number of local breeds and strains in the country. For proper maintenance of the animal genetic resources, information on their genetic diversity and structure is essential. At present, few data are available on the genetic purity and variability of the Hungarian common carp. In this study, we genetically analyzed 13 strains in Hungary and, in addition, the Amur wild carp, using 12 microsatellite markers. A total of 117 unique alleles were detected in 630 individuals. Low levels of genetic differentiation (Fst and Cavalli-Sforza and Edwards distance) were estimated among strains. The AMOVA showed the low but significant level of genetic differentiation among strains (3.79%). Bayesian clustering analysis using STRUCTURE classified the strains into 14 different clusters. The assignment test showed that 93.64% of the individuals could be assigned correctly into their original strain. Overall, our findings can be contributed to complementing scientific knowledge for conservation and management of threatened strains of common carp.

Evaluation of the Minimum Sampling Design for Population Genomic and Microsatellite Studies: An Analysis Based on Wild Maize

Massive parallel sequencing (MPS) is revolutionizing the field of molecular ecology by allowing us to understand better the evolutionary history of populations and species, and to detect genomic regions that could be under selection. However, the economic and computational resources needed generate a tradeoff between the amount of loci that can be obtained and the number of populations or individuals that can be sequenced. In this work, we analyzed and compared two simulated genomic datasets fitting a hierarchical structure, two extensive empirical genomic datasets, and a dataset comprising microsatellite information. For all datasets, we generated different subsampling designs by changing the number of loci, individuals, populations, and individuals per population to test for deviations in classic population genetics parameters (H_S, F_IS, F_ST). For the empirical datasets we also analyzed the effect of sampling design on landscape genetic tests (isolation by distance and environment, central abundance hypothesis). We also tested the effect of sampling a different number of populations in the detection of outlier SNPs. We found that the microsatellite dataset is very sensitive to the number of individuals sampled when obtaining summary statistics. F_IS was particularly sensitive to a low sampling of individuals in the simulated, genomic, and microsatellite datasets. For the empirical and simulated genomic datasets, we found that as long as many populations are sampled, few individuals and loci are needed. For the empirical datasets, we found that increasing the number of populations sampled was important in obtaining precise landscape genetic estimates. Finally, we corroborated that outlier tests are sensitive to the number of populations sampled. We conclude by proposing different sampling designs depending on the objectives.

Contrasting demographic trends and asymmetric migration rates in a spatially structured amphibian population

Natural populations often persist at the landscape scale as metapopulations, with breeding units (subpopulations) experiencing temporal extinction and recolonization events. Important parameters to forecast population viability in these systems include the ratio of the effective number of breeders (N_b ) to the total number of adults (N_a ) and migration rates among subpopulations. Here, we present the results of a 10-year integrative monitoring program of a metapopulation of the Iberian green frog (Pelophylax perezi) in central Spain. We characterized population dynamics at two main breeding ponds (Gravera and Laguna) using capture-mark-recapture data to estimate N_a in each breeding season, and multilocus genotypes to estimate the effective population size (N_e ), N_b , individual breeding success, and migration rates. Both ponds experienced population decline after a dry season, with Gravera subsequently recovering and Laguna suffering a bottleneck associated with genetic impoverishment. In this subpopulation, average allelic richness and private alleles decreased from 2010 (10.87 and 1.67, respectively) to 2018 (8.0 and 0.20). The N_b /N_a ratio in Laguna in 2018 was twice as high (0.95) than in Gravera (0.41) or in pre-bottleneck Laguna (0.50), suggesting plasticity or genetic compensation through increased individual breeding success. Migration rates were asymmetric between ponds, with a stronger contribution from Gravera to Laguna (29.9% vs. 16.2% in the opposite direction) that may result in a rescue effect. This study emphasizes the importance of integrative demographic approaches for the monitoring of natural populations based on a better understanding of their spatio-temporal dynamics, which provides valuable information for conservation actions.

'Mainland-island' population structure of a terrestrial salamander in a forest-bocage landscape with little evidence for in situ ecological speciation

Adaptation to different ecological environments can, through divergent selection, generate phenotypic and genetic differences between populations, and eventually give rise to new species. The fire salamander (Salamandra salamandra) has been proposed to represent an early stage of ecological speciation, driven by differential habitat adaptation through the deposition and development of larvae in streams versus ponds in the Kottenforst near Bonn (Germany). We set out to test this hypothesis of ecological speciation in an area different from the one where it was raised and we took the opportunity to explore for drivers of genetic differentiation at a landscape scale. A survey over 640 localities demonstrated the species' presence in ponds and streams across forests, hilly terrain and areas with hedgerows ('bocage'). Genetic variation at 14 microsatellite loci across 41 localities in and around two small deciduous forests showed that salamander effective population sizes were higher in forests than in the bocage, with panmixia in the forests (Fst < 0.010) versus genetic drift or founder effects in several of the small and more or less isolated bocage populations (Fst > 0.025). The system fits the 'mainland-island' metapopulation model rather than indicating adaptive genetic divergence in pond versus stream larval habitats. A reanalysis of the Kottenforst data indicated that microsatellite genetic variation fitted a geographical rather than an environmental axis, with a sharp transition from a western pond-breeding to an eastern, more frequently stream-breeding group of populations. A parallel changeover in mitochondrial DNA exists but remains to be well documented. The data support the existence of a hybrid zone following secondary contact of differentiated lineages, more so than speciation in situ.

The Genomic Makeup of Nine Horse Populations Sampled in the Netherlands

The spectrum of modern horse populations encompasses populations with a long history of development in isolation and relatively recently formed types. To increase our understanding of the evolutionary history and provide information on how to optimally conserve or improve these populations with varying development and background for the future, we analyzed genotype data of 184 horses from 9 Dutch or common horse populations in the Netherlands: The Belgian draft horse, Friesian horse, Shetland pony, Icelandic horse, Gelder horse, Groninger horse, harness horse, KWPN sport horse and the Lipizzaner horse population. Various parameters were estimated (e.g., runs of homozygosity and F_ST values) to gain insight into genetic diversity and relationships within and among these populations. The identified genomic makeup and quantified relationships did mostly conform to the development of these populations as well as past and current breeding practices. In general, populations that allow gene-flow showed less inbreeding and homozygosity. Also, recent bottlenecks (e.g., related to high selective pressure) caused a larger contribution of long ROHs to inbreeding. Maintaining genetic diversity through tailor-made breeding practices is crucial for a healthy continuation of the investigated, mostly inbred and (effectively) small sized horse populations, of which several already experience inbreeding related issues.

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

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

Reliable effective number of breeders/adult census size ratios in seasonal-breeding species: Opportunity for integrative demographic inferences based on capture-mark-recapture data and multilocus genotypes

The ratio of the effective number of breeders (Nb) to the adult census size (Na), Nb/Na, approximates the departure from the standard capacity of a population to maintain genetic diversity in one reproductive season. This information is relevant for assessing population status, understanding evolutionary processes operating at local scales, and unraveling how life-history traits affect these processes. However, our knowledge on Nb/Na ratios in nature is limited because estimation of both parameters is challenging. The sibship frequency (SF) method is adequate for reliable Nb estimation because it is based on sibship and parentage reconstruction from genetic marker data, thereby providing demographic inferences that can be compared with field-based information. In addition, capture-mark-recapture (CMR) robust design methods are well suited for Na estimation in seasonal-breeding species. We used tadpole genotypes of three pond-breeding amphibian species (Epidalea calamita, Hyla molleri, and Pelophylax perezi, n = 73-96 single-cohort tadpoles/species genotyped at 15-17 microsatellite loci) and candidate parental genotypes (n = 94-300 adults/species) to estimate Nb by the SF method. To assess the reliability of Nb estimates, we compared sibship and parentage inferences with field-based information and checked for the convergence of results in replicated subsampled analyses. Finally, we used CMR data from a 6-year monitoring program to estimate annual Na in the three species and calculate the Nb/Na ratio. Reliable ratios were obtained for E. calamita (Nb/Na = 0.18-0.28) and P. perezi (0.5), but in H. molleri, Na could not be estimated and genetic information proved insufficient for reliable Nb estimation. Integrative demographic studies taking full advantage of SF and CMR methods can provide accurate estimates of the Nb/Na ratio in seasonal-breeding species. Importantly, the SF method provides results that can be readily evaluated for reliability. This represents a good opportunity for obtaining robust demographic inferences with wide applications for evolutionary and conservation research.