Next-generation metrics for monitoring genetic erosion within populations of conservation concern

Testing the applicability of environmental DNA metabarcoding to landscape genetics

Landscape genetics is a field dealing with local genetic differences and contributes to strategic conservation planning. Recently, environmental DNA (eDNA) metabarcoding has proven useful not only for detecting species but also for assessing genetic diversity and genetic structure on a large scale such as in phylogeography. However, it remains unclear whether eDNA analysis also has sufficient power to perform the landscape genetics, which focuses on a local scale. To reveal the applicability of eDNA to landscape genetics, we conducted an eDNA metabarcoding analysis of the mitochondrial DNA D-loop region of the fluvial sculpin Cottus nozawae in the upper Sorachi River in Japan and compared the results with inferences based on traditional tissue-based approaches by the same D-loop region and genome-wide SNP data. As a result, the spatial distribution of haplotypes was generally consistent between the eDNA- and tissue-based approaches. In addition, the genetic differentiation statistics calculated using eDNA and tissue samples were highly correlated when comparing both in the D-loop region. The removal of low-frequency reads or the conversion to semi-quantitative rankings of eDNA data did not alter the correlation of genetic diversity and differentiation statistics with tissue-based approaches much. Finally, we confirmed that analyses using eDNA data can reveal patterns such as isolation-by-distance shown in previous studies on this species, indicating the applicability of eDNA to basic landscape genetics. Even though some limitations remain, eDNA may have great potential for conducting basic landscape genetics.

Discordance between taxonomy and population genomic data: An avian example relevant to the United States Endangered Species Act

Population genomics can reveal cryptic biological diversity that may impact fitness while simultaneously serving to delineate relevant conservation units. Here, we leverage the power of whole-genome resequencing for conservation by studying 433 individual lesser prairie-chicken (Tympanuchus pallidicinctus; LEPC, a federally endangered species of conservation concern in the United States) and greater prairie-chicken (Tympanuchus cupido; GRPC, a legally huntable species throughout much of its range). The genomic diversity of two formally recognized distinct population segments (DPSs) of LEPCs is similar, but they are genetically distinct. Neither DPS is depleted of its genomic diversity, neither is especially inbred, and temporal diversity is relatively stable in both conservation units. Interspecific differentiation between the two species was only slightly higher than that observed between LEPC DPSs, due largely to bidirectional introgression. The high resolution provided by our dataset identified a genomic continuum between the two species such that individuals sampled from the hybrid zone were imperfectly assigned to their presumptive species when considering only their physical characteristics. The admixture between the two species is reflected in the spectrum of individual ancestry coefficients, which has legal implications for the "take" of individuals under the Endangered Species Act. Overall, our data highlight the recurring dissonance between static policies and dynamic species boundaries that are increasingly obvious in the population genomic era.

Airports for the genetic rescue of a former agricultural pest

The intensification of agricultural practices and urbanisation are widespread causes of biodiversity loss. However, the role of artificial habitats in genetic rescue is an aspect that is not well understood. Implementing genetic rescue measures to improve gene flow and maintain a viable population of keystone species is a crucial prerequisite for promoting diverse and resilient ecosystems. Landscape fragmentation and modern agricultural methods have caused the decline and the isolation of the remnant colonies of the endangered European ground squirrel (Spermophilus citellus) throughout its range. However, the artificial habitat, such as airport fields with regular grass mowing, provides suitable conditions for this grassland specialist. We measured home range size and genetic variation of seven souslik colonies in western Slovakia. Based on the 6904 ddRAD SNPs, we found significantly higher individual heterozygosity in colonies on airports compared to colonies on pastures. This indicates a potential for higher fitness of individuals from airport colonies, which can serve as a source for evidence-based translocations. Such an intervention can preserve the genetic diversity of small and isolated populations in the region. We emphasize that conservation management strategies would be strengthened including a specific focus on human-made grassland habitats.

Estimates of heterozygosity from single nucleotide polymorphism markers are context-dependent and often wrong

Genetic diversity is frequently described using heterozygosity, particularly in a conservation context. Often, it is estimated using single nucleotide polymorphisms (SNPs); however, it has been shown that heterozygosity values calculated from SNPs can be biased by both study design and filtering parameters. Though solutions have been proposed to address these issues, our own work has found them to be inadequate in some circumstances. Here, we aimed to improve the reliability and comparability of heterozygosity estimates, specifically by investigating how sample size and missing data thresholds influenced the calculation of autosomal heterozygosity (heterozygosity calculated from across the genome, i.e. fixed and variable sites). We also explored how the standard practice of tri- and tetra-allelic site exclusion could bias heterozygosity estimates and influence eventual conclusions relating to genetic diversity. Across three distinct taxa (a frog, Litoria rubella; a tree, Eucalyptus microcarpa; and a grasshopper, Keyacris scurra), we found heterozygosity estimates to be meaningfully affected by sample size and missing data thresholds, partly due to the exclusion of tri- and tetra-allelic sites. These biases were inconsistent both between species and populations, with more diverse populations tending to have their estimates more severely affected, thus having potential to dramatically alter interpretations of genetic diversity. We propose a modified framework for calculating heterozygosity that reduces bias and improves the utility of heterozygosity as a measure of genetic diversity, whilst also highlighting the need for existing population genetic pipelines to be adjusted such that tri- and tetra-allelic sites be included in calculations.

Health of Holochilus chacarius (Rodentia: Cricetidae) in rice agroecosystem in a neotropical wetland assessed by histopathology

Small mammals have a short lifetime and are strictly associated with their environment. This work aimed to use histopathology to assess the health of Holochilus chacarius in a rice agroecosystem in the Pantanal of Mato Grosso do Sul. During necropsy, fragments of the lung, kidney, skin, liver, and reproductive system of 33 animals were collected and submitted to histological processing. Tissue damages were evaluated as mild, moderate, and severe and arranged in a matrix for further statistical analysis. Furthermore, we used generalized linear models to verify the influence of tissue changes on the body condition, obtained by a regression between body mass and length. In the lungs, we found an intense inflammatory infiltrate associated with anthracosis that had a negative influence on the body's condition. Also, we observed degenerative and inflammatory changes in the liver, kidneys, skin, and reproductive system that ranged from mild to moderate. The histopathological lesions observed in this study may be associated with environmental alterations of anthropic origin such as the exposure to soot from wildfires and heavy metals, evidenced by lesions in the lung, kidney, and liver. The present study provided a histopathological matrix as a new approach that allows to classify and quantify the tissue alterations. Tissue changes when associated with body condition demonstrated to be an effective tool to assess the health of small free-living mammals, showing that these animals can be used as bioindicators of environmental condition.

How Can Genomics Help or Hinder Wildlife Conservation?

Genomic data are becoming increasingly affordable and easy to collect, and new tools for their analysis are appearing rapidly. Conservation biologists are interested in using this information to assist in management and planning but are typically limited financially and by the lack of genomic resources available for non-model taxa. It is therefore important to be aware of the pitfalls as well as the benefits of applying genomic approaches. Here, we highlight recent methods aimed at standardizing population assessments of genetic variation, inbreeding, and forms of genetic load and methods that help identify past and ongoing patterns of genetic interchange between populations, including those subjected to recent disturbance. We emphasize challenges in applying some of these methods and the need for adequate bioinformatic support. We also consider the promises and challenges of applying genomic approaches to understand adaptive changes in natural populations to predict their future adaptive capacity.

Genomic insights into the Montseny brook newt (Calotriton arnoldi), a Critically Endangered glacial relict

The Montseny brook newt (Calotriton arnoldi), considered the most endangered amphibian in Europe, is a relict salamandrid species endemic to a small massif located in northeastern Spain. Although conservation efforts should always be guided by genomic studies, those are yet scarce among urodeles, hampered by the extreme sizes of their genomes. Here, we present the third available genome assembly for the order Caudata, and the first genomic study of the species and its sister taxon, the Pyrenean brook newt (Calotriton asper), combining whole-genome and ddRADseq data. Our results reveal significant demographic oscillations which accurately mirrored Europe's climatic history. Although severe bottlenecks have led to depauperate genomic diversity and long runs of homozygosity along a gigantic genome, inbreeding might have been avoided by assortative mating strategies. Other life history traits, however, seem to have been less advantageous, and the lack of land dispersal has driven to exceptional levels of population fragmentation.

Divergent mechanisms of reduced growth performance in Betula ermanii saplings from high-altitude and low-latitude range edges

The reduced growth performance of individuals from range edges is a common phenomenon in various taxa, and considered to be an evolutionary factor that limits the species' range. However, most studies did not distinguish between two mechanisms that can lead to this reduction: genetic load and adaptive selection to harsh conditions. To address this lack of understanding, we investigated the climatic and genetic factors underlying the growth performance of Betula ermanii saplings transplanted from 11 populations including high-altitude edge and low-latitude edge population. We estimated the climatic position of the populations within the overall B. ermanii's distribution, and the genetic composition and diversity using restriction-site associated DNA sequencing, and measured survival, growth rates and individual size of the saplings. The high-altitude edge population (APW) was located below the 95% significance interval for the mean annual temperature range, but did not show any distinctive genetic characteristics. In contrast, the low-latitude edge population (SHK) exhibited a high level of linkage disequilibrium, low genetic diversity, a distinct genetic composition from the other populations, and a high relatedness coefficient. Both APW and SHK saplings displayed lower survival rates, heights and diameters, while SHK saplings also exhibited lower growth rates than the other populations' saplings. The low heights and diameters of APW saplings was likely the result of adaptive selection to harsh conditions, while the low survival and growth rates of SHK saplings was likely the result of genetic load. Our findings shed light on the mechanisms underlying the reduced growth performance of range-edge populations.

Functional genomic diversity is correlated with neutral genomic diversity in populations of an endangered rattlesnake

Theory predicts that genetic erosion in small, isolated populations of endangered species can be assessed using estimates of neutral genetic variation, yet this widely used approach has recently been questioned in the genomics era. Here, we leverage a chromosome-level genome assembly of an endangered rattlesnake (Sistrurus catenatus) combined with whole genome resequencing data (N = 110 individuals) to evaluate the relationship between levels of genome-wide neutral and functional diversity over historical and future timescales. As predicted, we found positive correlations between genome-wide estimates of neutral genetic diversity (π) and inferred levels of adaptive variation and an estimate of inbreeding mutation load, and a negative relationship between neutral diversity and an estimate of drift mutation load. However, these correlations were half as strong for projected future levels of neutral diversity based on contemporary effective population sizes. Broadly, our results confirm that estimates of neutral genetic diversity provide an accurate measure of genetic erosion in populations of a threatened vertebrate. They also provide nuance to the neutral-functional diversity controversy by suggesting that while these correlations exist, anthropogenetic impacts may have weakened these associations in the recent past and into the future.

Evaluating the use of non-invasive hair sampling and ddRAD to characterize populations of endangered species: Application to a peripheral population of the European mink

The application of next-generation sequencing (NGS) to non-invasive samples is one of the most promising methods in conservation genomics, but these types of samples present significant challenges for NGS. The European mink (Mustela lutreola) is critically endangered throughout its range. However, important aspects such as census size and inbreeding remain still unknown in many populations, so it is crucial to develop new methods to monitor this species. In this work, we placed hair tubes along riverbanks in a border area of the Iberian population, which allowed the genetic identification of 76 European mink hair samples. We then applied a reduced representation genomic sequencing (ddRAD) technique to a subset of these samples to test whether we could extract sufficient genomic information from them. We show that several problems with the DNA, including contamination, fragmentation, oxidation, and possibly sample mixing, affected the samples. Using various bioinformatic techniques to reduce these problems, we were able to unambiguously genotype 19 hair samples belonging to six individuals. This small number of individuals showed that the demographic status of the species in this peripheral population is worse than expected. The data obtained also allowed us to perform preliminary analyses of relatedness and inbreeding. Although further improvements in sampling and analysis are needed, the application of the ddRAD technique to non-invasively obtained hairs represents a significant advance in the genomic study of endangered species.

How genomics can help biodiversity conservation

The availability of public genomic resources can greatly assist biodiversity assessment, conservation, and restoration efforts by providing evidence for scientifically informed management decisions. Here we survey the main approaches and applications in biodiversity and conservation genomics, considering practical factors, such as cost, time, prerequisite skills, and current shortcomings of applications. Most approaches perform best in combination with reference genomes from the target species or closely related species. We review case studies to illustrate how reference genomes can facilitate biodiversity research and conservation across the tree of life. We conclude that the time is ripe to view reference genomes as fundamental resources and to integrate their use as a best practice in conservation genomics.

Conservation genomics reveals fine-scale population structuring and recent declines in the Critically Endangered Australian Kuranda Treefrog

The Kuranda Treefrog occurs in tropical north-east Australia and is listed as Critically Endangered due to its small distribution and population size, with observed declines due to drought and human-associated impacts to habitat. Field surveys identified marked population declines in the mid-2000s, culminating in very low abundance at most sites in 2005 and 2006, followed by limited recovery. Here, samples from before (2001–2004) and after (2007–2009) this decline were analysed using 7132 neutral genome-wide SNPs to assess genetic connectivity among breeding sites, genetic erosion, and effective population size. We found a high level of genetic connectivity among breeding sites, but also structuring between the population at the eastern end of the distribution (Jumrum Creek) versus all other sites. Despite finding no detectable sign of genetic erosion between the two times periods, we observed a marked decrease in effective population size (Ne), from 1720 individuals pre-decline to 818 post-decline. This mirrors the decline detected in the field census data, but the magnitude of the decline suggested by the genetic data is greater. We conclude that the current effective population size for the Kuranda Treefrog remains around 800 adults, split equally between Jumrum Creek and all other sites combined. The Jumrum Creek habitat requires formal protection. Connectivity among all other sites must be maintained and improved through continued replanting of rainforest, and it is imperative that impacts to stream flow and water quality are carefully managed to maintain or increase population sizes and prevent genetic erosion.

Understanding Local Adaptation to Prepare Populations for Climate Change

Adaptation within species to local environments is widespread in nature. Better understanding this local adaptation is critical to conserving biodiversity. However, conservation practices can rely on species’ trait averages or can broadly assume homogeneity across the range to inform management. Recent methodological advances for studying local adaptation provide the opportunity to fine-tune efforts for managing and conserving species. The implementation of these advances will allow us to better identify populations at greatest risk of decline because of climate change, as well as highlighting possible strategies for improving the likelihood of population persistence amid climate change. In the present article, we review recent advances in the study of local adaptation and highlight ways these tools can be applied in conservation efforts. Cutting-edge tools are available to help better identify and characterize local adaptation. Indeed, increased incorporation of local adaptation in management decisions may help meet the imminent demands of managing species amid a rapidly changing world.

Challenges in quantifying genome erosion for conservation

Massive defaunation and high extinction rates have become characteristic of the Anthropocene. Genetic effects of population decline can lead populations into an extinction vortex, where declining populations show lower genetic fitness, in turn leading to lower populations still. The lower genetic fitness in a declining population due to a shrinking gene pool is known as genetic erosion. Three different types of genetic erosion are highlighted in this review: overall homozygosity, genetic load and runs of homozygosity (ROH), which are indicative of inbreeding. The ability to quantify genetic erosion could be a very helpful tool for conservationists, as it can provide them with an objective, quantifiable measure to use in the assessment of species at risk of extinction. The link between conservation status and genetic erosion should become more apparent. Currently, no clear correlation can be observed between the current conservation status and genetic erosion. However, the high quantities of genetic erosion in wild populations, especially in those species dealing with habitat fragmentation and habitat decline, may be early signs of deteriorating populations. Whole genome sequencing data is the way forward to quantify genetic erosion. Extra screening steps for genetic load and hybridization can be included, since they could potentially have great impact on population fitness. This way, the information yielded from genetic sequence data can provide conservationists with an objective genetic method in the assessment of species at risk of extinction. However, the great complexity of genome erosion quantification asks for consensus and bridging science and its applications, which remains challenging.

Population genomic monitoring provides insight into conservation status but no correlation with demographic estimates of extinction risk in a threatened trout

The current extinction crisis requires effective assessment and monitoring tools. Genetic approaches are appealing given the relative ease of field sampling required to estimate genetic diversity characteristics assumed related to population size, evolutionary potential, and extinction risk, and to evaluate hybridization with non‐native species simultaneously. However, linkages between population genetic metrics of diversity from survey‐style field collections and demographic estimates of population size and extinction risk are still in need of empirical examples, especially for remotely distributed species of conservation concern where the approach might be most beneficial. We capitalized on an exceptional opportunity to evaluate congruence between genetic diversity metrics and demographic‐based estimates of abundance and extinction risk from a comprehensive Multiple Population Viability Analysis (MPVA) in a threatened fish, the Lahontan cutthroat trout (LCT). We sequenced non‐native trout reference samples and recently collected and archived tissue samples of most remaining populations of LCT (N = 60) and estimated common genetic assessment metrics, predicting minimal hybridization with non‐native trout, low diversity, and declining diversity over time. We further hypothesized genetic metrics would correlate positively with MPVA‐estimated abundance and negatively with extinction probability. We uncovered several instances of hybridization that pointed to immediate management needs. After removing hybridized individuals, cautious interpretation of low effective population sizes (2–63) suggested reduced evolutionary potential for many LCT populations. Other genetic metrics did not decline over time nor correlate with MPVA‐based estimates of harmonic mean abundance or 30‐year extinction probability. Our results demonstrate benefits of genetic monitoring for efficiently detecting hybridization and, though genetic results were disconnected from demographic assessment of conservation status, they suggest reduced evolutionary potential and likely a higher conservation risk than currently recognized for this threatened fish. We emphasize that genetic information provides essential complementary insight, in addition to demographic information, for evaluating species status.

Global genetic diversity status and trends: towards a suite of Essential Biodiversity Variables (EBVs) for genetic composition

Biodiversity underlies ecosystem resilience, ecosystem function, sustainable economies, and human well‐being. Understanding how biodiversity sustains ecosystems under anthropogenic stressors and global environmental change will require new ways of deriving and applying biodiversity data. A major challenge is that biodiversity data and knowledge are scattered, biased, collected with numerous methods, and stored in inconsistent ways. The Group on Earth Observations Biodiversity Observation Network (GEO BON) has developed the Essential Biodiversity Variables (EBVs) as fundamental metrics to help aggregate, harmonize, and interpret biodiversity observation data from diverse sources. Mapping and analyzing EBVs can help to evaluate how aspects of biodiversity are distributed geographically and how they change over time. EBVs are also intended to serve as inputs and validation to forecast the status and trends of biodiversity, and to support policy and decision making. Here, we assess the feasibility of implementing Genetic Composition EBVs (Genetic EBVs), which are metrics of within‐species genetic variation. We review and bring together numerous areas of the field of genetics and evaluate how each contributes to global and regional genetic biodiversity monitoring with respect to theory, sampling logistics, metadata, archiving, data aggregation, modeling, and technological advances. We propose four Genetic EBVs: (i) Genetic Diversity; (ii) Genetic Differentiation; (iii) Inbreeding; and (iv) Effective Population Size (N_e). We rank Genetic EBVs according to their relevance, sensitivity to change, generalizability, scalability, feasibility and data availability. We outline the workflow for generating genetic data underlying the Genetic EBVs, and review advances and needs in archiving genetic composition data and metadata. We discuss how Genetic EBVs can be operationalized by visualizing EBVs in space and time across species and by forecasting Genetic EBVs beyond current observations using various modeling approaches. Our review then explores challenges of aggregation, standardization, and costs of operationalizing the Genetic EBVs, as well as future directions and opportunities to maximize their uptake globally in research and policy. The collection, annotation, and availability of genetic data has made major advances in the past decade, each of which contributes to the practical and standardized framework for large‐scale genetic observation reporting. Rapid advances in DNA sequencing technology present new opportunities, but also challenges for operationalizing Genetic EBVs for biodiversity monitoring regionally and globally. With these advances, genetic composition monitoring is starting to be integrated into global conservation policy, which can help support the foundation of all biodiversity and species' long‐term persistence in the face of environmental change. We conclude with a summary of concrete steps for researchers and policy makers for advancing operationalization of Genetic EBVs. The technical and analytical foundations of Genetic EBVs are well developed, and conservation practitioners should anticipate their increasing application as efforts emerge to scale up genetic biodiversity monitoring regionally and globally.

Conservation Genetics of Mediterranean Brown Trout in Central Italy (Latium): A Multi-Marker Approach

Brown trout is considered a complex of incipient species, including several phylogenetic lineages, whose natural distribution in the Mediterranean area has been altered, since the beginning of the 1900s, by massive introductions of domestic strains of Atlantic origin to support fisheries. Introduced trout naturalize in new suitable environments and extensively hybridize with native populations. Here, we characterized putatively neutral and adaptive genetic variability and population structure of Mediterranean brown trout from six river catchments in central peninsular Italy, as revealed by both mitochondrial (Control Region) and nuclear (microsatellites, LDH-C1, major histocompatibility complex) markers. We quantified the admixture of wild populations with hatchery strains and evaluated the effects of domestic trout introductions on shaping population genetics. Our analyses indicated: (1) a composite picture of genetic variability in the area, with the presence of all native Mediterranean trout mitochondrial lineages (“Adriatic”, “Mediterranean”, “marmoratus”), various frequencies of allochthonous genotypes and different rates of introgression among sampling sites; (2) asymmetric mito-nuclear introgression; (3) increasing nuclear marker diversity with increasing levels of admixture across populations; (4) strong population structure coupled with relatively low effective population size. Data allowed the identification of five management units and we propose specific actions to support ongoing and future conservation strategies within the examined area.

Supporting in situ conservation of the genetic diversity of crop wild relatives using genomic technologies

The last decade has witnessed huge technological advances in genomics, particularly in DNA sequencing. Here, we review the actual and potential application of genomics in supporting in situ conservation of crop wild relatives (CWRs). In addition to helping in prioritization of protection of CWR taxa and in situ conservation sites, genome analysis is allowing the identification of novel alleles that need to be prioritized for conservation. Genomics is enabling the identification of potential sources of important adaptive traits that can guide the establishment or enrichment of in situ genetic reserves. Genomic tools also have the potential for developing a robust framework for monitoring and reporting genome‐based indicators of genetic diversity changes associated with factors such as land use or climate change. These tools have been demonstrated to have an important role in managing the conservation of populations, supporting sustainable access and utilization of CWR diversity, enhancing accelerated domestication of new crops and forensic genomics thus preventing misappropriation of genetic resources. Despite this great potential, many policy makers and conservation managers have failed to recognize and appreciate the need to accelerate the application of genomics to support the conservation and management of biodiversity in CWRs to underpin global food security. Funding and inadequate genomic expertise among conservation practitioners also remain major hindrances to the widespread application of genomics in conservation.

THC and CBD Fingerprinting of an Elite Cannabis Collection from Iran: Quantifying Diversity to Underpin Future Cannabis Breeding

Cannabis (Cannabis sativa L.) has a rich history of human use, and the therapeutic importance of compounds produced by this species is recognized by the medical community. The active constituents of cannabis, collectively called cannabinoids, encompass hundreds of distinct molecules, the most well-characterized of which are tetrahydrocannabinol (THC) and cannabidiol (CBD), which have been used for centuries as recreational drugs and medicinal agents. As a first step to establish a cannabis breeding program, we initiated this study to describe the HPLC-measured quantity of THC and CBD biochemistry profiles of 161 feral pistillate cannabis plants from 20 geographical regions of Iran. Our data showed that Iran can be considered a new region of high potential for distribution of cannabis landraces with diverse THC and CBD content, predominantly falling into three groups, as Type I = THC-predominant, Type II = approximately equal proportions of THC and CBD (both CBD and THC in a ratio close to the unity), and Type III = CBD-predominant. Correlation analysis among two target cannabinoids and environmental and geographical variables indicated that both THC and CBD contents were strongly influenced by several environmental-geographical factors, such that THC and CBD contents were positively correlated with mean, min and max annual temperature and negatively correlated with latitude, elevation, and humidity. Additionally, a negative correlation was observed between THC and CBD concentrations, suggesting that further studies to unravel these genotype × environment interactions (G × E interactions) are warranted. The results of this study provide important pre-breeding information on a collection of cannabis that will underpin future breeding programs.

Crop genetic erosion: understanding and responding to loss of crop diversity

Crop diversity underpins the productivity, resilience and adaptive capacity of agriculture. Loss of this diversity, termed crop genetic erosion, is therefore concerning. While alarms regarding evident declines in crop diversity have been raised for over a century, the magnitude, trajectory, drivers and significance of these losses remain insufficiently understood. We outline the various definitions, measurements, scales and sources of information on crop genetic erosion. We then provide a synthesis of evidence regarding changes in the diversity of traditional crop landraces on farms, modern crop cultivars in agriculture, crop wild relatives in their natural habitats and crop genetic resources held in conservation repositories. This evidence indicates that marked losses, but also maintenance and increases in diversity, have occurred in all these contexts, the extent depending on species, taxonomic and geographic scale, and region, as well as analytical approach. We discuss steps needed to further advance knowledge around the agricultural and societal significance, as well as conservation implications, of crop genetic erosion. Finally, we propose actions to mitigate, stem and reverse further losses of crop diversity.

Genomic signatures of inbreeding in a critically endangered parrot, the kākāpō

Events of inbreeding are inevitable in critically endangered species. Reduced population sizes and unique life-history traits can increase the severity of inbreeding, leading to declines in fitness and increased risk of extinction. Here, we investigate levels of inbreeding in a critically endangered flightless parrot, the kākāpō (Strigops habroptilus), wherein a highly inbred island population and one individual from the mainland of New Zealand founded the entire extant population. Genotyping-by-sequencing (GBS), and a genotype calling approach using a chromosome-level genome assembly, identified a filtered set of 12,241 single-nucleotide polymorphisms (SNPs) among 161 kākāpō, which together encompass the total genetic potential of the extant population. Multiple molecular-based estimates of inbreeding were compared, including genome-wide estimates of heterozygosity (FH), the diagonal elements of a genomic-relatedness matrix (FGRM), and runs of homozygosity (RoH, FRoH). In addition, we compared levels of inbreeding in chicks from a recent breeding season to examine if inbreeding is associated with offspring survival. The density of SNPs generated with GBS was sufficient to identify chromosomes that were largely homozygous with RoH distributed in similar patterns to other inbred species. Measures of inbreeding were largely correlated and differed significantly between descendants of the two founding populations. However, neither inbreeding nor ancestry was found to be associated with reduced survivorship in chicks, owing to unexpected mortality in chicks exhibiting low levels of inbreeding. Our study highlights important considerations for estimating inbreeding in critically endangered species, such as the impacts of small population sizes and admixture between diverse lineages.

Genomic signatures of inbreeding and mutation load in a threatened rattlesnake

Theory predicts that threatened species living in small populations will experience high levels of inbreeding that will increase their genetic load, but recent work suggests that the impact of load may be minimized by purging resulting from long-term population bottlenecks. Empirical studies that examine this idea using genome-wide estimates of inbreeding and genetic load in threatened species are limited. Here we use individual genome resequencing data to compare levels of inbreeding, levels of genetic load (estimated as mutation load) and population history in threatened Eastern massasauga rattlesnakes (Sistrurus catenatus), which exist in small isolated populations, and closely related yet outbred Western massasauga rattlesnakes (Sistrurus tergeminus). In terms of inbreeding, S. catenatus genomes had a greater number of runs of homozygosity of varying sizes, indicating sustained inbreeding through repeated bottlenecks when compared to S. tergeminus. At the species level, outbred S. tergeminus had higher genome-wide levels of mutation load in the form of greater numbers of derived deleterious mutations compared to S. catenatus, presumably due to long-term purging of deleterious mutations in S. catenatus. In contrast, mutations that escaped species-level drift effects within S. catenatus populations were in general more frequent and more often found in homozygous genotypes than in S. tergeminus, suggesting a reduced efficiency of purifying selection in smaller S. catenatus populations for most mutations. Our results support an emerging idea that the historical demography of a threatened species has a significant impact on the type of genetic load present, which impacts implementation of conservation actions such as genetic rescue.

Influence of Holocene habitat availability on Pacific gray whale (Eschrichtius robustus) population dynamics as inferred from whole mitochondrial genome sequences and environmental niche modeling

Environmental changes since the Pleistocene and commercial whaling in the last few centuries have drastically reduced many whale populations, including gray whales in the North Pacific. Herein we use complete mitogenome sequences from 74 individuals to evaluate gray whale phylogeography and historical demography, then use environmental niche modeling to assess how habitat availability has changed through time for Pacific gray whales. We identify a large degree of haplotype sharing between gray whales sampled in Russian and Mexican waters, coupled with very limited matrilineal population structure. Confirming previous studies, our environmental niche models showed a decrease in available habitat during the Last Glacial Maximum, but we find no genetic signals of recent population declines in mitochondrial genomes despite both sustained habitat loss and a commercial whaling bottleneck. Our results illustrate the complex dynamics of baleen whale biogeography since the Holocene as well as the difficulty in detecting recent demographic bottlenecks from mitochondrial DNA sequences.

The genome of the Pyrenean desman and the effects of bottlenecks and inbreeding on the genomic landscape of an endangered species

The Pyrenean desman (Galemys pyrenaicus) is a small semiaquatic mammal endemic to the Iberian Peninsula. Despite its limited range, this species presents a strong genetic structure due to past isolation in glacial refugia and subsequent bottlenecks. Additionally, some populations are highly fragmented today as a consequence of river barriers, causing substantial levels of inbreeding. These features make the Pyrenean desman a unique model in which to study the genomic footprints of differentiation, bottlenecks and extreme isolation in an endangered species. To understand these processes, the complete genome of the Pyrenean desman was sequenced and assembled using a Bloom filter-based approach. An analysis of the 1.83 Gb reference genome and the sequencing of five additional individuals from different evolutionary units allowed us to detect its main genomic characteristics. The population differentiation of the species was reflected in highly distinctive demographic trajectories. In addition, a severe population bottleneck during the postglacial recolonization of the eastern Pyrenees created one of the lowest genomic heterozygosity values recorded in a mammal. Moreover, isolation and inbreeding gave rise to a high proportion of runs of homozygosity (ROH). Despite these extremely low levels of genetic diversity, two key multigene families from an eco-evolutionary perspective, the major histocompatibility complex and olfactory receptor genes, showed heterozygosity excess in the majority of individuals, revealing that functional diversity can be maintained up to a certain extent. Furthermore, these two classes of genes were significantly less abundant than expected within ROH. In conclusion, the genomic landscape of each analysed Pyrenean desman turned out to be strikingly distinctive and was a clear reflection of its recent ancestry and current conservation conditions. These results may help characterize the genomic health of each individual, and can be crucial for the conservation and management of the species.

The endangered White Sands pupfish (Cyprinodon tularosa) genome reveals low diversity and heterogenous patterns of differentiation

The White Sands pupfish (Cyprinodon tularosa), endemic to New Mexico in Southwestern North America, is of conservation concern due in part to invasive species, chemical pollution, and groundwater withdrawal. Herein, we developed a draft reference genome and use it to provide biological insights into the evolution and conservation of C. tularosa. We used our assembly to localize microsatellite markers previously used to demarcate evolutionary significant units (ESU), quantified genomic divergence and transposable element profiles between species, and compared C. tularosa genomic diversity related species. Our de novo assembly of PacBio Sequel II error-corrected reads resulted in a 1.08 Gb draft genome with a contig N50 of 1.4 Mb and 25,260 annotated protein coding genes, including 95% of the expected Actinopterygii conserved complete single-copy orthologues. Many of the C. tularosa microsatellite markers used for conservation assessments fell within, or near, genes and exhibited a pattern of increased heterozygosity near genic areas compared to those in intergenic regions. Nuclear alignments between these two species revealed 193 genes contained in rapidly diverging tracts; transposable element profiles were largely concordant and suggest a shared, rapid expansion of LINE and Gypsy elements. Genome-wide heterozygosity was markedly lower in C. tularosa compared to estimates from other related species, probably because of smaller long-term effective population sizes constrained by their isolated and limited habitat. Overall, these inferences provide new insights into C. tularosa that should help inform future management efforts.

Genetic load has potential in large populations but is realized in small inbred populations

Populations with higher genetic diversity and larger effective sizes have greater evolutionary capacity (i.e., adaptive potential) to respond to ecological stressors. We are interested in how the variation captured in protein-coding genes fluctuates relative to overall genomic diversity and whether smaller populations suffer greater costs due to their genetic load of deleterious mutations compared with larger populations. We analyzed individual whole-genome sequences (N = 74) from three different populations of Montezuma quail (Cyrtonyx montezumae), a small ground-dwelling bird that is sustainably harvested in some portions of its range but is of conservation concern elsewhere. Our historical demographic results indicate that Montezuma quail populations in the United States exhibit low levels of genomic diversity due in large part to long-term declines in effective population sizes over nearly a million years. The smaller and more isolated Texas population is significantly more inbred than the large Arizona and the intermediate-sized New Mexico populations we surveyed. The Texas gene pool has a significantly smaller proportion of strongly deleterious variants segregating in the population compared with the larger Arizona gene pool. Our results demonstrate that even in small populations, highly deleterious mutations are effectively purged and/or lost due to drift. However, we find that in small populations the realized genetic load is elevated because of inbreeding coupled with a higher frequency of slightly deleterious mutations that are manifested in homozygotes. Overall, our study illustrates how population genomics can be used to proactively assess both neutral and functional aspects of contemporary genetic diversity in a conservation framework while simultaneously considering deeper demographic histories.

New microsatellite markers for Ellipse, Venustaconcha ellipsiformis (Bivalvia: Unionidae), with notes on optimal sample size and cross-species amplification

Venustaconcha ellipsiformis (Unionidae) is a freshwater mussel species inhabiting small to medium streams of the Midwestern United States; however, its occurrence is rather sporadic and populations are often isolated. Due to anthropogenic habitat degradation and water pollution, this species is designated as some sort of conservation status in many states. To prioritize conservation strategies, highly variable genetic markers are necessary to assess population genetic structure and potential genetic erosion of V. ellipsiformis. Using whole genome sequence data, we developed and characterized microsatellite markers for V. ellipsiformis. Among 23 tetranucleotide loci tested, 14 loci were consistently amplified and showed polymorphism. Analyses performed on three populations in the upper Mississippi River basin showed that the number of alleles per locus ranged from 2 to 11 and the observed heterozygosity varied from 0.15 to 0.75. Based on genotypic and allelic rarefaction curves, these loci had adequate statistical power to genetically discriminate between individuals and the sample size was large enough to capture most alleles available in the populations at most loci. Finally, cross-species screening of the loci successfully amplified and showed polymorphism in six species in the tribe Lampsilini. The microsatellite loci developed in this study provide a valuable addition to extant genetic markers for freshwater mussels and can be useful to provide high-level resolution of population genetic parameters for V. ellipsiformis. Such information will be of great value for resource managers developing and prioritizing conservation strategies for imperiled mussel species.

Marine mammal conservation: over the horizon

Marine mammals can play important ecological roles in aquatic ecosystems, and their presence can be key to community structure and function. Consequently, marine mammals are often considered indicators of ecosystem health and flagship species. Yet, historical population declines caused by exploitation, and additional current threats, such as climate change, fisheries bycatch, pollution and maritime development, continue to impact many marine mammal species, and at least 25% are classified as threatened (Critically Endangered, Endangered or Vulnerable) on the IUCN Red List. Conversely, some species have experienced population increases/recoveries in recent decades, reflecting management interventions, and are heralded as conservation successes. To continue these successes and reverse the downward trajectories of at-risk species, it is necessary to evaluate the threats faced by marine mammals and the conservation mechanisms available to address them. Additionally, there is a need to identify evidence-based priorities of both research and conservation needs across a range of settings and taxa. To that effect we: (1) outline the key threats to marine mammals and their impacts, identify the associated knowledge gaps and recommend actions needed; (2) discuss the merits and downfalls of established and emerging conservation mechanisms; (3) outline the application of research and monitoring techniques; and (4) highlight particular taxa/populations that are in urgent need of focus.

Genomic scans for selective sweeps through haplotype homozygosity and allelic fixation in 14 indigenous sheep breeds from Middle East and South Asia

The performance and productivity of livestock have consistently improved by natural and artificial selection over the centuries. Both these selections are expected to leave patterns on the genome and lead to changes in allele frequencies, but natural selection has played the major role among indigenous populations. Detecting selective sweeps in livestock may assist in understanding the processes involved in domestication, genome evolution and discovery of genomic regions associated with economically important traits. We investigated population genetic diversity and selection signals in this study using SNP genotype data of 14 indigenous sheep breeds from Middle East and South Asia, including six breeds from Iran, namely Iranian Balochi, Afshari, Moghani, Qezel, Zel, and Lori-Bakhtiari, three breeds from Afghanistan, namely Afghan Balochi, Arabi, and Gadik, three breeds from India, namely Indian Garole, Changthangi, and Deccani, and two breeds from Bangladesh, namely Bangladeshi Garole and Bangladesh East. The SNP genotype data were generated by the Illumina OvineSNP50 Genotyping BeadChip array. To detect genetic diversity and population structure, we used principal component analysis (PCA), admixture, phylogenetic analyses, and Runs of homozygosity. We applied four complementary statistical tests, FST (fixation index), xp-EHH (cross-population extended haplotype homozygosity), Rsb (extended haplotype homozygosity between-populations), and FLK (the extension of the Lewontin and Krakauer) to detect selective sweeps. Our results not only confirm the previous studies but also provide a suite of novel candidate genes involved in different traits in sheep. On average, FST, xp-EHH, Rsb, and FLK detected 128, 207, 222, and 252 genomic regions as candidates for selective sweeps, respectively. Furthermore, nine overlapping candidate genes were detected by these four tests, especially TNIK, DOCK1, USH2A, and TYW1B which associate with resistance to diseases and climate adaptation. Knowledge of candidate genomic regions in sheep populations may facilitate the identification and potential exploitation of the underlying genes in sheep breeding.

The risk status of Polish local breeds under conservation programmes – new approach

Native breeds of farm animals, despite their unique characteristics and adaptability to diverse and often adverse environmental conditions, are at risk due to their lower profitability and corresponding declining use. According to FAO data, there are currently around 8800 livestock breeds in the world (http://www.fao.org/dad-is/en/ update: 09.11.2019), of which 17% have been assessed to be endangered. The objective of the study was to develop a model for estimating the risk status of Polish native breeds. The proposed new model is based on 3 factors: demographic, genetic, and socioeconomic. The last factor is composed of 6 elements: geographical concentration, demographic trend over the last 5 years, cultural value of a breed, parentage control, ex situ conservation and anthropogenic factors. Using the factors and elements above, among the examined 30 breeds of cattle, horses, pigs, sheep and goats, 1 was assessed as being critically endangered; 24 were determined to be endangered and in need of conservation; 5 were less endangered, but in need of constant monitoring. Out of 35 populations of poultry 29 were assessed as endangered as they were being kept in only single public collections. Six flocks of ducks were given a critical risk status, due to a very low number of females.

A selection signatures study among Middle Eastern and European sheep breeds

Selection, both natural and artificial, leaves patterns on the genome during domestication of animals and leads to changes in allele frequencies among populations. Detecting genomic regions influenced by selection in livestock may assist in understanding the processes involved in genome evolution and discovering genomic regions related to traits of economic and ecological interests. In the current study, genetic diversity analyses were conducted on 34,206 quality-filtered SNP positions from 450 individuals in 15 sheep breeds, including six indigenous breeds from the Middle East, namely Iranian Balouchi, Afshari, Moghani, Qezel, Karakas and Norduz, and nine breeds from Europe, namely East Friesian Sheep, Ile de France, Mourerous, Romane, Swiss Mirror, Spaelsau, Suffolk, Comisana and Engadine Red Sheep. The SNP genotype data generated by the Illumina OvineSNP50 Genotyping BeadChip array were used in this analysis. We applied two complementary statistical analyses, F_ST (fixation index) and xp-EHH (cross-population extended haplotype homozygosity), to detect selection signatures in Middle Eastern and European sheep populations. F_ST and xp-EHH detected 629 and 256 genes indicating signatures of selection, respectively. Genomic regions identified using F_ST and xp-EHH contained the CIDEA, HHATL, MGST1, FADS1, RTL1 and DGKG genes, which were reported earlier to influence a number of economic traits. Both F_ST and xp-EHH approaches identified 60 shared genes as the signatures of selection, including four candidate genes (NT5E, ADA2, C8A and C8B) that were enriched for two significant Gene Ontology (GO) terms associated with the adenosine metabolic procedure. Knowledge about the candidate genomic regions under selective pressure in sheep breeds may facilitate identification of the underlying genes and enhance our understanding on these genes role in local adaptation.

Genetic Rescue and the Plight of Ponui Hybrids

Long-term sustainable and resilient populations is a key goal of conservation. How to best achieve this is controversial. There are, for instance, polarized views concerning the fitness and conservation value of hybrid populations founded through multi-origin translocations. A classic example concerns Apteryx (kiwi) in New Zealand. The A. mantelli of Ponui Island constitute a hybrid population where the birds are highly successful in their island habitat. A key dilemma for managers is understanding the reason for this success. Are the hybrid birds of Ponui Island of “no future conservation value” as recently asserted, or do they represent an outstanding example of genetic rescue and an important resource for future translocations? There has been a paradigm shift in scientific thinking concerning hybrids, but the ecological significance of admixed genomes remains difficult to assess. This limits what we can currently predict in conservation science. New understanding from genome science challenges the sufficiency of population genetic models to inform decision making and suggests instead that the contrasting outcomes of hybridization, “outbreeding depression” and “heterosis,” require understanding additional factors that modulate gene and protein expression and how these factors are influenced by the environment. We discuss these findings and the investigations that might help us to better understand the birds of Ponui, inform conservation management of kiwi and provide insight relevant for the future survival of Apteryx.

An alternative method for determining the body condition index of the free-living South American coati

Assessing and monitoring the welfare of free-living mammals is not a usual process due to the logistical complications associated with their capture and sedation, collection and storage of biological samples and their release. In this context, non-invasive methods for monitoring wildlife constitute a good alternative approach for in situ conservation. Body condition index, as a measurement of health status, has been used in free-living mammals; its low value may be associated with negative effects on reproduction and survival. The present study aimed to generate an alternative and reliable non-invasive method and then determine the body condition index, based on previously-collected biometric measurements, without the need to capture and immobilise the animals. A total of 178 free-living Nasua nasua Linnaeus, 1766 were trapped, weighed and measured. Statistical methods were used, based on Boosted Regression Trees (BRT) using body mass, biometric measurements (body length, height and chest girth) and gender as explanatory variables. To assess the agreement between the real Body Condition Indices (BCIs) and the predicted values of BCIs, we explored the correlation between each model using the Bland-Altman method. This method showed a strong agreement between the predictive BRT models proposed (standardised residuals from a linear regression between body length and chest girth) and standardised residuals (linear regression between body mass and body length). The results obtained herein showed that BRT modelling, based on biometrical features, is an alternative way to verify the body conditions of coatis without the need to capture and immobilise the animals.

Examination of the efficacy of small genetic panels in genomic conservation of companion animal populations

In many ways, dogs are an ideal model for the study of genetic erosion and population recovery, problems of major concern in the field of conservation genetics. Genetic diversity in many dog breeds has been declining systematically since the beginning of the 1800s, when modern breeding practices came into fashion. As such, inbreeding in domestic dog breeds is substantial and widespread and has led to an increase in recessive deleterious mutations of high effect as well as general inbreeding depression. Pedigrees can in theory be used to guide breeding decisions, though are often incomplete and do not reflect the full history of inbreeding. Small microsatellite panels are also used in some cases to choose mating pairs to produce litters with low levels of inbreeding. However, the long-term impact of such practices has not been thoroughly evaluated. Here, we use forward simulation on a model of the dog genome to examine the impact of using limited marker panels to guide pairwise mating decisions on genome-wide population-level genetic diversity. Our results suggest that in unmanaged populations, where breeding decisions are made at the pairwise-rather than population-level, such panels can lead to accelerated loss of genetic diversity at genome regions unlinked to panel markers, compared to random mating. These results demonstrate the importance of genome-wide genetic panels for managing and conserving genetic diversity in dogs and other companion animals.

Population genomics for wildlife conservation and management

Biodiversity is under threat worldwide. Over the past decade, the field of population genomics has developed across nonmodel organisms, and the results of this research have begun to be applied in conservation and management of wildlife species. Genomics tools can provide precise estimates of basic features of wildlife populations, such as effective population size, inbreeding, demographic history and population structure, that are critical for conservation efforts. Moreover, population genomics studies can identify particular genetic loci and variants responsible for inbreeding depression or adaptation to changing environments, allowing for conservation efforts to estimate the capacity of populations to evolve and adapt in response to environmental change and to manage for adaptive variation. While connections from basic research to applied wildlife conservation have been slow to develop, these connections are increasingly strengthening. Here we review the primary areas in which population genomics approaches can be applied to wildlife conservation and management, highlight examples of how they have been used, and provide recommendations for building on the progress that has been made in this field.

Genetic Diversity in the Portuguese Mertolenga Cattle Breed Assessed by Pedigree Analysis

Simple Summary

The conservation and maintenance of genetic diversity is one of the priorities of the Convention of Biological Diversity and is included in the United Nations (UN’s) Sustainable Development Goals. The evaluation of the genetic variability of a breed is fundamental for its future use in a sustainable way, being indispensable to outline a successful conservation or improvement strategy. Preserving genetic diversity in a population is one of the main objectives for a breed conservation program. Nevertheless, the correct management of genetic diversity is also essential for the adaptation of a population to a new environment, production system or genetic improvement. For the purpose of population monitoring, assessing changes in genetic variability and genetic erosion in animal populations, many methodologies based on pedigree analyses of inbreeding and relationships, and on the probability of genetic origin from different herds, founders and ancestors, have been used. This study presents several genetic diversity indicators in a Portuguese native cattle breed, Mertolenga, assessed by pedigree analysis, and demonstrates the usefulness of these indicators and how they can be used in the genetic management of a breed.

Abstract

The Mertolenga beef cattle, currently with 27,000 breeding females in Portugal, is the largest Portuguese native breed, despite some variation in the breeding stock over the last years. The purpose of this study was to estimate parameters related to the population structure and genetic diversity and to investigate the major factors affecting genetic erosion in the breed, based on the pedigree herdbook information collected since the 1950s, including records on 221,567 animals from 425 herds. The mean generation intervals were 6.4 years for sires and 7.1 years for dams, respectively. The rate of inbreeding per year was 0.183% ± 0.020% and the correspondent effective population size was 38.83. In the reference population (35,017 calves born between 2015 and 2019), the average inbreeding and relatedness were 8.82% ± 10% and 2.05% ± 1.26%, respectively. The mean relationship among animals from the same and from different herds was 29.25% ± 9.36% and 1.87% ± 1.53%, respectively. The estimates for the effective number of founders, ancestors, founding herds and herds supplying sires were 87.9, 59.4, 21.4 and 73.5, respectively. Although the situation of the Mertolenga breed is not alarming, these results indicate the need to adopt measures to maintain the genetic variability of the population.

Whole-exome sequencing reveals a long-term decline in effective population size of red spruce (Picea rubens)

Understanding the factors influencing the current distribution of genetic diversity across a species range is one of the main questions of evolutionary biology, especially given the increasing threat to biodiversity posed by climate change. Historical demographic processes such as population expansion or bottlenecks and decline are known to exert a predominant influence on past and current levels of genetic diversity, and revealing this demo-genetic history can have immediate conservation implications. We used a whole-exome capture sequencing approach to analyze polymorphism across the gene space of red spruce (Picea rubens Sarg.), an endemic and emblematic tree species of eastern North America high-elevation forests that are facing the combined threat of global warming and increasing human activities. We sampled a total of 340 individuals, including populations from the current core of the range in northeastern USA and southeastern Canada and from the southern portions of its range along the Appalachian Mountains, where populations occur as highly fragmented mountaintop "sky islands." Exome capture baits were designed from the closely relative white spruce (P. glauca Voss) transcriptome, and sequencing successfully captured most regions on or near our target genes, resulting in the generation of a new and expansive genomic resource for studying standing genetic variation in red spruce applicable to its conservation. Our results, based on over 2 million exome-derived variants, indicate that red spruce is structured into three distinct ancestry groups that occupy different geographic regions of its highly fragmented range. Moreover, these groups show small Ne , with a temporal history of sustained population decline that has been ongoing for thousands (or even hundreds of thousands) of years. These results demonstrate the broad potential of genomic studies for revealing details of the demographic history that can inform management and conservation efforts of nonmodel species with active restoration programs, such as red spruce.

MHC Genotyping by SSCP and Amplicon-Based NGS Approach in Chamois

Genes of the major histocompatibility complex (MHC) code for cell surface proteins essential for adaptive immunity. They show the most outstanding genetic diversity in vertebrates, which has been connected with various fitness traits and thus with the long-term persistence of populations. In this study, polymorphism of the MHC class II DRB locus was investigated in chamois with Single-Strand Conformation Polymorphism (SSCP)/Sanger genotyping and Ion Torrent S5 next-generation sequencing (NGS). From eight identified DRB variants in 28 individuals, five had already been described, and three were new, undescribed alleles. With conventional SSCP/Sanger sequencing, we were able to detect seven alleles, all of which were also detected with NGS. We found inconsistencies in the individual genotypes between the two methods, which were mainly caused by allelic dropout in the SSCP/Sanger method. Six out of 28 individuals were falsely classified as homozygous with SSCP/Sanger analysis. Overall, 25% of the individuals were identified as genotyping discrepancies between the two methods. Our results show that NGS technologies are better performing in sequencing highly variable regions such as the MHC, and they also have a higher detection capacity, thus allowing a more accurate description of the genetic composition, which is crucial for evolutionary and population genetic studies.

Genetic diversity in Iranian indigenous sheep vis-à-vis selected exogenous sheep breeds and wild mouflon

The heterogeneity of climate and different agro-ecological conditions in Iran have resulted in development of 27 indigenous sheep breeds. Wild Asiatic mouflon (Ovis orientalis) is believed to be the ancestor of Iranian sheep. Evaluation of genetic diversity and population structure within and among domestic breeds has important implications for animal breeding programs and genetic resources management. Based on 50K SNP genotype data, we studied the genetic diversity of five indigenous Iranian sheep breeds: Afshari (n = 37), Moghani (n = 34), Qezel (n = 35), Zel (n = 46) and Lori-Bakhtiari (n = 46), and Asiatic mouflon (n = 8) sampled from Iran. Furthermore, genetic diversity and the breed admixture of Iranian sheep were assessed on a larger geographic scale using a reference panel comprising: three indigenous Afghan breeds - Arabi (n = 15), Balouchi (n = 15) and Gadik (n = 15); three indigenous breeds from Turkey and Cyprus - Cyprus Fat Tail (n = 30), Karakas (n = 18) and Norduz (n = 20); and three commercial European breeds - Suffolk (n = 19), Comisana (n = 24) and Engadine Red Sheep (n = 24). The results revealed that the investigated breeds are divided into five genetically distinct clusters according to their geographic origin. Afshari was closest to the local mouflon population and showed signs of mouflon admixture. Qezel was identified as a hybrid sheep breed. Much evidence supported the Afghan breeds being identical. Inbreeding values, which were estimated based on ROHs, were highest for Suffolk (F_ROH  = 0.0544) and lowest for Balouchi (F_ROH  = 0.0078). In conclusion, analysis of selected breeds from neighboring countries along with Asiatic mouflon gave a deeper insight into the evolutionary history and origin of Iranian sheep with important implications for future breed management.