A comparison of SNP and STR loci for delineating population structure and performing individual genetic assignment

Influence of molecular marker type on estimating effective population size and other genetic parameters in a critically endangered parrot

Genetics is a fast-moving field, and for conservation practitioners or ecologists, it can be bewildering. The choice of marker used in studies is fundamental; in the literature, preference has recently shifted from microsatellites to single nucleotide polymorphism (SNP) loci. Understanding how marker type affects estimates of population genetic parameters is important in the context of conservation, especially because the accuracy of estimates has a bearing on the actions taken to protect threatened species. We compare parameter estimates between seven microsatellites, 3761 SNP loci, and a random subset of 100 SNPs for the exact same 324 individual swift parrots, Lathamus discolor, and also use 457 additional samples from subsequent years to compare SNP estimates. Both marker types estimated a lower H O than H E. We show that microsatellites and SNPs mainly indicate a lack of spatial genetic structure, except when a priori collection locations were used on the SNP data in a discriminant analysis of principal components (DAPC). The 100-SNP subset gave comparable results to when the full dataset was used. Estimates of effective population size (N e) were comparable between markers when the same individuals were considered, but SNPs had narrower confidence intervals. This is reassuring because conservation assessments that rely on population genetic estimates based on a few microsatellites are unlikely to be nullified by the general shift toward SNPs in the literature. However, estimates between markers and datasets varied considerably when only adult samples were considered; hence, including samples of all age groups is recommended to be used when available. The estimated N e was higher for the full SNP dataset (2010-2019) than the smaller comparison data (2010-2015), which might be a better reflection of the species status. The lower precision of microsatellites may not necessarily be a barrier for most conservation applications; however, SNPs will improve confidence limits, which may be useful for practitioners.

Comparative Analysis of Microsatellite and SNP Markers for Genetic Management of Red Deer

The analysis of population genetic structure and individual multilocus heterozygosity are crucial for wildlife management and conservation. Microsatellite markers have traditionally been used to assess these genetic parameters. However, single-nucleotide polymorphisms (SNPs) are becoming increasingly popular. Our goal here was to determine to what extent SNPs can provide better insights than microsatellites into the overall genetic status and population genetic processes in the species. To this end, we genotyped 210 red deer (Cervus elaphus) in the Spanish wild population with both 11 microsatellites and 31,712 SNPs. We compared parameters related to population genetic structure and individual multilocus heterozygosity obtained with both types of markers. Our results showed correlations between parameters measured using both microsatellites and SNPs, particularly those related to the level of genetic diversity and genetic differentiation. However, we found notably lower precision of microsatellites in measuring the distribution of genetic diversity among individuals. We conclude that microsatellites can be used to monitor the overall genetic status and detect broad patterns in red deer populations. Nevertheless, the greater precision of SNPs in inferring genetic structure and multilocus heterozygosity leads us to encourage scientists and wildlife managers to prioritize their use whenever possible.

Design and validation of a 63K genome-wide SNP-genotyping platform for caribou/reindeer (Rangifer tarandus)

BACKGROUND

Development of large single nucleotide polymorphism (SNP) arrays can make genomic data promptly available for conservation problematic. Medium and high-density panels can be designed with sufficient coverage to offer a genome-wide perspective and the generated genotypes can be used to assess different genetic metrics related to population structure, relatedness, or inbreeding. SNP genotyping could also permit sexing samples with unknown associated metadata as it is often the case when using non-invasive sampling methods favored for endangered species. Genome sequencing of wild species provides the necessary information to design such SNP arrays. We report here the development of a SNP-array for endangered Rangifer tarandus using a multi-platform sequencing approach from animals found in diverse populations representing the entire circumpolar distribution of the species.

RESULTS

From a very large comprehensive catalog of SNPs detected over the entire sample set (N = 894), a total of 63,336 SNPs were selected. SNP selection accounted for SNPs evenly distributed across the entire genome (~ every 50Kb) with known minor alleles across populations world-wide. In addition, a subset of SNPs was selected to represent rare and local alleles found in Eastern Canada which could be used for ecotype and population assignments - information urgently needed for conservation planning. In addition, heterozygosity from SNPs located in the X-chromosome and genotyping call-rate of SNPs located into the SRY gene of the Y-chromosome yielded an accurate and robust sexing assessment. All SNPs were validated using a high-throughput SNP-genotyping chip.

CONCLUSION

This design is now integrated into the first genome-wide commercially available genotyping platform for Rangifer tarandus. This platform would pave the way to future genomic investigation of populations for this endangered species, including estimation of genetic diversity parameters, population assignments, as well as animal sexing from genetic SNP data for non-invasive samples.

Consequences of introgression and gene flow on the genetic structure and diversity of Lima bean (Phaseolus lunatus L.) in its Mesoamerican diversity area

We evaluated the role of gene flow and wild-crop introgression on the structure and genetic diversity of Lima bean (Phaseolus lunatus) in the Yucatan Peninsula, an important Mesoamerican diversity area for this crop, using a genotyping-by-sequencing approach (15,168 SNP markers) and two scales. At the local scale, STRUCTURE and NGSEP analyses showed predominantly crop-to-wild introgression, but also evidence of a bidirectional gene flow in the two wild-weedy-crop complexes studied (Itzinté and Dzitnup). The ABBA-BABA tests showed a higher introgression in Itzinté (the older complex) than in Dzitnup (the younger one); at the allelic level, the wild-crop introgression in Itzinté was similar in both directions, in Dzitnup it was higher from crop-to-wild; and at the chromosomal level, introgression in Itzinté was from wild-to-crop, whereas in Dzitnup it occured in the opposite direction. Also, we found H _E values slightly higher in the domesticated accessions than in the wild ones, in both complexes (Itzinté: wild = 0.31, domesticated = 0.34; Dzinup: wild = 0.27, domesticated = 0.36), but %P and π estimators were higher in the wild accessions than in the domesticated ones. At a regional scale, STRUCTURE and MIGRATE showed a low gene flow, predominantly from crop-to-wild; and STRUCTURE, Neighbor-Joining and PCoA analyses indicated the existence of two wild groups and one domesticated group, with a marked genetic structure based in the existence of domesticated MI and wild MII gene pools. Also, at the regional scale, we found a higher genetic diversity in the wild accessions than in the domesticated ones, in all estimators used (e.g., H _E = 0.27 and H _E = 0.17, respectively). Our results indicate that gene flow and introgression are playing an important role at the local scale, but its consequences on the structure and genetic diversity of the Lima bean are not clearly reflected at the regional scale, where diversity patterns between wild and domesticated populations could be reflecting historical events.

Genomic Data Reveals Profound Genetic Structure and Multiple Glacial Refugia in Lonicera oblata (Caprifoliaceae), a Threatened Montane Shrub Endemic to North China

Characterizing genetic diversity and structure and identifying conservation units are both crucial for the conservation and management of threatened species. The development of high-throughput sequencing technology provides exciting opportunities for conservation genetics. Here, we employed the powerful SuperGBS method to identify 33, 758 high-quality single-nucleotide polymorphisms (SNP) from 134 individuals of a critically endangered montane shrub endemic to North China, Lonicera oblata. A low level of genetic diversity and a high degree of genetic differentiation among populations were observed based on the SNP data. Both principal component and phylogenetic analyses detected seven clusters, which correspond exactly to the seven geographic populations. Under the optimal K = 7, Admixture suggested the combination of the two small and geographically neighboring populations in the Taihang Mountains, Dongling Mountains, and Lijiazhuang, while the division of the big population of Jiankou Great Wall in the Yan Mountains into two clusters. High population genetic diversity and a large number of private alleles were detected in the four large populations, while low diversity and non-private alleles were observed for the remaining three small populations, implying the importance of these large populations as conservation units in priority. Demographic history inference suggested two drastic contractions of population size events that occurred after the Middle Pleistocene Transition and the Last Glacial Maximum, respectively. Combining our previous ecological niche modeling results with the present genomic data, there was a possible presence of glacial refugia in the Taihang and Yan Mountains, North China. This study provides valuable data for the conservation and management of L. oblata and broadens the understanding of the high biodiversity in the Taihang and Yan Mountains.

SNPs in Sheep: Characterization of Lithuanian Sheep Populations

In Lithuania, there are two recognised native sheep breeds: old native Lithuanian Coarsewooled and Lithuanian Blackface. In addition, in 2005, primitive Heidschnucke-type Skudde sheep were imported to Lithuania and were argued to possibly represent a lost Lithuanian sheep type. The aim of the study was to investigate the genetic variation in the two Lithuanian native sheep breeds, compare them with the imported Skudde sheep and establish the historical patterns of admixture and the genetic relatedness of Lithuanian sheep to British, Central European and Nordic sheep breeds included in the SheepHapMap study. In total, 72 individuals, representing two Lithuanian native and imported Skudde sheep breeds, were genotyped using a Neogen 12K Illumina Infinium chip. The population analysis was carried out by model-based clustering, principal component analysis and neighbour net analysis, and showed similar patterns for the Lithuanian sheep populations. Lithuanian Coarsewooled and Skudde in Lithuania have unique divergence and possibly some shared ancestry, while the Lithuanian Blackface conforms to a modern synthetic breed. The study clearly showed that the Coarsewooled and the Skudde breeds are distinct from each other. Historical data strongly suggest that the Coarsewooled breed represents a local breed, while the Skudde origin is less directly linked to the geographical area of modern-day Lithuania. Within the modern-day Lithuanian context, the Lithuanian Coarsewooled sheep is very important historical sheep type for conservation.

The details are in the genome-wide SNPs: Fine scale evolution of the Malaysian weedy rice

Weedy rice (Oryza spp.) is a major nuisance to rice farmers from all over the world. Although the emergence of weedy rice in East Malaysia on the island of Borneo is very recent, the threat to rice yield has reached an alarming stage. Using 47,027 genotyping-by-sequencing (GBS)-derived SNPs and candidate gene analysis of the plant architecture domestication gene TAC1, we assessed the genetic variations and evolutionary origin of weedy rice in East Malaysia. Our findings revealed two major evolutionary paths for genetically distinct weedy rice types. Whilst the cultivar-like weedy rice are very likely to be the weedy descendant of local coexisting cultivars, the wild-like weedy rice appeared to have arisen through two possible routes: (i) accidental introduction from Peninsular Malaysia weedy rice populations, and (ii) weedy descendants of coexisting cultivars. The outcome of our genetic analyses supports the notion that Sabah cultivars and Peninsular Malaysia weedy rice are the potential progenitors of Sabah weedy rice. Similar TAC1 haplotypes were shared between Malaysian cultivated and weedy rice populations, which further supported the findings of our GBS-SNP analyses. These different strains of weedy rice have convergently evolved shared traits, such as seeds shattering and open tillers. A comparison with our previous simple-sequence repeat-based population genetic analyses highlights the strength of genome-wide SNPs, including detection of admixtures and low-level introgression events. These findings could inform better strategic management for controlling the spread of weedy rice in the region.

Hybridization with mountain hares increases the functional allelic repertoire in brown hares

Brown hares (Lepus europaeus Pallas) are able to hybridize with mountain hares (L. timidus Linnaeus) and produce fertile offspring, which results in cross-species gene flow. However, not much is known about the functional significance of this genetic introgression. Using targeted sequencing of candidate loci combined with mtDNA genotyping, we found the ancestral genetic diversity in the Finnish brown hare to be small, likely due to founder effect and range expansion, while gene flow from mountain hares constitutes an important source of functional genetic variability. Some of this variability, such as the alleles of the mountain hare thermogenin (uncoupling protein 1, UCP1), might have adaptive advantage for brown hares, whereas immunity-related MHC alleles are reciprocally exchanged and maintained via balancing selection. Our study offers a rare example where an expanding species can increase its allelic variability through hybridization with a congeneric native species, offering a route to shortcut evolutionary adaptation to the local environmental conditions.

Reproductive development and genetic structure of the mycoheterotrophic orchid Pogoniopsis schenckii Cogn

BACKGROUND

Pogoniopsis schenckii Cogn. is a mycoheterotrophic orchid that can be used as a model to understand the influence of mycoheterotrophy at different stages of the reproductive cycle. We aimed to verify the presence of endophytic and epiphytic fungi at each stage of the reproductive process and investigated how the breeding system may relate to genetic structure and diversity of populations. In this study we performed anatomical and ultrastructural analyses of the reproductive organs, field tests to confirm the breeding system, and molecular analysis to assess genetic diversity and structure of populations.

RESULTS

During the development of the pollen grain, embryo sac and embryogenesis, no fungal infestation was observed. The presence of endophytic fungal hyphae was observed just within floral stems and indehiscent fruit. Beyond assuring the presence of fungus that promote seed germination, specific fungi hyphae in the fruit may affect other process, such as fruit ripening. As other mycoheterotrophic orchids, P. schenckii is autogamous, which may explain the low genetic diversity and high genetic structure in populations.

CONCLUSIONS

We discuss an interesting interaction: fungal hyphae in the indehiscent fruit. These fungal hyphae seem to play different roles inside fruit tissues, such as acting in the fruit maturation process and increasing the proximity between fungi and plant seeds even before dispersion occurs. As other mycoheterotrophic orchids, P. schenckii is autogamous, which may explain the low genetic diversity and high genetic structure in populations. Altogether, our findings provide important novel information about the mechanisms shaping ecology and evolution of fragmented populations of mycoheterotrophic plant.

Accuracy of Imputation of Microsatellite Markers from a 50K SNP Chip in Spanish Assaf Sheep

Simple Summary

Parentage misassignments directly affect genetic gain in traditional breeding programs. The use of genetic markers facilitates parentage verification. In sheep, microsatellite markers and single nucleotide polymorphism (SNP) markers have been proposed by the International Society of Animal Sciences (ISAG) for parentage testing. Since the implementation of genomic selection, the microsatellite information used for parental testing in previous generations is gradually being replaced by SNPs. However, parentage verifications should all be performed using the same technology. A strategy for transitioning from microsatellites to SNP markers, while avoiding extra genotyping costs, is the imputation of microsatellite alleles from SNP haplotypes. This study aims to identify the optimum approach, using a minimum number of SNPs to accurately impute microsatellite markers and developing a low-density SNP chip for parentage verification in the Assaf sheep breed. The imputation approach described here reached high accuracies using a low number of SNP markers, which supports the development of a low-density SNP chip that could avoid the problems of genotyping with both technologies, being a cost-effective method for parentage testing. This study will help sheep breeders to perform parentage verification when different genotyping platforms have been used across generations.

Abstract

Transitioning from traditional to new genotyping technologies requires the development of bridging methodologies to avoid extra genotyping costs. This study aims to identify the optimum number of single nucleotide polymorphisms (SNPs) necessary to accurately impute microsatellite markers to develop a low-density SNP chip for parentage verification in the Assaf sheep breed. The accuracy of microsatellite marker imputation was assessed with three metrics: genotype concordance (C), genotype dosage (length r²), and allelic dosage (allelic r²), for all imputation scenarios tested (0.5–10 Mb microsatellite flanking SNP windows). The imputation accuracy for the three metrics analyzed for all haplotype lengths tested was higher than 0.90 (C), 0.80 (length r²), and 0.75 (allelic r²), indicating strong genotype concordance. The window with 2 Mb length provides the best accuracy for the imputation procedure and the design of an affordable low-density SNP chip for parentage testing. We additionally evaluated imputation performance under two null models, naive (imputing the most common allele) and random (imputing by randomly selecting the allele), which in comparison showed weak genotype concordances (0.41 and 0.15, respectively). Therefore, we describe a precise methodology in the present article to impute multiallelic microsatellite genotypes from a low-density SNP chip in sheep and solve the problem of parentage verification when different genotyping platforms have been used across generations.

Unlocking the potential of a validated single nucleotide polymorphism array for genomic monitoring of trade in cheetahs (Acinonyx jubatus)

Cheetahs (Acinonyx jubatus) are listed as vulnerable on the International Union for Conservation of Nature Red List of Threatened Species. Threats include loss of habitat, human-wildlife conflict and illegal wildlife trade. In South Africa, the export of wild cheetah is a restricted activity under the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES), however, limited legal trade is permitted of animals born to captive parents. To effectively monitor the legal and illegal trade in South Africa, it was thus essential to develop a validated molecular test. Here, we designed a single nucleotide polymorphism (SNP) array for cheetah from Double Digest Restriction Associated DNA sequencing data for individual identification and parentage testing. In order to validate the array, unrelated individuals and 16 family groups consisting of both parents and one to three offspring were genotyped using the Applied Biosystems™ QuantStudio™ 12K Flex Real-Time PCR System. In addition, parentage assignments were compared to microsatellite data. Cross-species amplification was tested in various felids and cheetah sub-species in order to determine the utility of the SNP array in other species. We obtained successful genotyping results for 218 SNPs in cheetah (A. j. jubatus) with an optimal DNA input concentration ranging from 10 to 30 ng/µl. The combination of SNPs had a higher resolving power for individual identification compared to microsatellites and provided high assignment accuracy in known pedigrees. Cross-species amplification in other felids was determined to be limited. However, the SNP array demonstrated a clear genetic discrimination of two cheetah subspecies tested here. We conclude that the described SNP array is suitable for accurate parentage assignment and provides an important traceability tool for forensic investigations of cheetah trade.

Comparative population genetic analyses suggest hybrid origin of R hododendron pubicostatum, an endangered plant species with extremely small populations endemic to Yunnan, China

Gene flow between sympatric congeneric plants is thought to be very common and may pose serious threats to endangered species. In the present study, we evaluate the genetic diversity and divergence of three sympatric Rhododendron species in Jiaozi Mountain using newly developed microsatellites through the Illumina MiSeq sequencing approach. Genetic diversity of all three Rhododendron species studied was moderate in comparison to genetic parameters previously reported from species of this genus. Interestingly, genetic structure analysis of the three species identified a possible hybrid origin of the threatened Rh. pubicostatum. This sympatry should be considered a unimodal hybrid zone, since R h. pubicostatum is predominant here. Unimodal hybrid zones are uncommon in Rhododendron, despite the fact that hybridization frequently occurs in the genus. Issues pertaining to the conservation of R h. pubicostatum resulting from admixture of genetic material from its parental species are discussed.

Variation in breeding practices and geographic isolation drive subpopulation differentiation, contributing to the loss of genetic diversity within dog breed lineages

Background

Discrete breed ideals are not restricted to delimiting dog breeds from another, but also are key drivers of subpopulation differentiation. As genetic differentiation due to population fragmentation results in increased rates of inbreeding and loss of genetic diversity, detecting and alleviating the reasons of population fragmentation can provide effective tools for the maintenance of healthy dog breeds.

Results

Using a genome-wide SNP array, we detected genetic differentiation to subpopulations in six breeds, Belgian Shepherd, English Greyhound, Finnish Lapphund, Italian Greyhound, Labrador Retriever and Shetland Sheepdog, either due to geographical isolation or as a result of differential breeding strategies. The subpopulation differentiation was strongest in show dog lineages.

Conclusions

Besides geographical differentiation caused by founder effect and lack of gene flow, selection on champion looks or restricted pedigrees is a strong driver of population fragmentation. Artificial barriers for gene flow between the different subpopulations should be recognized, their necessity evaluated critically and perhaps abolished in order to maintain genetic diversity within a breed. Subpopulation differentiation might also result in false positive signals in genome-wide association studies of different traits.

Lay summary

Purebred dogs are, by definition, reproductively isolated from other breeds. However, similar isolation can also occur within a breed due to conflicting breeder ideals and geographic distances between the dog populations. We show here that both of these examples can contribute to breed division, with subsequent loss of genetic variation in the resulting breed lineages. Breeders should avoid creating unnecessary boundaries between breed lineages and facilitate the exchange of dogs between countries.

Chromosome aberrations in pressure-induced triploid Atlantic salmon

BACKGROUND

Triploid organisms have three sets of chromosomes. In Atlantic salmon, hydrostatic pressure treatment of newly fertilized eggs has been extensively used to produce triploids which are functionally sterile due to their unpaired chromosomes. These fish often perform poorly on commercial farms, sometimes without explanation. Inheritance patterns in individuals subjected to pressure treatment have not been investigated in Atlantic salmon thus far. However, work on other species suggests that this treatment can result in aberrant inheritance. We therefore studied this in Atlantic salmon by genotyping 16 polymorphic microsatellites in eyed eggs and juveniles which had been subjected to pressure-induction of triploidy. Communally reared juveniles including fish subjected to pressure-induction of triploidy and their diploid siblings were included as a control.

RESULTS

No diploid offspring were detected in any of the eggs or juveniles which were subjected to hydrostatic pressure; therefore, the induction of triploidy was highly successful. Aberrant inheritance was nevertheless observed in 0.9% of the eggs and 0.9% of the juveniles that had been subjected to pressure treatment. In the communally reared fish, 0.3% of the fish subjected to pressure treatment displayed aberrant inheritance, while their diploid controls displayed 0% aberrant inheritance. Inheritance errors included two eyed eggs lacking maternal DNA across all microsatellites, and, examples in both eggs and juveniles of either the maternal or paternal allele lacking in one of the microsatellites. All individuals displaying chromosome aberrations were otherwise triploid.

CONCLUSIONS

This is the first study to document aberrant inheritance in Atlantic salmon that have been subjected to pressure-induction of triploidy. Our experiments unequivocally demonstrate that even when induction of triploidy is highly successful, this treatment can cause chromosome aberrations in this species. Based upon our novel data, and earlier studies in other organisms, we hypothesize that in batches of Atlantic salmon where low to modest triploid induction rates have been reported, aberrant inheritance is likely to be higher than the rates observed here. Therefore, we tentatively suggest that this could contribute to the unexplained poor performance of triploid salmon that is occasionally reported in commercial aquaculture. These hypotheses require further investigation.

Inference of Breed Structure in Farm Animals: Empirical Comparison between SNP and Microsatellite Performance

Knowledge of population structure is essential to improve the management and conservation of farm animal genetic resources. Microsatellites, which have long been popular for this type of analysis, are more and more neglected in favor of whole-genome single nucleotide polymorphism (SNP) chips that are now available for the main farmed animal species. In this study, we compared genetic patterns derived from microsatellites to that inferred by SNPs, considering three pairs of datasets of sheep and cattle. Population genetic differentiation analyses (Fixation index, F_ST), as well as STRUCTURE analyses showed a very strong consistency between the two types of markers. Microsatellites gave pictures that were largely concordant with SNPs, although less accurate. The best concordance was found in the most complex dataset, which included 17 French sheep breeds (with a Pearson correlation coefficient of 0.95 considering the 136 values of pairwise F_ST, obtained with both types of markers). The use of microsatellites reduces the cost and the related analyses do not require specific computer equipment (i.e., information technology (IT) infrastructure able to provide adequate computing and storage capacity). Therefore, this tool may still be a very appropriate solution to evaluate, in a first stage, the general state of livestock at national scales. At a time when local breeds are disappearing at an alarming rate, it is urgent to improve our knowledge of them, in particular by promoting tools accessible to the greatest number.

Ancestry-informative marker (AIM) SNP panel for the Malay population

Ancestry-informative markers (AIMs) can be used to infer the ancestry of an individual to minimize the inaccuracy of self-reported ethnicity in biomedical research. In this study, we describe three methods for selecting AIM SNPs for the Malay population (Malay AIM panel) using different approaches based on pairwise F_ST, informativeness for assignment (I_n), and PCA-correlated SNPs (PCAIMs). These Malay AIM panels were extracted from genotype data stored in SNP arrays hosted by the Malaysian node of the Human Variome Project (MyHVP) and the Singapore Genome Variation Project (SGVP). In particular, genotype data from a total of 165 Malay individuals were analyzed, comprising data on 117 individual genotypes from the Affymetrix SNP-6 SNP array platform and data on 48 individual genotypes from the OMNI 2.5 Illumina SNP array platform. The HapMap phase 3 database (1397 individuals from 11 populations) was used as a reference for comparison with the Malay genotype data. The accuracy of each resulting Malay AIM panel was evaluated using a machine learning "ancestry-predictive model" constructed by using WEKA, a comprehensive machine learning platform written in Java. A total of 1250 SNPs were finally selected, which successfully identified Malay individuals from other world populations with an accuracy of 90%, but the accuracy decreased to 80% using 157 SNPs according to the pairwise F_ST method, while a panel of 200 SNPs selected using I_n and PCAIMs could be used to identify Malay individuals with an accuracy of approximately 80%.

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

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

An extensive common-garden study with domesticated and wild Atlantic salmon in the wild reveals impact on smolt production and shifts in fitness traits

Interactions between domesticated escapees and wild conspecifics represent a threat to the genetic integrity and fitness of native populations. For Atlantic salmon, the recurrent presence of large numbers of domesticated escapees in the wild makes it necessary to better understand their impacts on native populations. We planted 254,400 eggs from 75 families of domesticated, F1-hybrid, and wild salmon in a river containing up- and downstream traps. Additionally, 41,630 hatchery smolts of the same pedigrees were released into the river. Over 8 years, 6,669 out-migrating smolts and 356 returning adults were recaptured and identified to their families of origin with DNA. In comparison with wild salmon, domesticated fish had substantially lower egg to smolt survival (1.8% vs. 3.8% across cohorts), they migrated earlier in the year (11.8 days earlier across years), but they only displayed marginally larger smolt sizes and marginally lower smolt ages. Upon return to freshwater, domesticated salmon were substantially larger at age than wild salmon (2.4 vs. 2.0, 4.8 vs. 3.2, and 8.5 vs. 5.6 kg across sexes for 1, 2, and 3 sea-winter fish) and displayed substantially lower released smolt to adult survival (0.41% vs. 0.94% across releases). Overall, egg-to-returning adult survival ratios were 1:0.76:0.30 and 1:0.44:0.21 for wild:F1-hybrid:domesticated salmon, respectively, using two different types of data. This study represents the most updated and extensive analysis of domesticated, hybrid, and wild salmon in the wild and provides the first documentation of a clear genetic difference in the timing of smolt migration-an adaptive trait presumed to be linked with optimal timing of entry to seawater. We conclude that spawning and hybridization of domesticated escapees can lead to (i) reduced wild smolt output and therefore wild adult abundance, through resource competition in freshwater, (ii) reduced total adult abundance due to freshwater competition and reduced marine survival of domesticated salmon, and (iii) maladaptive changes in phenotypic traits.

Landscape genomics of an obligate mutualism: Concordant and discordant population structures between the leafcutter ant Atta texana and its two main fungal symbiont types

To explore landscape genomics at the range limit of an obligate mutualism, we use genotyping-by-sequencing (ddRADseq) to quantify population structure and the effect of host-symbiont interactions between the northernmost fungus-farming leafcutter ant Atta texana and its two main types of cultivated fungus. Genome-wide differentiation between ants associated with either of the two fungal types is of the same order of magnitude as differentiation associated with temperature and precipitation across the ant's entire range, suggesting that specific ant-fungus genome-genome combinations may have been favoured by selection. For the ant hosts, we found a broad cline of genetic structure across the range, and a reduction of genetic diversity along the axis of range expansion towards the range margin. This population-genetic structure was concordant between the ants and one cultivar type (M-fungi, concordant clines) but discordant for the other cultivar type (T-fungi). Discordance in population-genetic structures between ant hosts and a fungal symbiont is surprising because the ant farmers codisperse with their vertically transmitted fungal symbionts. Discordance implies that (a) the fungi disperse also through between-nest horizontal transfer or other unknown mechanisms, and (b) genetic drift and gene flow can differ in magnitude between each partner and between different ant-fungus combinations. Together, these findings imply that variation in the strength of drift and gene flow experienced by each mutualistic partner affects adaptation to environmental stress at the range margin, and genome-genome interactions between host and symbiont influence adaptive genetic differentiation of the host during range evolution in this obligate mutualism.

Environmental extremes drive population structure at the northern range limit of Atlantic salmon in North America

Conservation of exploited species requires an understanding of both genetic diversity and the dominant structuring forces, particularly near range limits, where climatic variation can drive rapid expansions or contractions of geographic range. Here, we examine population structure and landscape associations in Atlantic salmon (Salmo salar) across a heterogeneous landscape near the northern range limit in Labrador, Canada. Analysis of two amplicon-based data sets containing 101 microsatellites and 376 single nucleotide polymorphisms (SNPs) from 35 locations revealed clear differentiation between populations spawning in rivers flowing into a large marine embayment (Lake Melville) compared to coastal populations. The mechanisms influencing the differentiation of embayment populations were investigated using both multivariate and machine-learning landscape genetic approaches. We identified temperature as the strongest correlate with genetic structure, particularly warm temperature extremes and wider annual temperature ranges. The genomic basis of this divergence was further explored using a subset of locations (n = 17) and a 220K SNP array. SNPs associated with spatial structuring and temperature mapped to a diverse set of genes and molecular pathways, including regulation of gene expression, immune response, and cell development and differentiation. The results spanning molecular marker types and both novel and established methods clearly show climate-associated, fine-scale population structure across an environmental gradient in Atlantic salmon near its range limit in North America, highlighting valuable approaches for predicting population responses to climate change and managing species sustainability.

Targeted resequencing of coding DNA sequences for SNP discovery in nonmodel species

Targeted capture coupled with high-throughput sequencing can be used to gain information about nuclear sequence variation at hundreds to thousands of loci. Divergent reference capture makes use of molecular data of one species to enrich target loci in other (related) species. This is particularly valuable for nonmodel organisms, for which often no a priori knowledge exists regarding these loci. Here, we have used targeted capture to obtain data for 809 nuclear coding DNA sequences (CDS) in a nonmodel organism, the Eurasian lynx Lynx lynx, using baits designed with the help of the published genome of a related model organism (the domestic cat Felis catus). Using this approach, we were able to survey intraspecific variation at hundreds of nuclear loci in L. lynx across the species' European range. A large set of biallelic candidate SNPs was then evaluated using a high-throughput SNP genotyping platform (Fluidigm), which we then reduced to a final 96 SNP-panel based on assay performance and reliability; validation was carried out with 100 additional Eurasian lynx samples not included in the SNP discovery phase. The 96 SNP-panel developed from CDS performed very successfully in the identification of individuals and in population genetic structure inference (including the assignment of individuals to their source population). In keeping with recent studies, our results show that genic SNPs can be valuable for genetic monitoring of wildlife species.

Fresh is best: Accurate SNP genotyping from koala scats

Maintaining genetic diversity is a crucial component in conserving threatened species. For the iconic Australian koala, there is little genetic information on wild populations that is not either skewed by biased sampling methods (e.g., sampling effort skewed toward urban areas) or of limited usefulness due to low numbers of microsatellites used. The ability to genotype DNA extracted from koala scats using next‐generation sequencing technology will not only help resolve location sample bias but also improve the accuracy and scope of genetic analyses (e.g., neutral vs. adaptive genetic diversity, inbreeding, and effective population size). Here, we present the successful SNP genotyping (1272 SNP loci) of koala DNA extracted from scat, using a proprietary DArTseq^™ protocol. We compare genotype results from two‐day‐old scat DNA and 14‐day‐old scat DNA to a blood DNA template, to test accuracy of scat genotyping. We find that DNA from fresher scat results in fewer loci with missing information than DNA from older scat; however, 14‐day‐old scat can still provide useful genetic information, depending on the research question. We also find that a subset of 209 conserved loci can accurately identify individual koalas, even from older scat samples. In addition, we find that DNA sequences identified from scat samples through the DArTseq^™ process can provide genetic identification of koala diet species, bacterial and viral pathogens, and parasitic organisms.

Genetic Diversity and Population Structure of Ethiopian Sheep Populations Revealed by High-Density SNP Markers

Sheep in Ethiopia are adapted to a wide range of environments, including extreme habitats. Elucidating their genetic diversity is critical for improving breeding strategies and mapping quantitative trait loci associated with productivity. To this end, the present study investigated the genetic diversity and population structure of five Ethiopian sheep populations exhibiting distinct phenotypes and sampled from distinct production environments, including arid lowlands and highlands. To investigate the genetic relationships in greater detail and infer population structure of Ethiopian sheep breeds at the continental and global levels, we analyzed genotypic data of selected sheep breeds from the Ovine SNP50K HapMap dataset. All Ethiopian sheep samples were genotyped with Ovine Infinium HD SNP BeadChip (600K). Mean genetic diversity ranged from 0.29 in Arsi-Bale to 0.32 in Menz sheep, while estimates of genetic differentiation among populations ranged from 0.02 to 0.07, indicating low to moderate differentiation. An analysis of molecular variance revealed that 94.62 and 5.38% of the genetic variation was attributable to differences within and among populations, respectively. Our population structure analysis revealed clustering of five Ethiopian sheep populations according to tail phenotype and geographic origin-i.e., short fat-tailed (very cool high-altitude), long fat-tailed (mid to high-altitude), and fat-rumped (arid low-altitude), with clear evidence of admixture between long fat-tailed populations. North African sheep breeds showed higher levels of within-breed diversity, but were less differentiated than breeds from Eastern and Southern Africa. When African breeds were grouped according to geographic origin (North, South, and East), statistically significant differences were detected among groups (regions). A comparison of population structure between Ethiopian and global sheep breeds showed that fat-tailed breeds from Eastern and Southern Africa clustered together, suggesting that these breeds were introduced to the African continent via the Horn and migrated further south.

Global genotype flow in Cercospora beticola populations confirmed through genotyping-by-sequencing

Genotyping-by-sequencing (GBS) was conducted on 333 Cercospora isolates collected from Beta vulgaris (sugar beet, table beet and swiss chard) in the USA and Europe. Cercospora beticola was confirmed as the species predominantly isolated from leaves with Cercospora leaf spot (CLS) symptoms. However, C. cf. flagellaris also was detected at a frequency of 3% in two table beet fields in New York. Resolution of the spatial structure and identification of clonal lineages in C. beticola populations using genome-wide single nucleotide polymorphisms (SNPs) obtained from GBS was compared to genotyping using microsatellites. Varying distance thresholds (bitwise distance = 0, 1.854599 × 10-4, and 1.298 × 10-3) were used for delineation of clonal lineages in C. beticola populations. Results supported previous reports of long distance dispersal of C. beticola through genotype flow. The GBS-SNP data set provided higher resolution in discriminating clonal lineages; however, genotype identification was impacted by filtering parameters and the distance threshold at which the multi-locus genotypes (MLGs) were contracted to multi-locus lineages. The type of marker or different filtering strategies did not impact estimates of population differentiation and structure. Results emphasize the importance of robust filtering strategies and designation of distance thresholds for delineating clonal lineages in population genomics analyses that depend on individual assignment and identification of clonal lineages. Detection of recurrent clonal lineages shared between the USA and Europe, even in the relaxed-filtered SNP data set and with a conservative distance threshold for contraction of MLGs, provided strong evidence for global genotype flow in C. beticola populations. The implications of intercontinental migration in C. beticola populations for CLS management are discussed.

Conservation priorities for endangered Indian tigers through a genomic lens

Tigers have lost 93% of their historical range worldwide. India plays a vital role in the conservation of tigers since nearly 60% of all wild tigers are currently found here. However, as protected areas are small (<300 km² on average), with only a few individuals in each, many of them may not be independently viable. It is thus important to identify and conserve genetically connected populations, as well as to maintain connectivity within them. We collected samples from wild tigers (Panthera tigris tigris) across India and used genome-wide SNPs to infer genetic connectivity. We genotyped 10,184 SNPs from 38 individuals across 17 protected areas and identified three genetically distinct clusters (corresponding to northwest, southern and central India). The northwest cluster was isolated with low variation and high relatedness. The geographically large central cluster included tigers from central, northeastern and northern India, and had the highest variation. Most genetic diversity (62%) was shared among clusters, while unique variation was highest in the central cluster (8.5%) and lowest in the northwestern one (2%). We did not detect signatures of differential selection or local adaptation. We highlight that the northwest population requires conservation attention to ensure persistence of these tigers.

SNP discovery and characterisation in White Rhino (Ceratotherium simum) with application to parentage assignment

The white rhino is one of the great success stories of modern wildlife conservation, growing from as few as 50-100 animals in the 1880s, to approximately 20,000 white rhinoceros remaining today. However, illegal trade in conservational rhinoceros horns is adding constant pressure on remaining populations. Captive management of ex situ populations of endangered species using molecular methods can contribute to improving the management of the species. Here we compare for the first time the utility of 33 Single Nucleotide Polymorphisms (SNPs) and nine microsatellites (MS) in isolation and in combination for assigning parentage in captive White Rhinoceros. We found that a combined dataset of SNPs and microsatellites was most informative with the highest confidence level. This study thus provided us with a useful set of SNP and MS markers for parentage and relatedness testing. Further assessment of the utility of these markers over multiple (> three) generations and the incorporation of a larger variety of relationships among individuals (e.g. half-siblings or cousins) is strongly suggested.

Estimating genomic diversity and population differentiation - an empirical comparison of microsatellite and SNP variation in Arabidopsis halleri

BACKGROUND

Microsatellite markers are widely used for estimating genetic diversity within and differentiation among populations. However, it has rarely been tested whether such estimates are useful proxies for genome-wide patterns of variation and differentiation. Here, we compared microsatellite variation with genome-wide single nucleotide polymorphisms (SNPs) to assess and quantify potential marker-specific biases and derive recommendations for future studies. Overall, we genotyped 180 Arabidopsis halleri individuals from nine populations using 20 microsatellite markers. Twelve of these markers were originally developed for Arabidopsis thaliana (cross-species markers) and eight for A. halleri (species-specific markers). We further characterized 2 million SNPs across the genome with a pooled whole-genome re-sequencing approach (Pool-Seq).

RESULTS

Our analyses revealed that estimates of genetic diversity and differentiation derived from cross-species and species-specific microsatellites differed substantially and that expected microsatellite heterozygosity (SSR-H e) was not significantly correlated with genome-wide SNP diversity estimates (SNP-H e and θ Watterson) in A. halleri. Instead, microsatellite allelic richness (A r) was a better proxy for genome-wide SNP diversity. Estimates of genetic differentiation among populations (F ST) based on both marker types were correlated, but microsatellite-based estimates were significantly larger than those from SNPs. Possible causes include the limited number of microsatellite markers used, marker ascertainment bias, as well as the high variance in microsatellite-derived estimates. In contrast, genome-wide SNP data provided unbiased estimates of genetic diversity independent of whether genome- or only exome-wide SNPs were used. Further, we inferred that a few thousand random SNPs are sufficient to reliably estimate genome-wide diversity and to distinguish among populations differing in genetic variation.

CONCLUSIONS

We recommend that future analyses of genetic diversity within and differentiation among populations use randomly selected high-throughput sequencing-based SNP data to draw conclusions on genome-wide diversity patterns. In species comparable to A. halleri, a few thousand SNPs are sufficient to achieve this goal.

SNPs selected by information content outperform randomly selected microsatellite loci for delineating genetic identification and introgression in the endangered dark European honeybee (Apis mellifera mellifera)

The honeybee (Apis mellifera) has been threatened by multiple factors including pests and pathogens, pesticides and loss of locally adapted gene complexes due to replacement and introgression. In western Europe, the genetic integrity of the native A. m. mellifera (M-lineage) is endangered due to trading and intensive queen breeding with commercial subspecies of eastern European ancestry (C-lineage). Effective conservation actions require reliable molecular tools to identify pure-bred A. m. mellifera colonies. Microsatellites have been preferred for identification of A. m. mellifera stocks across conservation centres. However, owing to high throughput, easy transferability between laboratories and low genotyping error, SNPs promise to become popular. Here, we compared the resolving power of a widely utilized microsatellite set to detect structure and introgression with that of different sets that combine a variable number of SNPs selected for their information content and genomic proximity to the microsatellite loci. Contrary to every SNP data set, microsatellites did not discriminate between the two lineages in the PCA space. Mean introgression proportions were identical across the two marker types, although at the individual level, microsatellites' performance was relatively poor at the upper range of Q-values, a result reflected by their lower precision. Our results suggest that SNPs are more accurate and powerful than microsatellites for identification of A. m. mellifera colonies, especially when they are selected by information content.

Next Gen Pop Gen: implementing a high-throughput approach to population genetics in boarfish (Capros aper)

The recently developed approach for microsatellite genotyping by sequencing (GBS) using individual combinatorial barcoding was further improved and used to assess the genetic population structure of boarfish (Capros aper) across the species' range. Microsatellite loci were developed de novo and genotyped by next-generation sequencing. Genetic analyses of the samples indicated that boarfish can be subdivided into at least seven biological units (populations) across the species' range. Furthermore, the recent apparent increase in abundance in the northeast Atlantic is better explained by demographic changes within this area than by influx from southern or insular populations. This study clearly shows that the microsatellite GBS approach is a generic, cost-effective, rapid and powerful method suitable for full-scale population genetic studies-a crucial element for assessment, sustainable management and conservation of valuable biological resources.

Accuracy of Assignment of Atlantic Salmon (Salmo salar L.) to Rivers and Regions in Scotland and Northeast England Based on Single Nucleotide Polymorphism (SNP) Markers

Understanding the habitat use patterns of migratory fish, such as Atlantic salmon (Salmo salar L.), and the natural and anthropogenic impacts on them, is aided by the ability to identify individuals to their stock of origin. Presented here are the results of an analysis of informative single nucleotide polymorphic (SNP) markers for detecting genetic structuring in Atlantic salmon in Scotland and NE England and their ability to allow accurate genetic stock identification. 3,787 fish from 147 sites covering 27 rivers were screened at 5,568 SNP markers. In order to identify a cost-effective subset of SNPs, they were ranked according to their ability to differentiate between fish from different rivers. A panel of 288 SNPs was used to examine both individual assignments and mixed stock fisheries and eighteen assignment units were defined. The results improved greatly on previously available methods and, for the first time, fish caught in the marine environment can be confidently assigned to geographically coherent units within Scotland and NE England, including individual rivers. As such, this SNP panel has the potential to aid understanding of the various influences acting upon Atlantic salmon on their marine migrations, be they natural environmental variations and/or anthropogenic impacts, such as mixed stock fisheries and interactions with marine power generation installations.

Population genetic inferences using immune gene SNPs mirror patterns inferred by microsatellites

Single nucleotide polymorphisms (SNPs) are replacing microsatellites for population genetic analyses, but it is not apparent how many SNPs are needed or how well SNPs correlate with microsatellites. We used data from the gopher tortoise, Gopherus polyphemus-a species with small populations, to compare SNPs and microsatellites to estimate population genetic parameters. Specifically, we compared one SNP data set (16 tortoises from four populations sequenced at 17 901 SNPs) to two microsatellite data sets, a full data set of 101 tortoises and a partial data set of 16 tortoises previously genotyped at 10 microsatellites. For the full microsatellite data set, observed heterozygosity, expected heterozygosity and F_ST were correlated between SNPs and microsatellites; however, allelic richness was not. The same was true for the partial microsatellite data set, except that allelic richness, but not observed heterozygosity, was correlated. The number of clusters estimated by structure differed for each data set (SNPs = 2; partial microsatellite = 3; full microsatellite = 4). Principle component analyses (PCA) showed four clusters for all data sets. More than 800 SNPs were needed to correlate with allelic richness, observed heterozygosity and expected heterozygosity, but only 100 were needed for F_ST . The number of SNPs typically obtained from next-generation sequencing (NGS) far exceeds the number needed to correlate with microsatellite parameter estimates. Our study illustrates that diversity, F_ST and PCA results from microsatellites can mirror those obtained with SNPs. These results may be generally applicable to small populations, a defining feature of endangered and threatened species, because theory predicts that genetic drift will tend to outweigh selection in small populations.

Population structure in Atlantic cod in the eastern North Sea-Skagerrak-Kattegat: early life stage dispersal and adult migration

Background

In marine fish species, where pelagic egg and larvae drift with ocean currents, population structure has been suggested to be maintained by larval retention due to hydrographic structuring and by homing of adult fish to natal areas. Whilst natal homing of adults has been demonstrated for anadromous and coral reef fishes, there are few documented examples of philopatric migration in temperate marine fish species.

Results

Here, we demonstrate temporally stable genetic differentiation among spawning populations of Atlantic cod (Gadus morhua L.), and present genetic and behavioural evidence for larval drift and philopatric migration in the eastern North Sea-Skagerrak-Kattegat area. We show that juvenile cod collected in the eastern Skagerrak and central Kattegat are genetically similar to cod from offshore spawning areas in the eastern North Sea. Genetic assignment of individual 2–5 year old fish indicates that cod residing at, or migrating towards, spawning areas in Kattegat and the North Sea display philopatric behaviours.

Conclusions

Together these findings suggest a loop between spawning, larval drift and adult return-migrations to spawning areas and underlines that both oceanographic processes and migratory behaviour in the adult phase may be important for stock separation and integrity in marine temperate fishes such as Atlantic cod.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-016-1878-9) contains supplementary material, which is available to authorized users.

A common garden design reveals population-specific variability in potential impacts of hybridization between populations of farmed and wild Atlantic salmon, Salmo salar L

Released individuals can have negative impacts on native populations through various mechanisms, including competition, disease transfer and introduction of maladapted gene complexes. Previous studies indicate that the level of farmed Atlantic salmon introgression in native populations is population specific. However, few studies have explored the potential role of population diversity or river characteristics, such as temperature, on the consequences of hybridization. We compared freshwater growth of multiple families derived from two farmed, five wild and two F1 hybrid salmon populations at three contrasting temperatures (7°C, 12°C and 16°C) in a common garden experiment. As expected, farmed salmon outgrew wild salmon at all temperatures, with hybrids displaying intermediate growth. However, differences in growth were population specific and some wild populations performed better than others relative to the hybrid and farmed populations at certain temperatures. Therefore, the competitive balance between farmed and wild salmon may depend both on the thermal profile of the river and on the genetic characteristics of the respective farmed and wild strains. While limited to F1 hybridization, this study shows the merits in adopting a more complex spatially resolved approach to risk management of local populations.

Effectiveness of a 95 SNP panel for the screening of breed label fraud in the Chinese meat market

Breed assignment has proved to be useful to control meat trade and protect the value of special productions. Meat-related frauds have been detected in China; therefore, 95 SNPs selected from the ISAG core panel were evaluated to develop an automated and technologically updated tool to screen breed label fraud in the Chinese meat market. A total of 271 animals from four Chinese yellow cattle (CYC) populations, six Bos taurus breeds, two Bos indicus and one composite were used. The allocation test distinguished European, Japanese and Zebu breeds, and two Chinese genetic components. It correctly allocated Japanese Black, Zebu and British breeds in 100, 90 and 89% of samples, respectively. CYC evidenced the Zebu, Holstein and Limousin introgression. The test did not detect CYC components in any of the 25 samples from Argentinean butchers. The method could be useful to certify Angus, Hereford and Japanese Black meat, but a modification in the panel would be needed to differentiate other breeds.

Single nucleotide polymorphism-based dispersal estimates using noninvasive sampling

Quantifying dispersal within wild populations is an important but challenging task. Here we present a method to estimate contemporary, individual-based dispersal distance from noninvasively collected samples using a specialized panel of 96 SNPs (single nucleotide polymorphisms). One main issue in conducting dispersal studies is the requirement for a high sampling resolution at a geographic scale appropriate for capturing the majority of dispersal events. In this study, fecal samples of brown bear (Ursus arctos) were collected by volunteer citizens, resulting in a high sampling resolution spanning over 45,000 km² in Gävleborg and Dalarna counties in Sweden. SNP genotypes were obtained for unique individuals sampled (n = 433) and subsequently used to reconstruct pedigrees. A Mantel test for isolation by distance suggests that the sampling scale was appropriate for females but not for males, which are known to disperse long distances. Euclidean distance was estimated between mother and offspring pairs identified through the reconstructed pedigrees. The mean dispersal distance was 12.9 km (SE 3.2) and 33.8 km (SE 6.8) for females and males, respectively. These results were significantly different (Wilcoxon’s rank-sum test: P-value = 0.02) and are in agreement with the previously identified pattern of male-biased dispersal. Our results illustrate the potential of using a combination of noninvasively collected samples at high resolution and specialized SNPs for pedigree-based dispersal models.

Toward a genome-wide approach for detecting hybrids: informative SNPs to detect introgression between domestic cats and European wildcats (Felis silvestris)

Endemic gene pools have been severely endangered by human-mediated hybridization, which is posing new challenges in the conservation of several vertebrate species. The endangered European wildcat is an example of this problem, as several natural populations are suffering introgression of genes from the domestic cat. The implementation of molecular methods for detecting hybridization is crucial for supporting appropriate conservation programs on the wildcat. In this study, genetic variation at 158 single-nucleotide polymorphisms (SNPs) was analyzed in 139 domestic cats, 130 putative European wildcats and 5 captive-bred hybrids (N=274). These SNPs were variable both in wild (HE=0.107) and domestic cats (HE=0.340). Although we did not find any SNP that was private in any population, 22 SNPs were monomorphic in wildcats and pairwise FCT values revealed marked differences between domestic and wildcats, with the most divergent 35 loci providing an average FCT>0.74. The power of all the loci to accurately identify admixture events and discriminate the different hybrid categories was evaluated. Results from simulated and real genotypes show that the 158 SNPs provide successful estimates of admixture, with 100% hybrid individuals (two to three generations in the past) being correctly identified in STRUCTURE and over 92% using the NEWHYBRIDS' algorithm. None of the unclassified cats were wrongly allocated to another hybrid class. Thirty-five SNPs, showing the highest FCT values, provided the most parsimonious panel for robust inferences of parental and first generations of admixed ancestries. This approach may be used to further reconstruct the evolution of wildcat populations and, hopefully, to develop sound conservation guidelines for its legal protection in Europe.

Discovery and characterization of single nucleotide polymorphisms in coho salmon, Oncorhynchus kisutch

Molecular population genetic analyses have become an integral part of ecological investigation and population monitoring for conservation and management. Microsatellites have been the molecular marker of choice for such applications over the last several decades, but single nucleotide polymorphism (SNP) markers are rapidly expanding beyond model organisms. Coho salmon (Oncorhynchus kisutch) is native to the north Pacific Ocean and its tributaries, where it is the focus of intensive fishery and conservation activities. As it is an anadromous species, coho salmon typically migrate across multiple jurisdictional boundaries, complicating management and requiring shared data collection methods. Here, we describe the discovery and validation of a suite of novel SNPs and associated genotyping assays which can be used in the genetic analyses of this species. These assays include 91 that are polymorphic in the species and one that discriminates it from a sister species, Chinook salmon. We demonstrate the utility of these SNPs for population assignment and phylogeographic analyses, and map them against the draft trout genome. The markers constitute a large majority of all SNP markers described for coho salmon and will enable both population- and pedigree-based analyses across the southern part of the species native range.

The frequency of spontaneous triploidy in farmed Atlantic salmon produced in Norway during the period 2007-2014

BACKGROUND

Spontaneous triploidy has been reported in a number of fish species, and is often linked with in vivo or in vitro ageing of eggs post ovulation. Here, we provide the first investigation into the frequency of spontaneous triploidy in farmed Atlantic salmon by analysing more than 4000 fish from 55 farms, and approximately 1000 recaptured escapees, all sampled in the period 2007-2014. In addition, we compare microsatellite genotyping against flow cytometry and red blood cell diameter in a set of 45 putatively diploid and 45 putatively triploid Atlantic salmon.

RESULTS

The three methods implemented for ploidy determination gave consistent results, thus validating the methods used here. Overall, 2.0% spontaneous triploids were observed in salmon sampled on farms. The frequency of spontaneous triploids varied greatly among sea cages (0-28%), but they were observed in similar frequencies among the three primary breeding companies (1.8-2.4%). Spontaneous triploids were observed in all farming regions in Norway, and in all years sampled. Spontaneous triploids were also observed among the escapees recaptured in both the marine environment and in rivers.

CONCLUSIONS

Spontaneous triploidy in commercially produced Atlantic salmon is likely to be a result of the practices employed by the industry. For logistical reasons, there is sometimes a pause of hours, and in some cases overnight, between killing the female broodfish, removal of her eggs, and fertilization. This gives the eggs time to age post ovulation, and increases the probability of duplication of the maternal chromosome set by inhibition of the second polar body release after normal meiosis II in the oocyte.

Challenges in analysis and interpretation of microsatellite data for population genetic studies

Advancing technologies have facilitated the ever-widening application of genetic markers such as microsatellites into new systems and research questions in biology. In light of the data and experience accumulated from several years of using microsatellites, we present here a literature review that synthesizes the limitations of microsatellites in population genetic studies. With a focus on population structure, we review the widely used fixation (F ST) statistics and Bayesian clustering algorithms and find that the former can be confusing and problematic for microsatellites and that the latter may be confounded by complex population models and lack power in certain cases. Clustering, multivariate analyses, and diversity-based statistics are increasingly being applied to infer population structure, but in some instances these methods lack formalization with microsatellites. Migration-specific methods perform well only under narrow constraints. We also examine the use of microsatellites for inferring effective population size, changes in population size, and deeper demographic history, and find that these methods are untested and/or highly context-dependent. Overall, each method possesses important weaknesses for use with microsatellites, and there are significant constraints on inferences commonly made using microsatellite markers in the areas of population structure, admixture, and effective population size. To ameliorate and better understand these constraints, researchers are encouraged to analyze simulated datasets both prior to and following data collection and analysis, the latter of which is formalized within the approximate Bayesian computation framework. We also examine trends in the literature and show that microsatellites continue to be widely used, especially in non-human subject areas. This review assists with study design and molecular marker selection, facilitates sound interpretation of microsatellite data while fostering respect for their practical limitations, and identifies lessons that could be applied toward emerging markers and high-throughput technologies in population genetics.

A low-density SNP array for analyzing differential selection in freshwater and marine populations of threespine stickleback (Gasterosteus aculeatus)

Background

The threespine stickleback (Gasterosteus aculeatus) has become an important model species for studying both contemporary and parallel evolution. In particular, differential adaptation to freshwater and marine environments has led to high differentiation between freshwater and marine stickleback populations at the phenotypic trait of lateral plate morphology and the underlying candidate gene Ectodysplacin (EDA). Many studies have focused on this trait and candidate gene, although other genes involved in marine-freshwater adaptation may be equally important. In order to develop a resource for rapid and cost efficient analysis of genetic divergence between freshwater and marine sticklebacks, we generated a low-density SNP (Single Nucleotide Polymorphism) array encompassing markers of chromosome regions under putative directional selection, along with neutral markers for background.

Results

RAD (Restriction site Associated DNA) sequencing of sixty individuals representing two freshwater and one marine population led to the identification of 33,993 SNP markers. Ninety-six of these were chosen for the low-density SNP array, among which 70 represented SNPs under putatively directional selection in freshwater vs. marine environments, whereas 26 SNPs were assumed to be neutral. Annotation of these regions revealed several genes that are candidates for affecting stickleback phenotypic variation, some of which have been observed in previous studies whereas others are new.

Conclusions

We have developed a cost-efficient low-density SNP array that allows for rapid screening of polymorphisms in threespine stickleback. The array provides a valuable tool for analyzing adaptive divergence between freshwater and marine stickleback populations beyond the well-established candidate gene Ectodysplacin (EDA).

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-867) contains supplementary material, which is available to authorized users.