An empirical comparison of population genetic analyses using microsatellite and SNP data for a species of conservation concern

Genomic insights into the conservation status of the Idle Crayfish Austropotamobius bihariensis Pârvulescu, 2019: low genetic diversity in the endemic crayfish species of the Apuseni Mountains

BACKGROUND

Biodiversity in freshwater ecosystems is declining due to an increased anthropogenic footprint. Freshwater crayfish are keystone species in freshwater ecosystems and play a crucial role in shaping the structure and function of their habitats. The Idle Crayfish Austropotamobius bihariensis is a native European species with a narrow distribution range, endemic to the Apuseni Mountains (Romania). Although its area is small, the populations are anthropogenically fragmented. In this context, the assessment of its conservation status is timely.

RESULTS

Using a reduced representation sequencing approach, we identified 4875 genomic SNPs from individuals belonging to 13 populations across the species distribution range. Subsequent population genomic analyses highlighted low heterozygosity levels, low number of private alleles and small effective population size. Our structuring analyses revealed that the genomic similarity of the populations is conserved within the river basins.

CONCLUSION

Genomic SNPs represented excellent tools to gain insights into intraspecific genomic diversity and population structure of the Idle Crayfish. Our study highlighted that the analysed populations are at risk due to their limited genetic diversity, which makes them extremely vulnerable to environmental alterations. Thus, our results emphasize the need for conservation measures and can be used as a baseline to establish species management programs.

Thirty years of conservation breeding: Assessing the genetic diversity of captive Livingstone's fruit bats

Fruit bats (genus Pteropus) are typically island-endemic species important in seed dispersal and reforestation that are vulnerable to increased extinction risk. An effective method of reducing extinction risk in vulnerable species that cannot be conserved in their native habitat is establishing an ex-situ captive breeding programme. Due to anthropogenic threats and low population numbers, in the early 1990s, a captive breeding programme was established at Jersey Zoo, British Isles, for Critically Endangered Livingstone's fruit bats (Pteropus livingstonii). Here we use six polymorphic microsatellite loci to assess genetic diversity in the captive breeding population of Livingstone's fruit bats (P. livingstonii), 30 years after the programme's establishment, investigating change over generations and comparing our findings with published data from the wild population. We found no significant difference between the genetic diversity in the captive and wild populations of Livingstone's fruit bats (P. livingstonii), in both expected heterozygosity and allelic richness. The captive population has retained a comparable level of genetic diversity to that documented in the wild, and there has been no significant decline in genetic diversity over the last 30 years. We advise that a full pedigree of the paternal lineage is created to improve the management of the captive breeding programme and further reduce the possibility of inbreeding. However, it appears that the captive breeding programme is currently effective at maintaining genetic diversity at levels comparable to those seen in the wild population, which suggests reintroductions could be viable if genetic diversity remains stable in captivity.

Fine-scale adaptive divergence and population genetic structure of Aedes aegypti in Metropolitan Manila, Philippines

BACKGROUND

The adaptive divergence of Aedes aegypti populations to heterogeneous environments can be a driving force behind the recent expansion of their habitat distribution and outbreaks of dengue disease in urbanized areas. In this study, we investigated the population genomics of Ae. aegypti at a regional scale in Metropolitan Manila, Philippines.

METHODS

We used the Pool-Seq double digestion restriction-site association DNA sequencing (ddRAD-Seq) approach to generate a high number of single nucleotide polymorphisms (SNPs), with the aim to determine local adaptation and compare the population structure with 11 microsatellite markers. A total of 217 Ae. aegypti individuals from seven female and seven male populations collected from Metropolitan Manila were used in the assays.

RESULTS

We detected 65,473 SNPs across the populations, of which 76 were non-neutral SNPs. Of these non-neutral SNPs, the multivariate regression test associated 50 with eight landscape variables (e.g. open space, forest, etc.) and 29 with five climate variables (e.g. air temperature, humidity, etc.) (P-value range 0.005-0.045) in female and male populations separately. Male and female populations exhibited contrasting spatial divergence, with males exhibiting greater divergence than females, most likely reflecting the different dispersal abilities of male and female mosquitoes. In the comparative analysis of the same Ae. aegypti individuals, the pairwise FST values of 11 microsatellite markers were lower than those of the neutral SNPs, indicating that the neutral SNPs generated via pool ddRAD-Seq were more sensitive in terms of detecting genetic differences between populations at fine-spatial scales.

CONCLUSIONS

Overall, our study demonstrates the utility of pool ddRAD-Seq for examining genetic differences in Ae. aegypti populations in areas at fine-spatial scales that could inform vector control programs such as Wolbachia-infected mosquito mass-release programs. This in turn would provide information on mosquito population dispersal patterns and the potential barriers to mosquito movement within and around the release area. In addition, the potential of environmental adaptability observed in Ae. aegypti could help population control efforts.

Evaluation of DNA yield from various tissue and sampling sources for use in single nucleotide polymorphism panels

Genetics studies are used by wildlife managers and researchers to gain inference into a population of a species of interest. To gain these insights, microsatellites have been the primary method; however, there currently is a shift from microsatellites to single nucleotide polymorphisms (SNPs). With the different DNA requirements between microsatellites and SNPs, an investigation into which samples can provide adequate DNA yield is warranted. Using samples that were collected from previous genetic projects from regions in the USA from 2014 to 2021, we investigated the DNA yield of eight sample categories to gain insights into which provided adequate DNA to be used in ddRADseq or already developed high- or medium-density SNP panels. We found seven sample categories that met the DNA requirements for use in all three panels, and one sample category that did not meet any of the three panels requirements; however, DNA integrity was highly variable and not all sample categories that met panel DNA requirements could be considered high quality DNA. Additionally, we used linear random-effects models to determine which covariates would have the greatest influence on DNA yield. We determined that all covariates (tissue type, storage method, preservative, DNA quality, time until DNA extraction and time after DNA extraction) could influence DNA yield.

Effects of landscape, resource use, and body size on genetic structure in bee populations

Quantifying genetic structure and levels of genetic variation are fundamentally important to predicting the ability of populations to persist in human-altered landscapes and adapt to future environmental changes. Genetic structure reflects the dispersal of individuals over generations, which can be mediated by species-level traits or environmental factors. Dispersal distances are commonly positively associated with body size and negatively associated with the amount of degraded habitat between sites, motivating the investigation of these potential drivers of dispersal concomitantly. We quantified genetic structure and genetic variability within populations of seven bee species from the genus Euglossa across fragmented landscapes. We genotyped bees at SNP loci and tested the following predictions: (1) deforested areas restrict gene flow; (2) larger species have lower genetic structure; (3) species with greater resource specialization have higher genetic structure; and (4) sites surrounded by more intact habitat have higher genetic diversity. Contrasting with previous work on bees, we found no associations between body size and genetic structure. Genetic structure was higher for species with greater resource specialization, and the amount of intact habitat between or surrounding sites was positively associated with parameters reflecting gene flow and genetic diversity. These results challenge the dominant paradigm that individuals of larger species disperse farther, and they suggest that landscape and resource requirements are important factors mediating dispersal.

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.

Ex situ conservation of two rare oak species using microsatellite and SNP markers

Plant collections held by botanic gardens and arboreta are key components of ex situ conservation. Maintaining genetic diversity in such collections allows them to be used as resources for supplementing wild populations. However, most recommended minimum sample sizes for sufficient ex situ genetic diversity are based on microsatellite markers, and it remains unknown whether these sample sizes remain valid in light of more recently developed next-generation sequencing (NGS) approaches. To address this knowledge gap, we examine how ex situ conservation status and sampling recommendations differ when derived from microsatellites and single nucleotide polymorphisms (SNPs) in garden and wild samples of two threatened oak species. For Quercus acerifolia, SNPs show lower ex situ representation of wild allelic diversity and slightly lower minimum sample size estimates than microsatellites, while results for each marker are largely similar for Q. boyntonii. The application of missing data filters tends to lead to higher ex situ representation, while the impact of different SNP calling approaches is dependent on the species being analyzed. Measures of population differentiation within species are broadly similar between markers, but larger numbers of SNP loci allow for greater resolution of population structure and clearer assignment of ex situ individuals to wild source populations. Our results offer guidance for future ex situ conservation assessments utilizing SNP data, such as the application of missing data filters and the usage of a reference genome, and illustrate that both microsatellites and SNPs remain viable options for botanic gardens and arboreta seeking to ensure the genetic diversity of their collections.

Analysis of Genetic Diversity in Romanian Carpatina Goats Using SNP Genotyping Data

Animal husbandry is one of man's oldest occupations. It began with the domestication of animals and developed continuously, in parallel with the evolution of human society. The selection and improvement of goats in Romania was not a clearly defined objective until around 1980. In recent years, with the increasing economic value given to goats, breeding programs are becoming established. In Romania, a few goat genetic studies using microsatellites and mtDNA have been carried out; however, a systematic characterization of the country's goat genomic resources remains missing. In this study, we analyzed the genetic variability of Carpatina goats from four distinct geographical areas (northern, north-eastern, eastern and southern Romania), using the Illumina OvineSNP60 (RefSeq ARS1) high-density chip for 67 goats. Heterozygosity values, inbreeding coefficients and effective population size across all autosomes were calculated for those populations that inhabit high- and low-altitude and high- and low-temperature environments. Diversity, as measured by expected heterozygosity (HE), ranged from 0.413 in the group from a low-temperature environment to 0.420 in the group from a high-temperature environment. Within studied groups, the HT (high temperature) goats were the only group with a positive but low average inbreeding coefficient value, which was 0.009. After quality control (QC) analysis, 46,965 SNPs remained for analysis (MAF < 0.01). LD was calculated for each chromosome separately. The Ne has been declining since the time of domestication, having recently reached 123, 125, 185 and 92 for the HA (high altitude), LA (low altitude), HT (high temperature) and LT (low temperature) group, respectively. Our study revealed a low impact of inbreeding in the Carpatina population, and the Ne trend also indicated a steep decline in the last hundred years. These results will contribute to the genetic improvement of the Carpatina breed.

Genomic insights into isolation of the threatened Florida crested caracara (Caracara plancus)

We conducted a population genomic study of the crested caracara (Caracara plancus) using samples (n = 290) collected from individuals in Florida, Texas, and Arizona, United States. Crested caracaras are non-migratory raptors ranging from the southern tip of South America to the southern United States, including a federally protected relict population in Florida long thought to have been isolated since the last ice age. Our objectives were to evaluate genetic diversity and population structure of Florida's apparently isolated population and to evaluate taxonomic relationships of crested caracaras at the northern edge of their range. Using DNA purified from blood samples, we conducted double-digest restriction site associated DNA sequencing and sequenced the mitochondrial ND2 gene. Analyses of population structure using over 9,000 SNPs suggest that two major clusters are best supported, one cluster including only Florida individuals and the other cluster including Arizona and Texas individuals. Both SNPs and mitochondrial haplotypes reveal the Florida population to be highly differentiated genetically from Arizona and Texas populations, whereas, Arizona and Texas populations are moderately differentiated from each other. The Florida population's mitochondrial haplotypes form a separate monophyletic group, while Arizona and Texas populations share mitochondrial haplotypes. Results of this study provide substantial genetic evidence that Florida's crested caracaras have experienced long-term isolation from caracaras in Arizona and Texas and thus, represent a distinct evolutionary lineage possibly warranting distinction as an Evolutionarily Significant Unit (ESU) or subspecies. This study will inform conservation strategies focused on long-term survival of Florida's distinct, panmictic population.

Development and evaluation of a novel single nucleotide polymorphism panel for North American bison

Genome-wide single nucleotide polymorphism (SNP) genotyping platforms have become increasingly popular in characterizing livestock and wildlife populations, replacing traditional methods such as microsatellite fragment analysis. Herein, we report the development and evaluation of a novel bison SNP panel for population management and conservation. Initially, 2474 autosomal SNPs were selected from existing bison whole-genome sequences and variable sites among bison on the GGSP bovine 50K Chip, based on minor allele frequency, data completeness, and chromosome location. Additionally, 20 mitochondrial SNPs were chosen to identify known mitochondrial haplotypes in bison according to previous research. The SNPs were further evaluated using genotyping-by-sequencing with 190 bison, representing the historical lineages that survived the major population crash of the late 1800s. Variants with high potential for genotyping error were filtered out, and the remaining SNPs were placed on a custom Illumina™ array. The final panel consisting of 798 autosomal and 13 mitochondrial SNPs was used to establish baseline genetic parameters, compare populations, and assign mitochondrial haplotypes in 995 bison across ten populations. These SNPs were also found to be highly informative for individual animal identification and parentage assignment. This SNP panel provides a powerful new method to establish a baseline for estimating genetic health of bison populations and a new tool for bison managers to make informed management decisions based on genetic information specific to their populations.

Ppe.RPT/SSC-1: from QTL mapping to a predictive KASP test for ripening time and soluble solids concentration in peach

Genomic regions associated with ripening time (RPT) and soluble solids concentration (SSC) were mapped using a pedigreed population including multiple F1 and F2 families from the Clemson University peach breeding program (CUPBP). RPT and SSC QTLs were consistently identified in two seasons (2011 and 2012) and the average datasets (average of two seasons). A target region spanning 10,981,971-11,298,736 bp on chromosome 4 of peach reference genome used for haplotype analysis revealed four haplotypes with significant differences in trait values among different diplotype combinations. Favorable alleles at the target region for both RPT and SSC were determined and a DNA test for predicting RPT and SSC was developed. Two Kompetitive Allele Specific PCR (KASP) assays were validated on 84 peach cultivars and 163 seedlings from the CUPBP, with only one assay (Ppe.RPT/SSC-1) needed to predict between early and late-season ripening cultivars and low and high SSC. These results advance our understanding of the genetic basis of RPT and SSC and facilitate selection of new peach cultivars with the desired RPT and SSC.

Contrasting patterns of sequence variation in steelhead populations reflect distinct evolutionary processes

Multiple evolutionary processes influence genome-wide allele frequencies and quantifying effects of genetic drift, and multiple forms of selection remain challenging in natural populations. Here, we investigate variation at major effect loci in contrast to patterns of neutral drift across a wide collection of steelhead (Oncorhynchus mykiss) populations that have declined in abundance due to anthropogenic impacts. Whole-genome resequencing of 74 populations of steelhead revealed genome-wide patterns (~8 million SNPs) consistent with expected neutral population structure. However, allelic variation at major effect loci associated with adult migration timing (chromosome 28: GREB1L/ROCK1) and age at maturity (chromosome 25: SIX6) reflected how selection has acted on phenotypic variation in contrast with neutral structure. Variation at major effect loci was influenced by evolutionary processes with differing signals between the strongly divergent Coastal and Inland lineages, while allele frequencies within and among populations within the Inland lineage have been driven by local natural selection as well as recent anthropogenic influences. Recent anthropogenic effects appeared to have influenced the frequency of major effect alleles including artificial selection for specific traits in hatchery stocks with subsequent gene flow into natural populations. Selection from environmental factors at various scales has also likely influenced variation for major effect alleles. These results reveal evolutionary mechanisms that influence allele frequencies at major effect loci that are critical for conservation of phenotypic traits and life history variation of this protected species.

Morphological and genetic mechanisms underlying the plasticity of the coral Porites astreoides across depths in Bermuda

The widespread decline of shallow-water coral reefs has fueled interest in assessing whether mesophotic reefs can act as refugia replenishing deteriorated shallower reefs through larval exchange. Here we explore the morphological and molecular basis facilitating survival of planulae and adults of the coral Porites astreoides (Lamarck, 1816; Hexacorallia: Poritidae) along the vertical depth gradient in Bermuda. We found differences in micro-skeletal features such as bigger calyxes and coarser surface of the skeletal spines in shallow corals. Yet, tomographic reconstructions reveal an analogous mineral distribution between shallow and mesophotic adults, pointing to similar skeleton growth dynamics. Our study reveals patterns of host genetic connectivity and minimal symbiont depth-zonation across a broader depth range than previously known for this species in Bermuda. Transcriptional variations across life stages showed different regulation of metabolism and stress response functions, unraveling molecular responses to environmental conditions at different depths. Overall, these findings increase our understanding of coral acclimatory capability across broad vertical gradients, ultimately allowing better evaluation of the refugia potential of mesophotic reefs.

Genomics Reveals Complex Population History and Unexpected Diversity of Eurasian Otters (Lutra lutra) in Britain Relative to Genetic Methods

Conservation genetic analyses of many endangered species have been based on genotyping of microsatellite loci and sequencing of short fragments of mtDNA. The increase in power and resolution afforded by whole genome approaches may challenge conclusions made on limited numbers of loci and maternally inherited haploid markers. Here, we provide a matched comparison of whole genome sequencing versus microsatellite and control region (CR) genotyping for Eurasian otters (Lutra lutra). Previous work identified four genetically differentiated "stronghold" populations of otter in Britain, derived from regional populations that survived the population crash of the 1950s-1980s. Using whole genome resequencing data from 45 samples from across the British stronghold populations, we confirmed some aspects of population structure derived from previous marker-driven studies. Importantly, we showed that genomic signals of the population crash bottlenecks matched evidence from otter population surveys. Unexpectedly, two strongly divergent mitochondrial lineages were identified that were undetectable using CR fragments, and otters in the east of England were genetically distinct and surprisingly variable. We hypothesize that this previously unsuspected variability may derive from past releases of Eurasian otters from other, non-British source populations in England around the time of the population bottleneck. Our work highlights that even reasonably well-studied species may harbor genetic surprises, if studied using modern high-throughput sequencing methods.

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.

Genome-wide variant analyses reveal new patterns of admixture and population structure in Australian dingoes

Admixture between species is a cause for concern in wildlife management. Canids are particularly vulnerable to interspecific hybridisation, and genetic admixture has shaped their evolutionary history. Microsatellite DNA testing, relying on a small number of genetic markers and geographically restricted reference populations, has identified extensive domestic dog admixture in Australian dingoes and driven conservation management policy. But there exists a concern that geographic variation in dingo genotypes could confound ancestry analyses that use a small number of genetic markers. Here, we apply genome-wide single-nucleotide polymorphism (SNP) genotyping to a set of 402 wild and captive dingoes collected from across Australia and then carry out comparisons to domestic dogs. We then perform ancestry modelling and biogeographic analyses to characterise population structure in dingoes and investigate the extent of admixture between dingoes and dogs in different regions of the continent. We show that there are at least five distinct dingo populations across Australia. We observed limited evidence of dog admixture in wild dingoes. Our work challenges previous reports regarding the occurrence and extent of dog admixture in dingoes, as our ancestry analyses show that previous assessments severely overestimate the degree of domestic dog admixture in dingo populations, particularly in south-eastern Australia. These findings strongly support the use of genome-wide SNP genotyping as a refined method for wildlife managers and policymakers to assess and inform dingo management policy and legislation moving forwards.

The role and risks of selective adaptation in extreme coral habitats

The alarming rate of climate change demands new management strategies to protect coral reefs. Environments such as mangrove lagoons, characterized by extreme variations in multiple abiotic factors, are viewed as potential sources of stress-tolerant corals for strategies such as assisted evolution and coral propagation. However, biological trade-offs for adaptation to such extremes are poorly known. Here, we investigate the reef-building coral Porites lutea thriving in both mangrove and reef sites and show that stress-tolerance comes with compromises in genetic and energetic mechanisms and skeletal characteristics. We observe reduced genetic diversity and gene expression variability in mangrove corals, a disadvantage under future harsher selective pressure. We find reduced density, thickness and higher porosity in coral skeletons from mangroves, symptoms of metabolic energy redirection to stress response functions. These findings demonstrate the need for caution when utilizing stress-tolerant corals in human interventions, as current survival in extremes may compromise future competitive fitness.

Multiple Introductions and Distinct Genetic Groups of Canada Goldenrod (Solidago canadensis) in China Revealed by Genomic Single-Nucleotide Polymorphisms

Despite numerous studies reported in the context of ecology, the introduction history of the infamous invasive plant Canada goldenrod (Solidago canadensis L.) remains elusive. In the present study, we explored the sources and the number of introduction events of this species from its native areas into China. Using the genotyping-by-sequencing approach, we identified 34,035 selectively neutral single-nucleotide polymorphism (SNP) markers to infer the evolutionary trajectories of 77 S. canadensis individuals. Both the principal component analysis and the ADMIXTURE analysis revealed two genetic groups that are sympatric to each other in China and suggested the absence of genetic admixtures. The phylogenetic analysis indicated three feasible introduction routes and multiple introduction events of Canada goldenrod into China. Specifically, the one from the USA directly into China, the other from the USA into China through Japan, and the third from the USA into China through Europe. Based on the site frequency spectrum of these identified SNPs, we inferred strong bottleneck events for both genetic groups, and that the multiple introductions did not rescue the decline of genetic diversity. To conclude, multiple introduction events, genetic bottlenecks, and potential human-mediated spread characterize the introduction history of Canada goldenrod in China. The present study harnesses the power of SNP data in deciphering the evolutionary trajectory of invasive plants and paves the way for future studies concerning the invasion mechanism of Canada goldenrod.

A genetic warning system for a hierarchically structured wildlife monitoring framework

Genetic variation is a well-known indicator of population fitness yet is not typically included in monitoring programs for sensitive species. Additionally, most programs monitor populations at one scale, which can lead to potential mismatches with ecological processes critical to species' conservation. Recently developed methods generating hierarchically nested population units (i.e., clusters of varying scales) for greater sage-grouse (Centrocercus urophasianus) have identified population trend declines across spatiotemporal scales to help managers target areas for conservation. The same clusters used as a proxy for spatial scale can alert managers to local units (i.e., neighborhood-scale) with low genetic diversity, further facilitating identification of management targets. We developed a genetic warning system utilizing previously developed hierarchical population units to identify management-relevant areas with low genetic diversity within the greater sage-grouse range. Within this warning system we characterized conservation concern thresholds based on values of genetic diversity and developed a statistical model for microsatellite data to robustly estimate these values for hierarchically nested populations. We found that 41 of 224 neighborhood-scale clusters had low genetic diversity, 23 of which were coupled with documented local population trend decline. We also found evidence of cross-scale low genetic diversity in the small and isolated Washington population, unlikely to be reversed through typical local management actions alone. The combination of low genetic diversity and a declining population suggests relatively high conservation concern. Our findings could further facilitate conservation action prioritization in combination with population trend assessments and (or) local information, and act as a base-line of genetic diversity for future comparison. Importantly, the approach we used is broadly applicable across taxa.

Monitoring changes in the genetic structure of Brown Tsaiya duck selected for feeding efficiency by microsatellite markers

OBJECTIVE

Few studies have genetically monitored chickens over time, and no research has been conducted on ducks. To ensure the sustainable management of key duck breeds, we used microsatellite markers to monitor Brown Tsaiya ducks over time genetically.

METHODS

The second, fourth, sixth to eighth generations of the Brown Tsaiya duck selected for feeding efficiency and control lines were included in this study to investigate the genetic variations, effective population size, population structure and the differentiation between populations over time with 11 microsatellite markers derived from Brown Tsaiya duck.

RESULTS

The results showed there were a slight decrease in the genetic variations and an increase in within-population inbreeding coefficient (FIS) in both lines, but no consistent increase in FIS was observed in each line. The effective population size in the second and eighth generations was 27.2 for the selected line and 23.9 for the control line. The change in allele richness showed a downward trend over time, and the selected line was slightly lower than the control line in each generation. The number of private alleles (Np) in the selected line were higher than in the control line. Moderate differentiation was observed between the second and eighth generations in the selected line (FST = 0.0510) and the control line (FST = 0.0606). Overall, differentiation tended to increase with each generation, but genetic variation and structure did not change considerably after six generations in the two lines.

CONCLUSION

This study provides a reference for poultry conservation and helps to implement cross-generation genetic monitoring and breeding plans in other duck breeds or lines to promote sustainable management.

Polyandry and non-random fertilisation maintain long-term genetic diversity in an isolated island population of adders (Vipera berus)

Conservation genetic theory suggests that small and isolated populations should be subject to reduced genetic diversity i.e., heterozygosity and allelic diversity. Our 34 years study of an isolated island population of adders (Vipera berus) in southern Sweden challenges this notion. Despite a lack of gene flow and a yearly mean estimated reproductive adult population size of only 65 adult adders (range 12-171), the population maintains high levels of heterozygosity and allelic diversity similar to that observed in two mainland populations. Even a 14-year major "bottleneck" i.e., a reduction in adult adder numbers, encompassing at least four adder generations, did not result in any reduction in the island adders' heterozygosity and allelic diversity. Female adders are polyandrous, and fertilisation is non-random, which our empirical data and modelling suggest are underpinning the maintenance of the population's high level of heterozygosity. Our empirical results and subsequent modelling suggest that the positive genetic effects of polyandry in combination with non-random fertilisation, often overlooked in conservation genetic analyses, deserve greater consideration when predicting long-term survival of small and isolated populations.

Inferring genetic structure when there is little: population genetics versus genomics of the threatened bat Miniopterus schreibersii across Europe

Despite their paramount importance in molecular ecology and conservation, genetic diversity and structure remain challenging to quantify with traditional genotyping methods. Next-generation sequencing holds great promises, but this has not been properly tested in highly mobile species. In this article, we compared microsatellite and RAD-sequencing (RAD-seq) analyses to investigate population structure in the declining bent-winged bat (Miniopterus schreibersii) across Europe. Both markers retrieved general patterns of weak range-wide differentiation, little sex-biased dispersal, and strong isolation by distance that associated with significant genetic structure between the three Mediterranean Peninsulas, which could have acted as glacial refugia. Microsatellites proved uninformative in individual-based analyses, but the resolution offered by genomic SNPs illuminated on regional substructures within several countries, with colonies sharing migrators of distinct ancestry without admixture. This finding is consistent with a marked philopatry and spatial partitioning between mating and rearing grounds in the species, which was suspected from marked-recaptured data. Our study advocates that genomic data are necessary to properly unveil the genetic footprints left by biogeographic processes and social organization in long-distant flyers, which are otherwise rapidly blurred by their high levels of gene flow.

Genetic diversity and connectivity of moose (Alces americanus americanus) in eastern North America

Genetic diversity is critical to a population’s ability to overcome gradual environment change. Large-bodied wildlife existing in regions with relatively high human population density are vulnerable to isolation-induced genetic drift, population bottlenecks, and loss of genetic diversity. Moose (Alces americanus americanus) in eastern North America have a complex history of drastic population changes. Current and potential threats to moose populations in this region could be exacerbated by loss of genetic diversity and connectivity among subpopulations. Existing genetic diversity, gene flow, and population clustering and fragmentation of eastern North American moose are not well quantified, while physical and anthropogenic barriers to population connectivity already exist. Here, single nucleotide polymorphism (SNP) genotyping of 507 moose spanning five northeastern U.S. states and one southeastern Canadian province indicated low diversity, with a high proportion of the genomes sharing identity-by-state, with no consistent evidence of non-random mating. Gene flow estimates indicated bidirectionality between all pairs of sampled areas, with magnitudes reflecting clustering and differentiation patterns. A Discriminant Analysis of Principal Components analysis indicated that these genotypic data were best described with four clusters and indicated connectivity across the Saint Lawrence River and Seaway, a potential physical barrier to gene flow. Tests for genetic differentiation indicated restricted gene flow between populations across the Saint Lawrence River and Seaway, and between many sampled areas facing expanding human activity. These results document current genetic variation and connectivity of moose populations in eastern North America, highlight potential challenges to current population connectivity, and identify areas for future research and conservation.

History, Phylogeny, Biodiversity, and New Computer-Based Tools for Efficient Micropropagation and Conservation of Pistachio (Pistacia spp.) Germplasm

The word "pstk" [pistag], used in the ancient Persian language, is the linguistic root from which the current name "pistachio", used worldwide, derives. The word pistachio is generally used to designate the plants and fruits of a single species: Pistacia vera L. Both the plant and its fruits have been used by mankind for thousands of years, specifically the consumption of its fruits by Neanderthals has been dated to about 300,000 years ago. Native to southern Central Asia (including northern Afghanistan and northeastern Iran), its domestication and cultivation occurred about 3000 years ago in this region, spreading to the rest of the Mediterranean basin during the Middle Ages and finally being exported to America and Australia at the end of the 19th century. The edible pistachio is an excellent source of unsaturated fatty acids, carbohydrates, proteins, dietary fiber, vitamins, minerals and bioactive phenolic compounds that help promote human health through their antioxidant capacity and biological activities. The distribution and genetic diversity of wild and domesticated pistachios have been declining due to increasing population pressure and climatic changes, which have destroyed natural pistachio habitats, and the monoculture of selected cultivars. As a result, the current world pistachio industry relies mainly on a very small number of commercial cultivars and rootstocks. In this review we discuss and summarize the current status of: etymology, origin, domestication, taxonomy and phylogeny by molecular analysis (RAPID, RFLP, AFLP, SSR, ISSR, IRAP, eSSR), main characteristics and world production, germplasm biodiversity, main cultivars and rootstocks, current conservation strategies of both conventional propagation (seeds, cutting, and grafting), and non-conventional propagation methods (cryopreservation, slow growth storage, synthetic seed techniques and micropropagation) and the application of computational tools (Design of Experiments (DoE) and Machine Learning: Artificial Neural Networks, Fuzzy logic and Genetic Algorithms) to design efficient micropropagation protocols for the genus Pistacia.

Genetic diversity of two native sheep breeds by genome-wide analysis of single nucleotide polymorphisms

Wallachian and Sumava sheep are autochthonous breeds that have undergone a significant bottleneck effect and subsequent restoration efforts. The first objective of this study was to evaluate the degree of genetic variability of both breeds and, therefore, the current management of the breeding. The second was to determine whether these two breeds still retain their genetic uniqueness in relation to each other and other breeds, despite regenerative interventions. Our data consisted of 48 individuals of Sumava and 37 individuals of Wallachian sheep. The comparison data contained 25 other breeds (primarily European) from the HapMap dataset generated by the International Sheep Genomics Consortium. When comparing all 27 breeds, the Czech breeds clustered with 15 other breeds and formed a single branch with them according to Nei's distances. At the same time, however, the clusters of both breeds were integral and easily distinguishable from the others when displayed with principal component analysis (PCA). Population substructure analysis did not show any common genetic ancestry of the Czech national breeds and breeds used for regeneration or, eventually, breeds whose ancestral population was used for regeneration. The average values of F_ST were higher in Wallachian sheep (F_ST = 0.14) than in Sumava sheep (F_ST = 0.08). The linkage disequilibrium (LD) extension per autosome was higher in Wallachian than in Sumava sheep. Consequently, the N_e estimates five generations ago were 68 for Sumava versus 34 for Wallachian sheep. Both native Czech breeds exhibit a wide range of inbreeding based on the excess of homozygosity (F_HOM) among individuals, from -0.04 to 0.16 in Sumava and from -0.13 to 0.12 in Wallachian. Average inbreeding based on runs of homozygosity was 0.21 in Sumava and 0.27 in Wallachian. Most detected runs of homozygosity (ROH) were less than 5 Mb long for both breeds. ROH segments longer than 15 Mb were absent in Wallachian sheep. Concerning putative selection signatures, a total of 471 candidate genes in Wallachian sheep within 11 hotspots and 653 genes within 13 hotspots in Sumava sheep were identified. Czech breeds appear to be well differentiated from each other and other European breeds. Their genetic diversity is low, especially in the case of the Wallachian breed. Sumava is not so threatened by low diversity but has a larger share of the non-native gene pool.

Temporal Monitoring of the Floreana Island Galapagos Giant Tortoise Captive Breeding Program

Captive breeding programs benefit from genetic analyses that identify relatedness between individuals, assign parentage to offspring, and track levels of genetic diversity. Monitoring these parameters across breeding cycles is critical to the success of a captive breeding program as it allows conservation managers to iteratively evaluate and adjust program structure. However, in practice, genetic tracking of breeding outcomes is rarely conducted. Here, we examined the first three offspring cohorts (2017-2020) of the genetically informed captive breeding program for the Floreana Island Galapagos giant tortoise, Chelonoidis niger. This captive breeding program is unique as the Floreana tortoise has been extinct since the 1800s, but its genome persists, in part, in the form of living hybrids with the extant Volcano Wolf tortoise, C. becki. Breeding over the study period took place at the Galapagos National Park Directorate breeding facility in four corrals, each containing three females and two males. Using 17 microsatellite markers, we were able to assign parentage to 94 of the 98 offspring produced over the study period. We observe that despite the addition of more founders since the pilot breeding program, the effective population size remains low, and changes to the arrangements of breeding corrals may be necessary to encourage more equal reproductive output from the males. This study demonstrates the value of hybrids for species restoration and the importance of continually reassessing the outcomes of captive breeding.

Genetic Diversity of Wild Boar and Deer

Genetic diversity provides the long-term capacity of species, communities, and the biosphere to persist under change [...].

Genetic Characterization and Alternative Preservation Ways of Locally Adapted Sheep Breeds: Cases of Private and Public Sheep Sectors in Tunisia and Italy

Simple Summary

The genetic characterization of native sheep breeds from the Tunisian-Italian Mediterranean corridor, the earliest westward introduction route of sheep breeding from the domestication center, was carried out using microsatellite markers in order to compare the genetic diversity level between both Mediterranean sides’ local sheep breeds, and to investigate their level of risk for conservation purposes. Sampling concerned both private and public/institutional farms of all the native Tunisian breeds in the case of the southern Mediterranean side, and one center Italian and the Venetian native breeds from the northern side. The results revealed that the genetic diversity level of the Tunisian native breeds was higher than that of the Italian ones, with a lower inbreeding level. The comparison between private and public farms in terms of genetic diversity, consanguinity, and conservation decisions highlighted the superiority of the public/institutional breeding management strategies over private ones in both Tunisian and Italian cases. The present study illustrated the efficiency of genomic characterization in making genetic diversity evaluations and preservation decisions of native and well-adapted sheep breeds in both developing and developed countries’ rearing conditions.

Abstract

Non-commercialized sheep breeds known as local or native breeds are well adapted to their environmental constraints and constitute precious genetic resources that need prioritization for genetic diversity characterization and preservation. The aim of the present study was to assess the genetic diversity level and the related preservation decisions of very old and traditional native Mediterranean sheep breeds from Tunisia and Italy using 17 microsatellite markers. In total, 975 sheep were sampled from five Tunisian, one Center Italian, and four Venetian native breeds. Both private and publicly available farms were considered for each breed for breeding strategies’ comparison purposes. The microsatellite set used was highly informative (PIC = 0.80 ± 0.08), with a total of 383 alleles. Moderate genetic differentiation was revealed between the native sheep of the two Mediterranean sides (global overall loci F_ST = 0.081). The genetic diversity level was higher in the case of the Tunisian native breeds compared to the Italian ones, as evidenced by higher mean allelic richness, higher expected and observed heterozygosities, and lower inbreeding levels. Priority for conservation suggestions was carried out for each private or public breed population based on the contribution of each population to the diversity of the whole data. The four Venetian breeds, already undergoing conservation, the Tunisian dairy breed, and the very ancient Maghrebian breed, would be favored for conservation. In conclusion, our results highlighted the importance of the analyzed Mediterranean native sheep breeds as valuable inherited genetic reservoirs and supported previous conservation decisions made for the threatened breeds.

African wild dogs (Lycaon pictus) from the Kruger National Park, South Africa are currently not inbred but have low genomic diversity

African wild dogs (Lycaon pictus) have undergone severe population reductions and are listed as endangered on the International Union for Conservation of Nature Red List. Small, isolated populations have the potential to suffer from threats to their genetic diversity that may impact species viability and future survival. This study provides the first set of population-wide genomic data to address conservation concerns for this endangered species. Whole genome sequencing data were generated for 71 free-ranging African wild dogs from the Kruger National Park (KNP), South Africa, and used to estimate important population genomic parameters. Genomic diversity metrics revealed that variation levels were low; however, this African wild dog population showed low levels of inbreeding. Very few first- and second-order relationships were observed in this cohort, with most relationships falling into the third-order or distant category. Patterns of homozygosity could have resulted from historical inbreeding or a loss in genome variation due to a population bottleneck. Although the results suggest that this stronghold African wild dog population maintains low levels of inbreeding, likely due to their cooperative breeding system, it may lead to a continuous population decline when a reduced number of suitable mates are available. Consequently, the low genomic variation may influence species viability over time. This study highlights the importance of assessing population genomic parameters to set conservation priorities. Future studies should include the investigation of the potential of this endangered species to adapt to environmental changes considering the low genomic diversity in this population.

Tunisian camel casein gene characterization reveals similarities and differences with Sudanese and Nigerian populations

Milk is a primary protein source that has always played a role in mammalian health. Despite the intensification of research projects on dromedary and the knowledge of the genetic diversity at the casein loci, the genetic structure of the Tunisian camel population still needs exploration. This study sought to determine the genetic diversity of 3 casein gene variants in 5 Tunisian camel ecotypes: c.150G>T at CSN1S1 (α_S1-casein), g.2126A>G at CSN2 (β-casein), and g.1029T>C at CSN3 (κ-casein). The obtained results were compared with data published on Sudanese and Nigerian camels to establish the level of differentiation within and between populations. A total of 159 blood samples were collected from 5 Tunisian camel ecotypes and the extracted DNA was genotyped by PCR-RFLP. A streamlined genotyping protocol was also developed for CSN3. Results indicated that allele T was quite rare (0.06) at CSN1S1 for all ecotypes. Minor allele frequency was found for G (0.462) in CSN2 except for Ardhaoui Medenine ecotype who deviated from the average CSN2 allele frequency of the total population. Allele C showed minor allele frequency of 0.384 in CSN3. Among the Tunisian population, GAT (0.343) was the most represented haplotype in all ecotypes except for Ardhaoui Medenine, where GGC (0.322) was the most frequent one. Significant differences in heterozygosity and local inbreeding were observed across the Tunisian, Sudanese, and Nigerian populations, although the global fixation index indicated that only 2.2% of the genetic variance is related to ecotype differences. Instead, phylogenetic analysis revealed a closer link between the Tunisian and Sudanese populations through a clade subdivision with 3 main branches among the ecotypes. This study represents the first attempt to understand casein gene variability in Tunisian camels; with further study, milk traits and genetic differentiation among populations can be associated with the history of camel domestication.

Molecular Markers: An Overview of Data Published for Fungi over the Last Ten Years

Fungi are amongst the most abundant and diverse organisms. Despite being widely known for their adverse role in food spoilage or as pathogens for humans, animals, or plants, they also present several beneficial effects. Fungi contribute to human well-being due to their role as decomposers, degrading decay matter into smaller molecules which can be easily used by other ecosystem members. These organisms can produce medicinal compounds or modulate protective immune responses in human intestine. Fungi intervene in diverse food processes or act as a food supply. Due to fungal diversity, the unequivocal identification of these organisms is crucial to increasing their practical applications and decreasing their adverse effects. The process of identification could be achieved through the integral sequencing of fungi genomes. However, this procedure would be time-consuming and rather cost-inefficient. Therefore, several molecular markers have been developed to overcome these limitations. The chronology of DNA-based molecular markers development can be divided into three main steps: (1) prior to the development of the PCR technique (RFLP); (2) after the development of the PCR technique (RAPD, AFLP, ISSR, VNTR, SNP, InDels, and DNA barcoding); (3) after the development of the massive parallel sequencing technique (Metabarcoding and WGS). Therefore, the present review covers an overview of the most recently developed molecular markers used for fungal detection and identification.

2b-RAD Genotyping of the Seagrass Cymodocea nodosa Along a Latitudinal Cline Identifies Candidate Genes for Environmental Adaptation

Plant populations distributed along broad latitudinal gradients often show patterns of clinal variation in genotype and phenotype. Differences in photoperiod and temperature cues across latitudes influence major phenological events, such as timing of flowering or seed dormancy. Here, we used an array of 4,941 SNPs derived from 2b-RAD genotyping to characterize population differentiation and levels of genetic and genotypic diversity of three populations of the seagrass Cymodocea nodosa along a latitudinal gradient extending across the Atlantic-Mediterranean boundary (i.e., Gran Canaria—Canary Islands, Faro—Portugal, and Ebro Delta—Spain). Our main goal was to search for potential outlier loci that could underlie adaptive differentiation of populations across the latitudinal distribution of the species. We hypothesized that such polymorphisms could be related to variation in photoperiod-temperature regime occurring across latitudes. The three populations were clearly differentiated and exhibited diverse levels of clonality and genetic diversity. Cymodocea nodosa from the Mediterranean displayed the highest genotypic richness, while the Portuguese population had the highest clonality values. Gran Canaria exhibited the lowest genetic diversity (as observed heterozygosity). Nine SNPs were reliably identified as outliers across the three sites by two different methods (i.e., BayeScan and pcadapt), and three SNPs could be associated to specific protein-coding genes by screening available C. nodosa transcriptomes. Two SNPs-carrying contigs encoded for transcription factors, while the other one encoded for an enzyme specifically involved in the regulation of flowering time, namely Lysine-specific histone demethylase 1 homolog 2. When analyzing biological processes enriched within the whole dataset of outlier SNPs identified by at least one method, “regulation of transcription” and “signalling” were among the most represented. Our results highlight the fundamental importance signal integration and gene-regulatory networks, as well as epigenetic regulation via DNA (de)methylation, could have for enabling adaptation of seagrass populations along environmental gradients.

Genetic diversity and population structure of six autochthonous pig breeds from Croatia, Serbia, and Slovenia

BACKGROUND

The importance of local breeds as genetic reservoirs of valuable genetic variation is well established. Pig breeding in Central and South-Eastern Europe has a long tradition that led to the formation of several local pig breeds. In the present study, genetic diversity parameters were analysed in six autochthonous pig breeds from Slovenia, Croatia and Serbia (Banija spotted, Black Slavonian, Turopolje pig, Swallow-bellied Mangalitsa, Moravka and Krskopolje pig). Animals from each of these breeds were genotyped using microsatellites and single nucleotide polymorphisms (SNPs). The results obtained with these two marker systems and those based on pedigree data were compared. In addition, we estimated inbreeding levels based on the distribution of runs of homozygosity (ROH) and identified genomic regions under selection pressure using ROH islands and the integrated haplotype score (iHS).

RESULTS

The lowest heterozygosity values calculated from microsatellite and SNP data were observed in the Turopolje pig. The observed heterozygosity was higher than the expected heterozygosity in the Black Slavonian, Moravka and Turopolje pig. Both types of markers allowed us to distinguish clusters of individuals belonging to each breed. The analysis of admixture between breeds revealed potential gene flow between the Mangalitsa and Moravka, and between the Mangalitsa and Black Slavonian, but no introgression events were detected in the Banija spotted and Turopolje pig. The distribution of ROH across the genome was not uniform. Analysis of the ROH islands identified genomic regions with an extremely high frequency of shared ROH within the Swallow-bellied Mangalitsa, which harboured genes associated with cholesterol biosynthesis, fatty acid metabolism and daily weight gain. The iHS approach to detect signatures of selection revealed candidate regions containing genes with potential roles in reproduction traits and disease resistance.

CONCLUSIONS

Based on the estimation of population parameters obtained from three data sets, we showed the existence of relationships among the six pig breeds analysed here. Analysis of the distribution of ROH allowed us to estimate the level of inbreeding and the extent of homozygous regions in these breeds. The iHS analysis revealed genomic regions potentially associated with phenotypic traits and allowed the detection of genomic regions under selection pressure.

How methodological changes have influenced our understanding of population structure in threatened species: insights from tiger populations across India

Unprecedented advances in sequencing technology in the past decade allow a better understanding of genetic variation and its partitioning in natural populations. Such inference is critical to conservation: to understand species biology and identify isolated populations. We review empirical population genetics studies of Endangered Bengal tigers within India, where 60–70% of wild tigers live. We assess how changes in marker type and sampling strategy have impacted inferences by reviewing past studies, and presenting three novel analyses including a single-nucleotide polymorphism (SNP) panel, genome-wide SNP markers, and a whole-mitochondrial genome network. At a broad spatial scale, less than 100 SNPs revealed the same patterns of population clustering as whole genomes (with the exception of one additional population sampled only in the SNP panel). Mitochondrial DNA indicates a strong structure between the northeast and other regions. Two studies with more populations sampled revealed further substructure within Central India. Overall, the comparison of studies with varied marker types and sample sets allows more rigorous inference of population structure. Yet sampling of some populations is limited across all studies, and these should be the focus of future sampling efforts. We discuss challenges in our understanding of population structure, and how to further address relevant questions in conservation genetics.

This article is part of the theme issue ‘Celebrating 50 years since Lewontin's apportionment of human diversity’.

Genome-wide SNPs detect fine-scale genetic structure in threatened populations of squirrel glider Petaurus norfolcensis

Australian arboreal mammals are experiencing significant population declines, particularly due to land clearing and resulting habitat fragmentation. The squirrel glider, Petaurus norfolcensis, is a threatened species in New South Wales, with a stronghold population in the Lake Macquarie Local Government Area (LGA) where fragmentation due to urbanization is an ongoing problem for the species conservation. Here we report on the use of squirrel glider mitochondrial (385 bp cytochrome b gene, 70 individuals) and nuclear DNA (6,834 SNPs, 87 individuals) markers to assess their population genetic structure and connectivity across 14 locations sampled in the Lake Macquarie LGA. The mitochondrial DNA sequences detected evidence of a historical genetic bottleneck, while the genome-wide SNPs detected significant population structure in the Lake Macquarie squirrel glider populations at scales as fine as one kilometer. There was no evidence of inbreeding within patches, however there were clear effects of habitat fragmentation and biogeographical barriers on gene flow. A least cost path analysis identified thin linear corridors that have high priority for conservation. These areas should be protected to avoid further isolation of squirrel glider populations and the loss of genetic diversity through genetic drift.

A reduced SNP panel to trace gene flow across southern European wolf populations and detect hybridization with other Canis taxa

Intra- and inter-specific gene flow are natural evolutionary processes. However, human-induced hybridization is a global conservation concern across taxa, and the development of discriminant genetic markers to differentiate among gene flow processes is essential. Wolves (Canis lupus) are affected by hybridization, particularly in southern Europe, where ongoing recolonization of historic ranges is augmenting gene flow among divergent populations. Our aim was to provide diagnostic canid markers focused on the long-divergent Iberian, Italian and Dinaric wolf populations, based on existing genomic resources. We used 158 canid samples to select a panel of highly informative single nucleotide polymorphisms (SNPs) to (i) distinguish wolves in the three regions from domestic dogs (C. l. familiaris) and golden jackals (C. aureus), and (ii) identify their first two hybrid generations. The resulting 192 SNPs correctly identified the five canid groups, all simulated first-generation (F1) hybrids (0.482 ≤ Q_i ≤ 0.512 between their respective parental groups) and all first backcross (BC1) individuals (0.723 ≤ Q_i ≤ 0.827 to parental groups). An assay design and test with invasive and non-invasive canid samples performed successfully for 178 SNPs. By separating natural population admixture from inter-specific hybridization, our reduced panel can help advance evolutionary research, monitoring, and timely conservation management.

Weaving place‐based knowledge for culturally significant species in the age of genomics: Looking to the past to navigate the future

Relationships with place provide critical context for characterizing biocultural diversity. Yet, genetic and genomic studies are rarely informed by Indigenous or local knowledge, processes, and practices, including the movement of culturally significant species. Here, we show how place‐based knowledge can better reveal the biocultural complexities of genetic or genomic data derived from culturally significant species. As a case study, we focus on culturally significant southern freshwater kōura (crayfish) in Aotearoa me Te Waipounamu (New Zealand, herein Aotearoa NZ). Our results, based on genotyping‐by‐sequencing markers, reveal strong population genetic structure along with signatures of population admixture in 19 genetically depauperate populations across the east coast of Te Waipounamu. Environment association and differentiation analyses for local adaptation also indicate a role for hydroclimatic variables—including temperature, precipitation, and water flow regimes—in shaping local adaptation in kōura. Through trusted partnerships between community and researchers, weaving genomic markers with place‐based knowledge has both provided invaluable context for the interpretation of data and created opportunities to reconnect people and place. We envisage such trusted partnerships guiding future genomic research for culturally significant species in Aotearoa NZ and beyond.

Population-level genome-wide STR discovery and validation for population structure and genetic diversity assessment of Plasmodium species

Short tandem repeats (STRs) are highly informative genetic markers that have been used extensively in population genetics analysis. They are an important source of genetic diversity and can also have functional impact. Despite the availability of bioinformatic methods that permit large-scale genome-wide genotyping of STRs from whole genome sequencing data, they have not previously been applied to sequencing data from large collections of malaria parasite field samples. Here, we have genotyped STRs using HipSTR in more than 3,000 Plasmodium falciparum and 174 Plasmodium vivax published whole-genome sequence data from samples collected across the globe. High levels of noise and variability in the resultant callset necessitated the development of a novel method for quality control of STR genotype calls. A set of high-quality STR loci (6,768 from P. falciparum and 3,496 from P. vivax) were used to study Plasmodium genetic diversity, population structures and genomic signatures of selection and these were compared to genome-wide single nucleotide polymorphism (SNP) genotyping data. In addition, the genome-wide information about genetic variation and other characteristics of STRs in P. falciparum and P. vivax have been available in an interactive web-based R Shiny application PlasmoSTR (https://github.com/bahlolab/PlasmoSTR).

Aspects of Molecular Genetics in Dromedary Camel

Dromedary camels are unique in their morphological and physiological characteristics and are capable of providing milk and meat even under extreme environmental conditions. Like other species, the dromedary camel has also benefitted from the development of the molecular genetics to increase the knowledge about different aspect in camel genetics (genetic variation, molecular marker, parentage control, gene of interest, whole genome, dating…etc.). In this paper we review the different molecular genetic technics used in this particular species and future prospects. Dromedary genetic studies started in the end of the 1980s with phenotypic evaluation and the attempts to highlight the protein and biochemical diversity. In the 2000s, with the development of molecular markers such as microsatellites, genetic diversity of different types in several countries were estimated and microsatellites were also used for parentage control. In terms of genetic characterization, microsatellites revealed a defined global structure, differentiating East African and South Arabian dromedaries from North African, North Arabian, and South Asian individuals, respectively. Also, mitochondrialDNA sequence analysis of ancient DNA proved to be crucial in resolving domestication processes in dromedaries. Ancient and modern DNA revealed dynamics of domestication and cross-continental dispersion of the dromedary. Nuclear SNPs, single nucleotide polymorphisms changes that occur approximately each 1000 bps in the mammalian genome were also applied in some studies in dromedary. These markers are a very useful alternative to microsatellites and have been employed in some studies on genetic diversity and relevant phenotypic traits in livestock. Finally, thanks to the use of Next Generation Sequencing (NGS) the whole-genome assemblies of the dromedary (Camelus dromedarius) and a work to establish the organization of the dromedary genome at chromosome level were recently published.

Waste not, want not: Microsatellites remain an economical and informative technology for conservation genetics

Comparisons of microsatellites and single-nucleotide polymorphisms (SNPs) have found that SNPs outperform microsatellites in population genetic analyses, questioning the continued utility of microsatellites in population and landscape genetics. Yet, highly polymorphic markers may be of value in species that have reduced genetic variation. This study repeated previous analyses that used microsatellites with SNPs developed from ddRAD sequencing in the black-capped vireo source-sink system. SNPs provided greater resolution of genetic diversity, population differentiation, and migrant detection but could not reconstruct parentage relationships due to insufficient heterozygosities. The biological inferences made by both sets of markers were similar: asymmetrical gene flow from source sites to the remaining sink sites. With the landscape genetic analyses, we found different results between the two molecular markers, but associations of the top environmental features (riparian, open habitat, agriculture, and human development) with dispersal estimates were shared between marker types. Despite the higher precision of SNPs, we find that microsatellites effectively uncover population processes and patterns and are superior for parentage analyses in this species with reduced genetic diversity. This study illustrates the continued applicability and relevance of microsatellites in population genetic research.

A Schistosoma mansoni tri- and tetramer microsatellite catalog for genetic population diversity and differentiation

All Schistosoma mansoni tri- and tetranucleotide repeat microsatellites published as of December 2018 were identified. All 52 were evaluated for autosomal location, strength of amplification, scorability and behavior as single-copy loci by polyacrylamide and capillary gel electrophoresis. Of these, 27 were unique, autosomal, polymorphic, easily scored and single copy as assessed on pooled adult worm DNA from two different continental origins and adult worm clones. These microsatellites were distributed across all seven autosomal chromosomes. On laboratory strains their heterozygosity ranged from 0.22 to 0.77. Individual markers had 5-13 alleles, allelic richness of 2-10 and an effective allele number of 1.3-8.14. Those infected by Schistosoma mansoni carry many genetically distinct, sexually reproducing parasites, therefore, for an individual infection the complete allele frequency profile of their progeny consists of a pool of DNA from multiple diploid eggs. Using a set of 25 microsatellites, we calculated allele frequency profiles of eggs in fecal samples from people in two Brazilian communities separated by 6 km: Jenipapo (n = 80) and Volta do Rio (n = 38). There were no a priori characteristics that could predict the performance of markers in natural infections based on their performance with laboratory strains. Increasing marker number did not change accuracy for differentiation and diversity but did improve precision. Our data suggest that using a random set of 10-20 microsatellites appears to result in values that exhibit low standard deviations for diversity and differentiation indices. All identified microsatellites as well as PCR conditions, allele size, primer sequences and references for all tri- and tetramer microsatellites markers presented in this work are available at: https://sites.google.com/case.edu/cwru-and-fiocruz-wdrc/home.

Harnessing the power of multi-omics data for predicting climate change response

Predicting how species will respond to future climate change is of central importance in the midst of the global biodiversity crisis, and recent work has demonstrated the utility of population genomics for improving these predictions. Here, we suggest a broadening of the approach to include other types of genomic variants that play an important role in adaptation, like structural (e.g. copy number variants) and epigenetic variants (e.g. DNA methylation). These data could provide additional power for forecasting response, especially in weakly structured or panmictic species. Incorporating structural and epigenetic variation into estimates of climate change vulnerability, or maladaptation, may not only improve prediction power but also provide insight into the molecular mechanisms underpinning species' response to climate change.

Genome-wide SNP analysis of three moose subspecies at the southern range limit in the contiguous United States

Genome-wide evaluations of genetic diversity and population structure are important for informing management and conservation of trailing-edge populations. North American moose (Alces alces) are declining along portions of the southern edge of their range due to disease, species interactions, and marginal habitat, all of which may be exacerbated by climate change. We employed a genotyping by sequencing (GBS) approach in an effort to collect baseline information on the genetic variation of moose inhabiting the species’ southern range periphery in the contiguous United States. We identified 1920 single nucleotide polymorphisms (SNPs) from 155 moose representing three subspecies from five states: A. a. americana (New Hampshire), A. a. andersoni (Minnesota), and A. a. shirasi (Idaho, Montana, and Wyoming). Molecular analyses supported three geographically isolated clusters, congruent with currently recognized subspecies. Additionally, while moderately low genetic diversity was observed, there was little evidence of inbreeding. Results also indicated > 20% shared ancestry proportions between A. a. shirasi samples from northern Montana and A. a. andersoni samples from Minnesota, indicating a putative hybrid zone warranting further investigation. GBS has proven to be a simple and effective method for genome-wide SNP discovery in moose and provides robust data for informing herd management and conservation priorities. With increasing disease, predation, and climate related pressure on range edge moose populations in the United States, the use of SNP data to identify gene flow between subspecies may prove a powerful tool for moose management and recovery, particularly if hybrid moose are more able to adapt.

One Step Away From Extinction: A Population Genomic Analysis of A Narrow Endemic, Tropical Plant Species

Intraspecific genetic variation plays a fundamental role in maintaining the evolutionary potential of wild populations. Hence, the assessment of genetic diversity patterns becomes essential to guide biodiversity conservation policies, particularly for threatened species. To inform management strategies for conservation of Mimosa catharinensis - a narrow endemic, critically endangered plant species - we identified 1,497 unlinked SNP markers derived from a reduced representation sequencing method (i.e., double digest restriction site associated DNA sequencing, or ddRADseq). This set of molecular markers was employed to assess intrapopulation genetic parameters and the demographic history of one extremely small population of M. catharinensis (N=33) located in the Brazilian Atlantic Forest. Contrary to what is expected for narrow endemic and threatened species with small population sizes, we observed a moderate level of genetic diversity for M. catharinensis [uH E(0%missing data)=0.205, 95% CI (0.160, 0.250); uH E(30%missing data)=0.233, 95% CI (0.174, 0.292)]. Interestingly, M. catharinensis, which is a lianescent shrub with no indication of seed production for at least two decades, presented high levels of outcrossing [t (0%missing data)=0.883, SE±0.0483; t (30%missing data)=0.909, SE±0.011] and an apparent absence of inbreeding [F (0%missing data)=-0.145, 95% CI (-0.189, -0.101); F (30%missing data)=-0.105, 95% CI (-0.199, -0.011)]. However, the reconstruction of demographic history of M. catharinensis indicated that the population should be suffered a recent bottleneck. Our population genomic study tackles a central issue in evolution and conservation biology and we expect that it will be useful to help safeguard the remaining genetic diversity reported for this unique genetic resource.

Chromosome-Level Genome Assemblies Expand Capabilities of Genomics for Conservation Biology

Genome assemblies are in the process of becoming an increasingly important tool for understanding genetic diversity in threatened species. Unfortunately, due to limited budgets typical for the area of conservation biology, genome assemblies of threatened species, when available, tend to be highly fragmented, represented by tens of thousands of scaffolds not assigned to chromosomal locations. The recent advent of high-throughput chromosome conformation capture (Hi-C) enables more contiguous assemblies containing scaffolds spanning the length of entire chromosomes for little additional cost. These inexpensive contiguous assemblies can be generated using Hi-C scaffolding of existing short-read draft assemblies, where N50 of the draft contigs is larger than 0.1% of the estimated genome size and can greatly improve analyses and facilitate visualization of genome-wide features including distribution of genetic diversity in markers along chromosomes or chromosome-length scaffolds. We compared distribution of genetic diversity along chromosomes of eight mammalian species, including six listed as threatened by IUCN, where both draft genome assemblies and newer chromosome-level assemblies were available. The chromosome-level assemblies showed marked improvement in localization and visualization of genetic diversity, especially where the distribution of low heterozygosity across the genomes of threatened species was not uniform.

Genetic structure of urban and non-urban populations differs between two common parid species

Landscape conversions induced by human activities can affect dispersal patterns of various bird species and, as a result, affect genetic structure of their populations. Genetic differentiation of bird populations may be enhanced by habitat variation, especially in urban-non-urban systems. The majority of population genetic studies focus on single species, which inflicts limitations for direct comparisons of genetic responses of avian populations to urbanization. Here, we used a set of microsatellite markers to examine genetic diversity, gene flow and population structure in two common parid species, great tits Parus major and blue tits Cyanistes caeruleus occupying three sites in habitats with contrasting urbanization level in central Poland. We found low but significant divergence of urban park population with both suburban and non-urban forest great tit populations, while no differentiation was found between suburban forest and non-urban forest populations. In contrast, no evidence for genetic differentiation was found between blue tit populations from the urban park, suburban forest and non-urban forest sites. We conclude that great tits and blue tits respond to urbanization-related changes in a different way, which may be a result of different rates of migration and/or dispersal, likely higher in blue tits. Some impact may be also induced by interspecific competition. We suggest that changing the focus of urban genetic research from single to multiple species may provide novel insights into how natural populations respond to the processes of urbanization.