High connectivity among argali sheep from Afghanistan and adjacent countries: Inferences from neutral and candidate gene microsatellites

Genome-wide diversity, population structure and signatures of inbreeding in the African buffalo in Mozambique

The African buffalo, Syncerus caffer, is a key species in African ecosystems. Like other large herbivores, it plays a fundamental role in its habitat acting as an ecosystem engineer. Over the last few centuries, African buffalo populations have declined because of range contraction and demographic decline caused by direct or indirect human activities. In Mozambique, historically home to large buffalo herds, the combined effect of colonialism and subsequent civil wars has created a critical situation that urgently needs to be addressed. In this study, we focused on the analysis of genetic diversity of Syncerus caffer caffer populations from six areas of Mozambique. Using genome-wide SNPs obtained from ddRAD sequencing, we examined the population structure across the country, estimated gene flow between areas under conservation management, including national reserves, and assessed the inbreeding coefficients. Our results indicate that all studied populations of Syncerus caffer caffer are genetically depauperate, with a high level of inbreeding. Moreover, buffaloes in Mozambique present a significant population differentiation between southern and central areas. We found an unexpected genotype in the Gorongosa National Park, where buffaloes experienced a dramatic population size reduction, that shares a common ancestry with southern populations of Catuane and Namaacha. This could suggest the past occurrence of a connection between southern and central Mozambique and that the observed population structuring could reflect recent events of anthropogenic origin. All the populations analysed showed high levels of homozygosity, likely due to extensive inbreeding over the last few decades, which could have increased the frequency of recessive deleterious alleles. Improving the resilience of Syncerus caffer caffer in Mozambique is essential for preserving the ecosystem integrity. The most viable approach appears to be facilitating translocations and re-establishing connectivity between isolated herds. However, our results also highlight the importance of assessing intraspecific genetic diversity when considering interventions aimed at enhancing population viability such as selecting suitable source populations.

Border fences reduce potential for transboundary migration of Marco Polo Sheep (Ovis ammon polii) in the Pamir Plateau

Border fences have severely impeded the transboundary migration of a number of large mammals worldwide, with central Asia being one of the most impacted. The Marco Polo sheep (Ovis ammon polii), an iconic species of Pamir, is threatened in its transboundary movement by increasing border fencing among their five distributed countries, including Tajikistan, Kyrgyzstan, China, Afghanistan, and Pakistan. In this study, by building ensemble species distribution models, we found that eastern Tajikistan had the largest suitable Macro Polo sheep habitat (about 42 % of the total suitable habitat), followed by China (about 32 %). We used least-cost paths to identify 51 ecological corridors including 5 transboundary ecological corridors, which may be important to maintain connectivity in both domestic and transboundary regions. To assess the potential barrier effect of border fences, we assessed four scenarios (30, 40, 50 and 60°) corresponding to the upper limit of the slope for the construction of fences. In areas too steep for fencing, these could be used by wild sheep to cross barriers or borders and may represent migration or movement routes, defined as natural passages. In the most pessimistic Scenario 60, only 25 migratory passages along the border fences were identified, compared to 997 in the most optimistic scenario (Scenario 30), indicating a strong negative effect of intensive border fencing on the transboundary movement of Marco Polo sheep. The establishment of transnational conservation parks, and ensuring permeability is maintained in key areas, could have a positive impact on the connectivity and persistence of Marco Polo sheep populations, and provide important lessons for other large migratory mammals in transboundary regions.

The European Ground Squirrel’s Genetic Diversity in Its Ancestral Land: Landscape Insights and Conservation Implications

The European ground squirrel is an endangered rodent whose populations are declining throughout its range. Only in Bulgaria, the genetic hotspot of the species, are some abundant populations still present. We employed 12 microsatellite loci in ten Bulgarian populations to look at population structure, gene flow and recent bottlenecks. We found that the populations are in good condition in terms of heterozygosity, where values ranged from 0.55 to 0.78. However, the inbreeding index (FIS) was significant for most populations. A recent bottleneck was detected in only one population. Based on Bayesian clustering methods, the populations in Bulgaria were attributed to two groups, northern and southern, with admixture in the northern one. The AMOVA test between these groups showed no differentiation in genetic diversity. The mean value of FST was 0.184, which shows strong diversification among all populations. Hence, gene flow is probably limited. All these results indicate that Bulgaria is the main area to focus the efforts for conservation of the species by ensuring that the complex and rich genetic structure of Bulgarian populations is preserved.

Genetic Structure of the Ca Rater Mallorquí Dog Breed Inferred by Microsatellite Markers

Ca Rater Mallorquí is a dog breed from the Island of Mallorca (Spain) traditionally used as a hunting and ratting dog to prevent disease spread and economic losses related to rodent activities on farms. However, the census data shows a population decline that should be addressed by implementing a conservation program. The first step to implementing a conservation plan is knowing the genetic situation of the Ca Rater Mallorquí population. Therefore, we aimed to genetically characterise the breed in our study. We analysed 33 microsatellites recommended by the International Society of Animal Genetics (ISAG) in 77 samples. Data were obtained from 13 samples of Balearic, Spanish, and international dog breeds to study the genetic diversity among breeds. The population did not significantly deviate from the Hardy-Weinberg equilibrium with heterozygosity (Ho) of 0.655 and expected heterozygosity (He) of 0.685. The Wright's fixation indices, the Factorial Correspondence Analysis (FCA), a dendrogram representing Reynolds genetic distance between populations, and the pairwise F_ST values establish the Ca Rater Mallorquí as an independent breed distinct from the Balearic, Spanish, and international breeds.

On the origins of American Criollo pigs: A common genetic background with a lasting Iberian signature

American Criollo pigs are thought to descend mainly from those imported from the Iberian Peninsula starting in the late 15th century. Criollo pigs subsequently expanded throughout the Americas, adapting to very diverse environments, and possibly receiving influences from other origins. With the intensification of agriculture in the mid-20th century, cosmopolitan breeds largely replaced Criollo pigs, and the few remaining are mostly maintained by rural communities in marginal areas where they still play an important socio-economic and cultural role. In this study, we used 24 microsatellite markers in samples from 1715 pigs representing 46 breeds with worldwide distribution, including 17 American Criollo breeds, with the major focus of investigating their genetic diversity, structure and breed relationships. We also included representatives of the Iberian, Local British, Hungarian, Chinese and Commercial breeds, as well as Wild Boar, in order to investigate their possible influence in the genetic composition of Criollos. Our results show that, when compared with the other breeds, Criollo pigs present higher levels of genetic diversity, both in terms of allelic diversity and expected heterozygosity. The various analyses indicate that breed differentiation overall explains nearly 21% of the total genetic diversity. Criollo breeds showed their own identity and shared a common genetic background, tending to cluster together in various analyses, even though they differ from each other. A close relationship of Criollos with Iberian breeds was revealed by all the different analyses, and the contribution of Iberian breeds, particularly of the Celtic breeds, is still present in various Criollo breeds. No influence of Chinese breeds was detected on Criollos, but a few were influenced by Commercial breeds or by wild pigs. Our results confirm the uniqueness of American Criollo pigs and the role that Iberian breeds have played in their development.

Spatial genetic structure of Rocky Mountain bighorn sheep (Ovis canadensis canadensis) at the northern limit of their native range

The Canadian Rocky Mountains are one of the few places on Earth where the spatial genetic structure of wide-ranging species has been relatively unaffected by anthropogenic disturbance. We characterized the spatial genetic structure of Rocky Mountain bighorn sheep (Ovis canadensis canadensis Shaw, 1804) in the northern portion of their range. Using microsatellites from 1495 individuals and mitochondrial DNA sequences from 188 individuals, we examined both broad- and fine-scale spatial genetic structure, assessed sex-biased gene flow within the northern portion of the species range, and identified geographic patterns of genetic diversity. We found that broad-scale spatial genetic structure was consistent with barriers to movement created by major river valleys. The fine-scale spatial genetic structure was characterized by a strong isolation-by-distance pattern, and analysis of neighborhood size using spatial autocorrelation indicated gene flow frequently occurred over distances of up to 100 km. However, analysis of sex-specific spatial autocorrelation and analysis of mitochondrial haplotype distributions failed to detect any evidence of sex-biased gene flow. Finally, our analyses reveal decreasing genetic diversity with increasing latitude, consistent with patterns of post-glacial recolonization of the Rocky Mountains.

Bighorn sheep gut microbiomes associate with genetic and spatial structure across a metapopulation

Studies in laboratory animals demonstrate important relationships between environment, host traits, and microbiome composition. However, host-microbiome relationships in natural systems are understudied. Here, we investigate metapopulation-scale microbiome variation in a wild mammalian host, the desert bighorn sheep (Ovis canadensis nelsoni). We sought to identify over-represented microbial clades and understand how landscape variables and host traits influence microbiome composition across the host metapopulation. To address these questions, we performed 16S sequencing on fecal DNA samples from thirty-nine bighorn sheep across seven loosely connected populations in the Mojave Desert and assessed relationships between microbiome composition, environmental variation, geographic distribution, and microsatellite-derived host population structure and heterozygosity. We first used a phylogenetically-informed algorithm to identify bacterial clades conserved across the metapopulation. Members of genus Ruminococcaceae, genus Lachnospiraceae, and family Christensenellaceae R7 group were among the clades over-represented across the metapopulation, consistent with their known roles as rumen symbionts in domestic livestock. Additionally, compositional variation among hosts correlated with individual-level geographic and genetic structure, and with population-level differences in genetic heterozygosity. This study identifies microbiome community variation across a mammalian metapopulation, potentially associated with genetic and geographic population structure. Our results imply that microbiome composition may diverge in accordance with landscape-scale environmental and host population characteristics.

Assessing changes in functional connectivity in a desert bighorn sheep metapopulation after two generations

Determining how species move across complex and fragmented landscapes and interact with human-made barriers is a major research focus in conservation. Studies estimating functional connectivity from movement, dispersal or gene flow usually rely on a single study period and rarely consider variation over time. We contrasted genetic structure and gene flow across barriers for a metapopulation of desert bighorn sheep (Ovis canadensis nelsoni) using genotypes collected 2000-2003 and 2013-2015. Based on the recently observed but unexpected spread of a respiratory pathogen across an interstate highway previously identified as a barrier to gene flow, we hypothesized that bighorn sheep changed how they interacted with that barrier, and that shifts in metapopulation structure influenced gene flow, genetic diversity and connectivity. Population assignment tests, genetic structure and genetic recapture demonstrated that bighorn sheep crossed the interstate highway in at least one location in 2013-2015, sharply reducing genetic structure between two populations, but supported conclusions of an earlier study that such crossings were very infrequent or unknown in 2000-2003. A recently expanded population established new links and caused decreases in genetic structure among multiple populations. Genetic diversity showed only slight increases in populations linked by new connections. Genetic structure and assignments revealed other previously undetected changes in movements and distribution, but much was consistent. Thus, we observed changes in both structural and functional connectivity over just two generations, but only in specific locations. Movement patterns of species should be revisited periodically to enable informed management, particularly in dynamic and fragmented systems.

Do Marco Polo argali Ovis ammon polii persist in Pakistan?

The distribution range of the Near Threatened Marco Polo argali, or Marco Polo sheep, Ovis ammon polii is restricted to the Pamir Mountains, spanning Afghanistan, Tajikistan, China and Pakistan. Until the early 1970s the Marco Polo argali was abundant in northern areas of Pakistan, particularly in the Khunjerab and Misgar Valleys around the Pamir Knot, bordering China, Afghanistan and Tajikistan. In Pakistan the subspecies now occurs only in one small watershed, in Khunjerab National Park, bordering China, which it visits sporadically during summer. We used map-based questionnaire surveys, double-observer surveys and camera trapping in a search for Marco Polo argali in the Pakistani Pamirs. We observed a herd of 19 individuals in Karachanai Nallah, in Khunjerab National Park, in 2012. The population that was formerly present in the Misgar Valley appears to have been extirpated, presumably as a result of anthropogenic factors such as poaching, competition with livestock, habitat disturbance, and the construction of a fence along the international border with China. Transboundary migration and range expansion into Pakistan could be facilitated by removal of the border fence adjacent to the Kilik–Mintika area and by limiting livestock grazing in former lambing areas.

SNP discovery in candidate adaptive genes using exon capture in a free-ranging alpine ungulate

Identification of genes underlying genomic signatures of natural selection is key to understanding adaptation to local conditions. We used targeted resequencing to identify SNP markers in 5321 candidate adaptive genes associated with known immunological, metabolic and growth functions in ovids and other ungulates. We selectively targeted 8161 exons in protein-coding and nearby 5' and 3' untranslated regions of chosen candidate genes. Targeted sequences were taken from bighorn sheep (Ovis canadensis) exon capture data and directly from the domestic sheep genome (Ovis aries v. 3; oviAri3). The bighorn sheep sequences used in the Dall's sheep (Ovis dalli dalli) exon capture aligned to 2350 genes on the oviAri3 genome with an average of 2 exons each. We developed a microfluidic qPCR-based SNP chip to genotype 476 Dall's sheep from locations across their range and test for patterns of selection. Using multiple corroborating approaches (lositan and bayescan), we detected 28 SNP loci potentially under selection. We additionally identified candidate loci significantly associated with latitude, longitude, precipitation and temperature, suggesting local environmental adaptation. The three methods demonstrated consistent support for natural selection on nine genes with immune and disease-regulating functions (e.g. Ovar-DRA, APC, BATF2, MAGEB18), cell regulation signalling pathways (e.g. KRIT1, PI3K, ORRC3), and respiratory health (CYSLTR1). Characterizing adaptive allele distributions from novel genetic techniques will facilitate investigation of the influence of environmental variation on local adaptation of a northern alpine ungulate throughout its range. This research demonstrated the utility of exon capture for gene-targeted SNP discovery and subsequent SNP chip genotyping using low-quality samples in a nonmodel species.

Why and how might genetic and phylogenetic diversity be reflected in the identification of key biodiversity areas?

‘Key biodiversity areas' are defined as sites contributing significantly to the global persistence of biodiversity. The identification of these sites builds from existing approaches based on measures of species and ecosystem diversity and process. Here, we therefore build from the work of Sgró et al. (2011 Evol. Appl.4, 326–337. (doi:10.1111/j.1752-4571.2010.00157.x)) to extend a framework for how components of genetic diversity might be considered in the identification of key biodiversity areas. We make three recommendations to inform the ongoing process of consolidating a key biodiversity areas standard: (i) thresholds for the threatened species criterion currently consider a site's share of a threatened species' population; expand these to include the proportion of the species' genetic diversity unique to a site; (ii) expand criterion for ‘threatened species' to consider ‘threatened taxa’ and (iii) expand the centre of endemism criterion to identify as key biodiversity areas those sites holding a threshold proportion of the compositional or phylogenetic diversity of species (within a taxonomic group) whose restricted ranges collectively define a centre of endemism. We also recommend consideration of occurrence of EDGE species (i.e. threatened phylogenetic diversity) in key biodiversity areas to prioritize species-specific conservation actions among sites.