Interannual variation in effective number of breeders and estimation of effective population size in long‐lived iteroparous lake sturgeon (Acipenser fulvescens)

Quantifying interannual variation in effective adult breeding number (N_b) and relationships between N_b, effective population size (N_e), adult census size (N) and population demographic characteristics are important to predict genetic changes in populations of conservation concern. Such relationships are rarely available for long‐lived iteroparous species like lake sturgeon (Acipenser fulvescens). We estimated annual N_b and generational N_e using genotypes from 12 microsatellite loci for lake sturgeon adults (n = 796) captured during ten spawning seasons and offspring (n = 3925) collected during larval dispersal in a closed population over 8 years. Inbreeding and variance N_b estimated using mean and variance in individual reproductive success derived from genetically identified parentage and using linkage disequilibrium (LD) were similar within and among years (interannual range of N_b across estimators: 41–205). Variance in reproductive success and unequal sex ratios reduced N_b relative to N on average 36.8% and 16.3%, respectively. Interannual variation in N_b/N ratios (0.27–0.86) resulted from stable N and low standardized variance in reproductive success due to high proportions of adults breeding and the species' polygamous mating system, despite a 40‐fold difference in annual larval production across years (437–16 417). Results indicated environmental conditions and features of the species' reproductive ecology interact to affect demographic parameters and N_b/N. Estimates of N_e based on three single‐sample estimators, including LD, approximate Bayesian computation and sibship assignment, were similar to annual estimates of N_b. Findings have important implications concerning applications of genetic monitoring in conservation planning for lake sturgeon and other species with similar life histories and mating systems.     

Genetic population structure and demography of an apex predator,the tiger shark Galeocerdo cuvier

Population genetics has been increasingly applied to study large sharks over the last decade. Whilst large shark species are often difficult to study with direct methods, improved knowledge is needed for both population management and conservation, especially for species vulnerable to anthropogenic and climatic impacts. The tiger shark, Galeocerdo cuvier, is an apex predator known to play important direct and indirect roles in tropical and subtropical marine ecosystems. While the global and Indo‐West Pacific population genetic structure of this species has recently been investigated, questions remain over population structure and demographic history within the western Indian (WIO) and within the western Pacific Oceans (WPO). To address the knowledge gap in tiger shark regional population structures, the genetic diversity of 286 individuals sampled in seven localities was investigated using 27 microsatellite loci and three mitochondrial genes (CR, COI, and cytb). A weak genetic differentiation was observed between the WIO and the WPO, suggesting high genetic connectivity. This result agrees with previous studies and highlights the importance of the pelagic behavior of this species to ensure gene flow. Using approximate Bayesian computation to couple information from both nuclear and mitochondrial markers, evidence of a recent bottleneck in the Holocene (2,000–3,000 years ago) was found, which is the most probable cause for the low genetic diversity observed. A contemporary effective population size as low as 111 [43,369] was estimated during the bottleneck. Together, these results indicate low genetic diversity that may reflect a vulnerable population sensitive to regional pressures. Conservation measures are thus needed to protect a species that is classified as Near Threatened.     

New approaches to species delimitation and population structure of anthozoans: Two case studies of octocorals using ultraconserved elements and exons

As coral populations decline worldwide in the face of ongoing environmental change, documenting their distribution, diversity and conservation status is now more imperative than ever. Accurate delimitation and identification of species is a critical first step. This task, however, is not trivial as morphological variation and slowly evolving molecular markers confound species identification. New approaches to species delimitation in corals are needed to overcome these challenges. Here, we test whether target enrichment of ultraconserved elements (UCEs) and exons can be used for delimiting species boundaries and population structure within species of corals by focusing on two octocoral genera, Alcyonium and Sinularia, as exemplary case studies. We designed an updated bait set (29,181 baits) to target‐capture 3,023 UCE and exon loci, recovering a mean of 1,910 ± 168 SD per sample with a mean length of 1,055 ± 208 bp. Similar numbers of loci were recovered from Sinularia (1,946 ± 227 SD) and Alcyonium (1,863 ± 177 SD). Species‐level phylogenies were highly supported for both genera. Clustering methods based on filtered single nucleotide polymorphisms delimited species and populations that are congruent with previous allozyme, DNA barcoding, reproductive and ecological data for Alcyonium, and offered further evidence of hybridization among species. For Sinularia, results were congruent with those obtained from a previous study using restriction site associated DNA sequencing. Both case studies demonstrate the utility of target‐enrichment of UCEs and exons to address a wide range of evolutionary and taxonomic questions across deep to shallow timescales in corals.     

Mature male parr contribution to the effective size of an anadromous Atlantic salmon (Salmo salar) population over 30 years

We describe temporal changes in the genetic composition of a small anadromous Atlantic salmon (Salmo salar) population from South Newfoundland, an area where salmon populations are considered threatened (COSEWIC 2010). We examined the genetic variability (13 microsatellite loci) in 869 out‐migrating smolt and post‐spawning kelt samples, collected from 1985 to 2011 for a total of 22 annual collections and a 30 year span of assigned cohorts. We estimated the annual effective number of breeders (N_b) and the generational effective population size (N_e) through genetic methods and demographically using the adult sex ratio. Comparisons between genetic and demographic estimates show that the adult spawners inadequately explain the observed N_e estimates, suggesting that mature male parr are significantly increasing N_b and N_e over the study period. Spawning as parr appears to be a viable and important strategy in the near absence of adult males.     

A genetic demographic analysis of Lake Malawi rock‐dwelling cichlids using spatio‐temporal sampling

We estimated the effective population sizes (N_e) and tested for short‐term temporal demographic stability of populations of two Lake Malawi cichlids: Maylandia benetos, a micro‐endemic, and Maylandia zebra, a widespread species found across the lake. We sampled a total of 351 individuals, genotyped them at 13 microsatellite loci and sequenced their mitochondrial D‐loop to estimate genetic diversity, population structure, demographic history and effective population sizes. At the microsatellite loci, genetic diversity was high in all populations. Yet, genetic diversity was relatively low for the sequence data. Microsatellites yielded mean N_e estimates of 481 individuals (±99 SD) for M. benetos and between 597 (±106.3 SD) and 1524 (±483.9 SD) individuals for local populations of M. zebra. The microsatellite data indicated no deviations from mutation–drift equilibrium. Maylandia zebra was further found to be in migration–drift equilibrium. Temporal fluctuations in allele frequencies were limited across the sampling period for both species. Bayesian Skyline analyses suggested a recent expansion of M. zebra populations in line with lake‐level fluctuations, whereas the demographic history of M. benetos could only be estimated for the very recent past. Divergence time estimates placed the origin of M. benetos within the last 100 ka after the refilling of the lake and suggested that it split off the sympatric M. zebra population. Overall, our data indicate that micro‐endemics and populations in less favourable habitats have smaller N_e, indicating that drift may play an important role driving their divergence. Yet, despite small population sizes, high genetic variation can be maintained.     

Disentangling the impact of demographic factors on population differentiation of an endangered freshwater crayfish (Austropotamobius pallipes) using population density and microsatellite data

1. Habitat fragmentation of stream ecosystems often results in decreased connectivity between populations and lower population sizes. Hence, understanding how habitat fragmentation affects genetic erosion is important for the preservation of freshwater biodiversity, in particular, as small populations suffer from loss of genetic diversity through genetic drift and loss of fitness because of inbreeding, increasing the risk of extinction. 2. Here, we assess the impact of demographic factors on population differentiation in the endangered freshwater crayfish Austropotamobius pallipes by analysing population genetic structure, estimating effective population sizes and comparing levels of polymorphism at five microsatellite loci with density estimates of 10 populations within a small French catchment that has become progressively confined to headwaters over the last six decades. 3. Levels of expected heterozygosity and allelic richness per population were relatively low (0.214–0.396 and 1.6–2.6, respectively). We found strong genetic differentiation between these geographically close populations (F_ST = 0.283), with weak statistical evidence for a pattern of isolation by distance. Estimates of effective population size were low (<150) in most populations, but potentially reached several thousands in three populations. 4. Population density and allelic richness were strongly positively correlated. A robust relationship between population density and heterozygosity values was also noted, but only after discarding two populations for which significant genetic signatures of a recent bottleneck were found; these two populations displayed high expected heterozygosity compared with a very low density. Populations with the highest densities of individuals had the highest effective population size estimates and vice versa. 5. Our results clearly show the importance of demographic factors and genetic drift on A. pallipes populations. Furthermore, analysis of genetic and population density data is a pragmatic and efficient approach to corroborate inferences from genetic data and can be particularly useful in the identification of populations experiencing a bottleneck and therefore in conservation genetics studies aiming at identifying priority populations for conservation.     

Nucleotide diversity patterns of three divergent soybean populations: evidences for population‐dependent linkage disequilibrium and taxonomic status of Glycine gracilis

The level of linkage disequilibrium (LD) is a major factor to determine DNA polymorphism pattern of a population and to construct high‐resolution maps useful in localizing and gene cloning of complicated traits. Here, we investigated LD level of three soybean populations with different genetic backgrounds and taxonomic status of G. gracilis by comparing the DNA polymorphism patterns of four high‐diversity single‐copy nuclear genes. A total of 152, 22, and 77 accessions of G. soja, G. gracilis, and G. max were observed. The results indicated that G. max retained only 75.3 (π) and 39% (θ) of the nucleotide polymorphism found in G. soja. Four gene loci evolved according to neutrality in both G. max and G. gracilis populations, and three gene loci evolved according to neutrality in G. soja population by Tajima's and Fu and Li's test. However, one gene locus deviated from neutrality by Fu and Li's test in the G. soja population. Further, medial level of LD (average r² = 0.2426) was found in intragene in G. max and G. gracilis populations, but unexpected low level of LD (r² ≤ 0.0539) was found in G. soja population. Significant genetic differentiation was detected between G. max and G. soja populations and also between G. max and G. gracilis populations; however, nonsignificant genetic differentiation was found between G. gracilis and G. soja populations. The results suggest that LD level depends on genetic background of soybean population, and implicit that G. gracilis should be regarded as the variant of G. soja, not as an independent species.     

Considerations for monitoring population trends of colonial waterbirds using the effective number of breeders and census estimates

Detecting trends in population size fluctuations is a major focus in ecology, evolution, and conservation biology. Populations of colonial waterbirds have been monitored using demographic approaches to determine annual census size (N_a). We propose the addition of genetic estimates of the effective number of breeders (N_b) as indirect measures of the risk of loss of genetic diversity to improve the evaluation of demographics and increase the accuracy of trend estimates in breeding colonies. Here, we investigated which methods of the estimation of N_b are more precise under conditions of moderate genetic diversity, limited sample sizes and few microsatellite loci, as often occurs with natural populations. We used the wood stork as a model species and we offered a workflow that researchers can follow for monitoring bird breeding colonies. Our approach started with simulations using five estimators of N_b and the theoretical results were validated with empirical data collected from breeding colonies settled in the Brazilian Pantanal wetland. In parallel, we estimated census size using a corrected method based on counting active nests. Both in simulations and in natural populations, the approximate Bayesian computation (ABC) and sibship assignment (SA) methods yielded more precise estimates than the linkage disequilibrium, heterozygosity excess, and molecular coancestry methods. In particular, the ABC method performed best with few loci and small sample sizes, while the other estimators required larger sample sizes and at least 13 loci to not underestimate N_b. Moreover, according to our N_b/N_a estimates (values were often ≤0.1), the wood stork colonies evaluated could be facing the loss of genetic diversity. We demonstrate that the combination of genetic and census estimates is a useful approach for monitoring natural breeding bird populations. This methodology has been recommended for populations of rare species or with a known history of population decline to support conservation efforts.     

Environmental and topographic variables shape genetic structure and effective population sizes in the endangered Yosemite toad

Aim Describing the landscape variables that accurately reflect how environmental and topographic variations affect population connectivity and demography is a major goal of landscape genetics and conservation biology. However, few landscape genetics studies have quantified the relationships between landscape variables and effective population size (N_e), although N_e is a key conservation and population genetics parameter. In this study, I estimated genetic structure and effective population sizes in the Yosemite toad (Bufo canorus) and tested for associations with environmental and geographic variables. Location Yosemite National Park, California, USA. Methods I estimated F_ST, D_ps and N_e using 10 microsatellite loci amplified from 781 individuals from 24 populations. I used three landscape variables (environmental variation, topography and slope) to generate geographic distance models and a series of regression analyses to identify the variables that contributed to genetic structure in this species. I also tested for correlations between N_e and a suite of variables, including geographic and genetic isolation, habitat suitability, elevation, temperature and precipitation. Results I found substantial variation in genetic distances between populations (F_ST = 0.004–0.396, D_ps = 0.045–0.839) and in effective population sizes (N_e = 9–52). Environmental variation and slope played important roles in explaining variation in genetic distances, and precipitation variables were significantly correlated with N_e. Main conclusions These results show that environmental and topographic variables are both important for understanding population connectivity in B. canorus and provide some of the first evidence, in any species, for a link between environmental variables and effective population size.     

Temporal dynamics of linkage disequilibrium in two populations of bighorn sheep

Linkage disequilibrium (LD) is the nonrandom association of alleles at two markers. Patterns of LD have biological implications as well as practical ones when designing association studies or conservation programs aimed at identifying the genetic basis of fitness differences within and among populations. However, the temporal dynamics of LD in wild populations has received little empirical attention. In this study, we examined the overall extent of LD, the effect of sample size on the accuracy and precision of LD estimates, and the temporal dynamics of LD in two populations of bighorn sheep (Ovis canadensis) with different demographic histories. Using over 200 microsatellite loci, we assessed two metrics of multi‐allelic LD, D′, and χ^′2. We found that both populations exhibited high levels of LD, although the extent was much shorter in a native population than one that was founded via translocation, experienced a prolonged bottleneck post founding, followed by recent admixture. In addition, we observed significant variation in LD in relation to the sample size used, with small sample sizes leading to depressed estimates of the extent of LD but inflated estimates of background levels of LD. In contrast, there was not much variation in LD among yearly cross‐sections within either population once sample size was accounted for. Lack of pronounced interannual variability suggests that researchers may not have to worry about interannual variation when estimating LD in a population and can instead focus on obtaining the largest sample size possible.     

Accounting for cryptic population substructure enhances detection of inbreeding depression with genomic inbreeding coefficients: an example from a critically endangered marsupial

Characterizing inbreeding depression in wildlife populations can be critical to their conservation. Coefficients of individual inbreeding can be estimated from genome‐wide marker data. The degree to which sensitivity of inbreeding coefficients to population genetic substructure alters estimates of inbreeding depression in wild populations is not well understood. Using generalized linear models, we tested the power of two frequently used inbreeding coefficients that are calculated from genome‐wide SNP markers, F_H and F^^III, to predict four fitness traits estimated over two decades in an isolated population of the critically endangered Leadbeater's possum. F_H estimates inbreeding as excess observed homozygotes relative to equilibrium expectations, whereas F^^III quantifies allelic similarity between the gametes that formed an individual, and upweights rare homozygotes. We estimated F_H and F^^III from 1,575 genome‐wide SNP loci in individuals with fitness trait data (N = 179–237 per trait), and computed revised coefficients, F_H^{by group} and F^^{IIIby group}, adjusted for population genetic substructure by calculating them separately within two different genetic groups of individuals identified in the population. Using F_H or F^^III in the models, inbreeding depression was detected for survival to sexual maturity, longevity and whether individuals bred during their lifetime. F^^{IIIby group} (but not F_H^{by group}) additionally revealed significant inbreeding depression for lifetime reproductive output (total offspring assigned to each individual). Estimates of numbers of lethal equivalents indicated substantial inbreeding load, but differing between inbreeding estimators. Inbreeding depression, declining population size, and low and declining genetic diversity suggest that genetic rescue may assist in preventing extinction of this unique Leadbeater's possum population.     

Estimating effective population size from linkage disequilibrium: severe bias in small samples

Effective population size (N _e) is a central concept in evolutionary biology and conservation genetics. It predicts rates of loss of neutral genetic variation, fixation of deleterious and favourable alleles, and the increase of inbreeding experienced by a population. A method exists for the estimation of N _e from the observed linkage disequilibrium between unlinked loci in a population sample. While an increasing number of studies have applied this method in natural and managed populations, its reliability has not yet been evaluated. We developed a computer program to calculate this estimator of N _e using the most widely used linkage disequilibrium algorithm and used simulations to show that this estimator is strongly biased when the sample size is small (<‰100) and below the true N _e. This is probably due to the linkage disequilibrium generated by the sampling process itself and the inadequate correction for this phenomenon in the method. Results suggest that N _e estimates derived using this method should be regarded with caution in many cases. To improve the method’s reliability and usefulness we propose a way to determine whether a given sample size exceeds the population N _e and can therefore be used for the computation of an unbiased estimate.     

DIRECT AND INDIRECT ESTIMATES OF NEIGHBORHOOD AND EFFECTIVE POPULATION SIZE IN A TROPICAL PALM,ASTROCARYUM MEXICANUM

To estimate the relative importance of genetic drift, the effective population size ???(N_e) can be used. Here we present estimates of the effective population size and related measures in Astrocaryum mexicanum, a tropical palm from Los Tuxtlas rain forest, Veracruz, Mexico. Seed and pollen dispersal were measured. Seeds are primarily dispersed by gravity and secondarily dispersed by small mammals. Mean primary and secondary dispersal distances for seeds were found to be small (0.78 m and 2.35 m, respectively). A. mexicanum is beetle pollinated and pollen movements were measured by different methods: a) using fluorescent dyes, b) as the minimum distance between active female and male inflorescences, and c) using rare allozyme alleles as genetic markers. All three estimates of pollen dispersal were similar, with a mean of approximately 20 m. Using the seed and pollen dispersal data, the genetic neighborhood area (A) was estimated to be 2,551 m². To obtain the effective population size, three different overlapping generation methods were used to estimate an effective density with demographic data from six permanent plots. The effective density ranged from 0.040 to 0.351 individuals per m². The product of effective density and neighborhood area yields a direct estimate of the neighborhood effective population size (N_b). N_b ranged from 102 to 895 individuals. Indirect estimates of population size and migration rate (Nm) were obtained using F_st for five different allozymic loci for both adults and seeds. We obtained a range of Nm from 1.2 to 19.7 in adults and a range of Nm from 4.0 to 82.6 for seeds. We discuss possible causes of the smaller indirect estimates of Nm relative to the direct and compare our estimates with values from other plant populations. Gene dispersal distances, neighborhood size, and effective population size in A. mexicanum are relatively high, suggesting that natural selection, rather than genetic drift, may play a dominant role in patterning the genetic variation in this tropical palm.     

SNP panel development for genetic management of wild and domesticated white bass (Morone chrysops)

White bass (Morone chrysops), striped bass and their interspecific hybrid are important game fishes, whereas the hybrid striped bass is an important aquaculture species in the US. Numerous state, federal and private hatcheries, therefore, rear these species for stocking purposes as well as for food fish. Although striped bass populations (both wild and domesticated) have been extensively evaluated, relatively little effort has been directed toward the study and improvement of white bass. In this study, we developed SNP resources to examine the genetic relationships among a long‐term domesticated white bass line and five potential founder stocks for selective breeding collected from drainages in Arkansas, Texas and Alabama. Using genotyping‐by‐sequencing, we generated 13 872 genome‐wide SNP loci across the six populations. Stringent filtering of SNP‐calling parameters identified 426 informative SNP loci. Population genetic and structure analyses using these loci revealed only moderate genetic differentiation between populations (global F_st = 0.083) and indicated two major genetic clusters. A final 57‐SNP assay was successfully designed and validated using the MassARRAY system. The developed SNP panel assigned 96 additional genotyped individuals to their population of origin with 100% accuracy. The SNP resources developed in this study should facilitate ongoing efforts in selective breeding and conservation of white bass.     

Genetic structure and gene flow among Komodo dragon populations inferredby microsatellite loci analysis

A general concern for the conservation of endangered species is the maintenance of genetic variation within populations, particularly when they become isolated and reduced in size. Estimates of gene flow and effective population size are therefore important for any conservation initiative directed to the long-term persistence of a species in its natural habitat. In the present study, 10 microsatellite loci were used to assess the level of genetic variability among populations of the Komodo dragon Varanus komodoensis. Effective population size was calculated and gene flow estimates were compared with palaeogeographic data in order to assess the degree of vulnerability of four island populations. Rinca and Flores, currently separated by an isthmus of about 200 m, retained a high level of genetic diversity and showed a high degree of genetic similarity, with gene flow values close to one migrant per generation. The island of Komodo showed by far the highest levels of genetic divergence, and its allelic distinctiveness was considered of great importance in the maintenance of genetic variability within the species. A lack of distinct alleles and low levels of gene flow and genetic variability were found for the small population of Gili Motang island, which was identified as vulnerable to stochastic threats. Our results are potentially important for both the short- and long-term management of the Komodo dragon, and are critical in view of future re-introduction or augmentation in areas where the species is now extinct or depleted.     

Global population genetic dynamics of a highly migratory,apex predator shark

Knowledge of genetic connectivity dynamics in the world's large‐bodied, highly migratory, apex predator sharks across their global ranges is limited. One such species, the tiger shark (Galeocerdo cuvier), occurs worldwide in warm temperate and tropical waters, uses remarkably diverse habitats (nearshore to pelagic) and possesses a generalist diet that can structure marine ecosystems through top‐down processes. We investigated the phylogeography and the global population structure of this exploited, phylogenetically enigmatic shark by using 10 nuclear microsatellites (n = 380) and sequences from the mitochondrial control region (CR, n = 340) and cytochrome oxidase I gene (n = 100). All three marker classes showed the genetic differentiation between tiger sharks from the western Atlantic and Indo‐Pacific ocean basins (microsatellite F_ST > 0.129; CR Φ_ST > 0.497), the presence of North vs. southwestern Atlantic differentiation and the isolation of tiger sharks sampled from Hawaii from other surveyed locations. Furthermore, mitochondrial DNA revealed high levels of intraocean basin matrilineal population structure, suggesting female philopatry and sex‐biased gene flow. Coalescent‐ and genetic distance‐based estimates of divergence from CR sequences were largely congruent (d_corr = 0.0015–0.0050), indicating a separation of Indo‐Pacific and western Atlantic tiger sharks <1 million years ago. Mitochondrial haplotype relationships suggested that the western South Atlantic Ocean was likely a historical connection for interocean basin linkages via the dispersal around South Africa. Together, the results reveal unexpectedly high levels of population structure in a highly migratory, behaviourally generalist, cosmopolitan ocean predator, calling for management and conservation on smaller‐than‐anticipated spatial scales.     

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction‐site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small N_e (range = 2.1–89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome‐wide divergence. Nonetheless, outlier tests identified 3.6–6.6% of loci as high F_ST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high F_ST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small N_e in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.     

Genetic structure and demographic history of the endangered tree species Dysoxylum malabaricum (Meliaceae) in Western Ghats,India: implications for conservation in a biodiversity hotspot

The impact of fragmentation by human activities on genetic diversity of forest trees is an important concern in forest conservation, especially in tropical forests. Dysoxylum malabaricum (white cedar) is an economically important tree species, endemic to the Western Ghats, India, one of the world's eight most important biodiversity hotspots. As D. malabaricum is under pressure of disturbance and fragmentation together with overharvesting, conservation efforts are required in this species. In this study, range‐wide genetic structure of twelve D. malabaricum populations was evaluated to assess the impact of human activities on genetic diversity and infer the species’ evolutionary history, using both nuclear and chloroplast (cp) DNA simple sequence repeats (SSR). As genetic diversity and population structure did not differ among seedling, juvenile and adult age classes, reproductive success among the old‐growth trees and long distance seed dispersal by hornbills were suggested to contribute to maintain genetic diversity. The fixation index (F_IS) was significantly correlated with latitude, with a higher level of inbreeding in the northern populations, possibly reflecting a more severe ecosystem disturbance in those populations. Both nuclear and cpSSRs revealed northern and southern genetic groups with some discordance of their distributions; however, they did not correlate with any of the two geographic gaps known as genetic barriers to animals. Approximate Bayesian computation‐based inference from nuclear SSRs suggested that population divergence occurred before the last glacial maximum. Finally we discussed the implications of these results, in particular the presence of a clear pattern of historical genetic subdivision, on conservation policies.     

Using heterozygosity–fitness correlations to study inbreeding depression in an isolated population of white‐tailed deer founded by few individuals

A heterozygosity–fitness correlations (HFCs) may reflect inbreeding depression, but the extent to which they do so is debated. HFCs are particularly likely to occur after demographic disturbances such as population bottleneck or admixture. We here study HFC in an introduced and isolated ungulate population of white‐tailed deer Odocoileus virginianus in Finland founded in 1934 by four individuals. A total of 422 ≥ 1‐year‐old white‐tailed deer were collected in the 2012 hunting season in southern Finland and genotyped for 14 microsatellite loci. We find significant identity disequilibrium as estimated by g₂. Heterozygosity was positively associated with size‐ and age‐corrected body mass, but not with jaw size or (in males) antler score. Because of the relatively high identity disequilibrium, heterozygosity of the marker panel explained 51% of variation in inbreeding. Inbreeding explained approximately 4% of the variation in body mass and is thus a minor, although significant source of variation in body mass in this population. The study of HFC is attractive for game‐ and conservation‐oriented wildlife management because it presents an affordable and readily used approach for genetic monitoring that allowing identification of fitness costs associated with genetic substructuring in what may seem like a homogeneous population.     

Genetic diversity,admixture, and hatchery influence in Brook Trout (Salvelinus fontinalis) throughout western New York State

Although Brook Trout are distributed across most of eastern North America, population numbers have declined in many regions due to habitat loss, climate change, and competition with non‐native species. In New York State, Brook Trout habitat has been substantially reduced, with many areas showing complete extirpation of Brook Trout populations, predominantly in the western portion of the state. Small, fragmented populations are at risk of genetic diversity loss, inbreeding depression, and reduced fitness, leading to a greater potential for local extirpation. Genetic monitoring is a practical tool that can facilitate further conservation‐decision making regarding small populations. In this study, we used 12 microsatellite loci to examine 3,436 sampled Brook Trout, representing 75 sites from the Allegheny, Erie/Niagara, Genesee, Oswego, Lake Ontario, and Susquehanna drainage basins throughout western New York State. Three Brook Trout hatchery strains were also genetically characterized to evaluate the degree of hatchery introgression between wild populations and hatchery strains stocked in the region. Overall, estimates of genetic diversity varied widely: Allelic richness ranged from 2.23 to 7.485, and expected heterozygosity ranged from 0.402 to 0.766. As observed for Brook Trout in other regions, we found a high degree of genetic differentiation among populations, with all comparisons except one showing significant F_ST values. Hatchery introgression was found to be minimal, with estimates ranging from 1.96% to 3.10% of wild individuals exhibiting membership proportions to a hatchery strain cluster exceeding 10% (q ≥ 0.10). Results from this investigation can be used to prioritize management efforts for Brook Trout in western New York State and act as a baseline to monitor future population trends.