Relatedness-based mate choice and female philopatry: inbreeding trends of wolf packs in a human-dominated landscape.

Inbreeding can reduce offspring fitness and has substantial implications for the genetic diversity and long-term viability of populations. In social cooperative canids, inbreeding is conditioned by the geographic proximity between opposite-sex kin outside natal groups and the presence of related individuals in neighbouring groups. Consequently, challenges in moving into other regions where the species is present can also affect inbreeding rates. These can be particularly problematic in areas of high human density, where movement can be restricted, even for highly vagile species. In this study, we investigate the socio-ecological dynamics of Iberian wolf packs in the human-dominated landscape of Alto Minho, in northwest Portugal, where wolves exhibit a high prevalence of short-distance dispersal and limited gene flow with neighbouring regions. We hypothesise that mating occurs regardless of relatedness, resulting in recurrent inbreeding due to high kin encounter rates. Using data from a 10-year non-invasive genetic monitoring programme and a combination of relatedness estimates and genealogical reconstructions, we describe genetic diversity, mate choice, and dispersal strategies among Alto Minho packs. In contrast with expectations, our findings reveal relatedness-based mate choice, low kin encounter rates, and a reduced number of inbreeding events. We observed a high prevalence of philopatry, particularly among female breeders, with the most common breeding strategy involving the pairing of a philopatric female with an unrelated immigrant male. Overall, wolves were not inbred, and temporal changes in genetic diversity were not significant. Our findings are discussed, considering the demographic trend of wolves in Alto Minho and its human-dominated landscape.

The population bottleneck of the Iberian wolf impacted genetic diversity but not admixture with domestic dogs: A temporal genomic approach.

After decades of intense persecution, the Iberian wolf subspecies faced a severe bottleneck in the 1970s that considerably reduced its range and population size, nearly leading to its extinction in central and southern Iberian Peninsula. Such population decline could have impacted the genetic diversity of Iberian wolves through different processes, namely genetic drift and dynamics of hybridization with domestic dogs. By contrasting the genomes of 68 contemporary with 54 historical samples spanning the periods before and immediately after the 1970s bottleneck, we found evidence of its impact on genetic diversity and dynamics of wolf-dog hybridization. Our genome-wide assessment revealed that wolves and dogs form two well-differentiated genetic groups in Iberia and that hybridization rates did not increase during the bottleneck. However, an increased number of hybrid individuals was found over time during the population re-expansion, particularly at the edge of the wolf range. We estimated a low percentage of dog ancestry (~1.4%) in historical samples, suggesting that dog introgression was not a key driver for wolf extinction in central and southern Iberia. Our findings also unveil a significant decline in genetic diversity in contemporary samples, with the highest proportion of homozygous segments in the genome being recently inherited. Overall, our study provides unprecedented insight into the impact of a sharp decline on the Iberian wolf genome and refines our understanding of the ecological and evolutionary drivers of wolf-dog hybridization in the wild.

Evolutionary history of an island endemic, the Azorean common quail.

Oceanic islands are characterized by conditions that favour diversification into endemic lineages that can be very different from their mainland counterparts. This can be the result of fast phenotypic divergence due to drift or the result of slower adaptation to local conditions. This uniqueness can obscure their evolutionary history. Here we used morphological, stable isotope, genetic and genomic data to characterize common quails (Coturnix coturnix) in the Azores archipelago and assess the divergence from neighbouring common quail populations. Historical documents suggested that these quails could have a recent origin associated with the arrival of humans in the last centuries. Our results show that Azorean quails constitute a well-differentiated lineage with small size and dark throat pigmentation that has lost the migratory ability and that diverged from mainland quail lineages more than 0.8 mya, contrary to the notion of a recent human-mediated arrival. Even though some Azorean quails carry an inversion that affects 115 Mbp of chromosome 1 and that has been associated with the loss of the migratory behaviour in other common quail populations, half of the analysed individuals do not have that inversion and still do not migrate. The long coexistence and evolution in isolation in the Azores of two chromosomal variants (with and without the inversion) is best explained by balancing selection. Thus, a unique and long evolutionary history led to the island endemic that we know today, C. c. conturbans.

African climate and geomorphology drive evolution and ghost introgression in sable antelope.

The evolutionary history of African ungulates has been explained largely in the light of Pleistocene climatic oscillations and the way these influenced the distribution of vegetation types, leading to range expansions and/or isolation in refugia. In contrast, comparatively fewer studies have addressed the continent's environmental heterogeneity and the role played by its geomorphological barriers. In this study, we performed a range-wide analysis of complete mitogenomes of sable antelope (Hippotragus niger) to explore how these different factors may have contributed as drivers of evolution in southcentral Africa. Our results supported two sympatric and deeply divergent mitochondrial lineages in west Tanzanian sables, which can be explained as the result of introgressive hybridization of a mitochondrial ghost lineage from an archaic, as-yet-undefined, congener. Phylogeographical subdivisions into three main lineages suggest that sable diversification may not have been driven solely by climatic events affecting populations differently across a continental scale. Often in interplay with climate, geomorphological features have also clearly shaped the species' patterns of vicariance, where the East Africa Rift System and the Eastern Arc Mountains acted as geological barriers. Subsequent splits among southern populations may be linked to rearrangements in the Zambezi system, possibly framing the most recent time when the river attained its current drainage profile. This work underlines how the use of comprehensive mitogenomic data sets on a model species with a wide geographical distribution can contribute to a much-enhanced understanding of environmental, geomorphological and evolutionary patterns in Africa throughout the Quaternary.

Demography and evolutionary history of grey wolf populations around the Bering Strait.

Glacial and interglacial periods throughout the Pleistocene have been substantial drivers of change in species distributions. Earlier analyses suggested that modern grey wolves (Canis lupus) trace their origin to a single Late Pleistocene Beringian population that expanded east and westwards, starting c. 25,000 years ago (ya). Here, we examined the demographic and phylogeographic histories of extant populations around the Bering Strait with wolves from two inland regions of the Russian Far East (RFE) and one coastal and two inland regions of North-western North America (NNA), genotyped for 91,327 single nucleotide polymorphisms. Our results indicated that RFE and NNA wolves had a common ancestry until c. 34,400 ya, suggesting that these populations started to diverge before the previously proposed expansion out of Beringia. Coastal and inland NNA populations diverged c. 16,000 ya, concordant with the minimum proposed date for the ecological viability of the migration route along the Pacific Northwest coast. Demographic reconstructions for inland RFE and NNA populations reveal spatial and temporal synchrony, with large historical effective population sizes that declined throughout the Pleistocene, possibly reflecting the influence of broadscale climatic changes across continents. In contrast, coastal NNA wolves displayed a consistently lower effective population size than the inland populations. Differences between the demographic history of inland and coastal wolves may have been driven by multiple ecological factors, including historical gene flow patterns, natural landscape fragmentation, and more recent anthropogenic disturbance.

Reliable wolf-dog hybrid detection in Europe using a reduced SNP panel developed for non-invasively collected samples.

BACKGROUND: Understanding the processes that lead to hybridization of wolves and dogs is of scientific and management importance, particularly over large geographical scales, as wolves can disperse great distances. However, a method to efficiently detect hybrids in routine wolf monitoring is lacking. Microsatellites offer only limited resolution due to the low number of markers showing distinctive allele frequencies between wolves and dogs. Moreover, calibration across laboratories is time-consuming and costly. In this study, we selected a panel of 96 ancestry informative markers for wolves and dogs, derived from the Illumina CanineHD Whole-Genome BeadChip (174 K). We designed very short amplicons for genotyping on a microfluidic array, thus making the method suitable also for non-invasively collected samples. RESULTS: Genotypes based on 93 SNPs from wolves sampled throughout Europe, purebred and non-pedigree dogs, and suspected hybrids showed that the new panel accurately identifies parental individuals, first-generation hybrids and first-generation backcrosses to wolves, while second- and third-generation backcrosses to wolves were identified as advanced hybrids in almost all cases. Our results support the hybrid identity of suspect individuals and the non-hybrid status of individuals regarded as wolves. We also show the adequacy of these markers to assess hybridization at a European-wide scale and the importance of including samples from reference populations. CONCLUSIONS: We showed that the proposed SNP panel is an efficient tool for detecting hybrids up to the third-generation backcrosses to wolves across Europe. Notably, the proposed genotyping method is suitable for a variety of samples, including non-invasive and museum samples, making this panel useful for wolf-dog hybrid assessments and wolf monitoring at both continental and different temporal scales.

Genomic evidence for the Old divergence of Southern European wolf populations.

The grey wolf (Canis lupus) is one of the most widely distributed mammals in which a variety of distinct populations have been described. However, given their currently fragmented distribution and recent history of human-induced population decline, little is known about the events that led to their differentiation. Based on the analysis of whole canid genomes, we examined the divergence times between Southern European wolf populations and their ancient demographic history. We found that all present-day Eurasian wolves share a common ancestor ca 36 000 years ago, supporting the hypothesis that all extant wolves derive from a single population that subsequently expanded after the Last Glacial Maximum. We also estimated that the currently isolated European populations of the Iberian Peninsula, Italy and the Dinarics-Balkans diverged very closely in time, ca 10 500 years ago, and maintained negligible gene flow ever since. This indicates that the current genetic and morphological distinctiveness of Iberian and Italian wolves can be attributed to their isolation dating back to the end of the Pleistocene, predating the recent human-induced extinction of wolves in Central Europe by several millennia.

Demographic expansion of an African opportunistic carnivore during the Neolithic revolution.

The diffusion of Neolithic technology together with the Holocene Climatic Optimum fostered the spread of human settlements and pastoral activities in North Africa, resulting in profound and enduring consequences for the dynamics of species, communities and landscapes. Here, we investigate the demographic history of the African wolf (Canis lupaster), a recently recognized canid species, to understand if demographic trends of this generalist and opportunistic carnivore reflect the increase in food availability that emerged after the arrival of the Neolithic economy in North Africa. We screened nuclear and mitochondrial DNA in samples collected throughout Algeria and Tunisia, and implemented coalescent approaches to estimate the variation of effective population sizes from present to ancestral time. We have found consistent evidence supporting the hypothesis that the African wolf population experienced a meaningful expansion concurring with a period of rapid population expansion of domesticates linked to the advent of agricultural practices.

Alternated selection mechanisms maintain adaptive diversity in different demographic scenarios of a large carnivore.

BACKGROUND: Different population trajectories are expected to impact the signature of neutral and adaptive processes at multiple levels, challenging the assessment of the relative roles of different microevolutionary forces. Here, we integrate adaptive and neutral variability patterns to disentangle how adaptive diversity is driven under different demographic scenarios within the Iberian wolf (Canis lupus) range. We studied the persistent, the expanding and a small, isolated group within the Iberian wolf population, using 3 MHC class II genes (DRB1, DQA1, and DQB1), which diversity was compared with 39 microsatellite loci. RESULTS: Both the persistent and the expanding groups show evidence of balancing selection, revealed by a significant departure from neutrality at MHC loci, significant higher observed and expected heterozygosity and lower differentiation at MHC than at neutral loci, and signs of positive selection. However, despite exhibiting a significantly higher genetic diversity than the isolated group, the persistent group did not show significant excess of MHC heterozygotes. The expanding group, while showing a similar level of genetic diversity than the persistent group, displays by contrast a significant excess of MHC heterozygotes, which is compatible with the heterozygote advantage mechanism. Results are not clear regarding the role of drift and selection in the isolated group due to the small size of this population. Although diversity indices of MHC loci correspond to neutral expectations in the isolated group, accelerated MHC divergence, revealed by a higher differentiation at MHC than neutral loci, may indicate diversifying selection. CONCLUSION: Different selective pressures were observed in the three different demographic scenarios, which are possibly driven by different selection mechanisms to maintain adaptive diversity.

Worldwide patterns of genomic variation and admixture in gray wolves.

The gray wolf (Canis lupus) is a widely distributed top predator and ancestor of the domestic dog. To address questions about wolf relationships to each other and dogs, we assembled and analyzed a data set of 34 canine genomes. The divergence between New and Old World wolves is the earliest branching event and is followed by the divergence of Old World wolves and dogs, confirming that the dog was domesticated in the Old World. However, no single wolf population is more closely related to dogs, supporting the hypothesis that dogs were derived from an extinct wolf population. All extant wolves have a surprisingly recent common ancestry and experienced a dramatic population decline beginning at least ∼30 thousand years ago (kya). We suggest this crisis was related to the colonization of Eurasia by modern human hunter-gatherers, who competed with wolves for limited prey but also domesticated them, leading to a compensatory population expansion of dogs. We found extensive admixture between dogs and wolves, with up to 25% of Eurasian wolf genomes showing signs of dog ancestry. Dogs have influenced the recent history of wolves through admixture and vice versa, potentially enhancing adaptation. Simple scenarios of dog domestication are confounded by admixture, and studies that do not take admixture into account with specific demographic models are problematic.

Unravelling population processes over the Late Pleistocene driving contemporary genetic divergence in Palearctic buzzards.

Population range expansions and contractions as a response to climate and habitat change throughout the Quaternary are known to have contributed to complex phylogenetic and population genetic events. Speciation patterns and processes in Palearctic buzzards (genus Buteo) are a long-standing example of morphological and genetic data incongruence, attributed to panmixia, habitat range shifts, contact zones, and climate change. Here we assess the systematics, phylogeography and population genetic structure of three nominal species of Palearctic buzzards, Buteo buteo (including B. b. vulpinus), B. rufinus (including B. r. cirtensis) and B. hemilasius. Phylogenetic analyses inferred from mitochondrial data recover B. hemilasius as sister to the sister clades B. r. rufinus and B. buteo complex (B. b. buteo, B. b. vulpinus, but also including B. r. cirtensis). In contrast, we find an unresolved genetic delimitation inferred from four nuclear loci, suggesting an ancestral genetic pool for all species. Time-trees suggest population contractions and expansions throughout the Pleistocene, which likely reflect habitat change and contrasting ecological niche requirements between species. Microsatellite-based extended Bayesian skyline plots reveal relatively constant population sizes for B. hemilasius, B. r. rufinus, and B. b. vulpinus, in contrast to a dramatic population expansion in B. r. cirtensis within the last 3 kya. Overall, our study illustrates how complex population processes over the Late Pleistocene have shaped the patterns of genetic divergence in Palearctic buzzards, due to the joint effects of shared ancestral polymorphisms, population expansions and contractions, with hybridization at contact zones leading to admixture and introgression.

Reassortment among picobirnaviruses found in wolves.

We conducted a viral metagenomics study in diarrheic free-ranging wolves in Portugal, revealing for the first time the presence of reassortant picobirnaviruses. These viruses shared identical capsid segments together with diverse RNA-dependent RNA polymerase segments. Even though causality between these picobirnaviruses and diarrhea could not be established, the study nonetheless confirms for the first time that wolves are a potential reservoir for picobirnaviruses, which might play a role as enteric pathogens.

Hybridization at an ecotone: ecological and genetic barriers between three Iberian vipers.

The formation of stable genetic boundaries between emerging species is often diagnosed by reduced hybrid fitness relative to parental taxa. This reduced fitness can arise from endogenous and/or exogenous barriers to gene flow. Although detecting exogenous barriers in nature is difficult, we can estimate the role of ecological divergence in driving species boundaries by integrating molecular and ecological niche modelling tools. Here, we focus on a three-way secondary contact zone between three viper species (Vipera aspis, V. latastei and V. seoanei) to test for the contribution of ecological divergence to the development of reproductive barriers at several species traits (morphology, nuclear DNA and mitochondrial DNA). Both the nuclear and mitochondrial data show that all taxa are genetically distinct and that the sister species V. aspis and V. latastei hybridize frequently and backcross over several generations. We find that the three taxa have diverged ecologically and meet at a hybrid zone coincident with a steep ecotone between the Atlantic and Mediterranean biogeographical provinces. Integrating landscape and genetic approaches, we show that hybridization is spatially restricted to habitats that are suboptimal for parental taxa. Together, these results suggest that niche separation and adaptation to an ecological gradient confer an important barrier to gene flow among taxa that have not achieved complete reproductive isolation.

Low signs of territorial behavior in the Eurasian otter during low-water conditions in a Mediterranean river.

The Eurasian otter Lutra lutra is a territorial semi-aquatic carnivore usually found at low densities in rivers, coastal areas, and wetlands. Its diet is based on prey associated with aquatic environments. Mediterranean rivers are highly seasonal, and suffer reduced flow during the summer, resulting in isolated river sections (pools) that sometimes can be left with a minimal amount of water, leading to concentrations of food for otters. To our knowledge, this process, which was known to field naturalists, has not been accurately described, nor have otter densities been estimated under these conditions. In this study, we describe the population size and movements of an aggregation of otters in an isolated pool in the Guadiana River in the Tablas de Daimiel National Park (central Spain), which progressively dried out during the spring-summer of 2022, in a context of low connectivity due to the absence of circulating water in the Guadiana and Gigüela rivers. Using non-invasive genetic sampling of 120 spraints collected along 79.4 km of sampling transects and spatial capture-recapture methods, we estimated the otter density at 1.71 individuals/km of river channel length (4.21 individuals/km2) in a progressively drying river pool, up to five times higher than previously described in the Iberian Peninsula. The movement patterns obtained with the spatial capture-recapture model are not quite different from those described in low density, which seems to indicate a wide home range overlap, with low signs of territoriality.

North African fox genomes show signatures of repeated introgression and adaptation to life in deserts.

Elucidating the evolutionary process of animal adaptation to deserts is key to understanding adaptive responses to climate change. Here we generated 82 individual whole genomes of four fox species (genus Vulpes) inhabiting the Sahara Desert at different evolutionary times. We show that adaptation of new colonizing species to a hot arid environment has probably been facilitated by introgression and trans-species polymorphisms shared with older desert resident species, including a putatively adaptive 25 Mb genomic region. Scans for signatures of selection implicated genes affecting temperature perception, non-renal water loss and heat production in the recent adaptation of North African red foxes (Vulpes vulpes), after divergence from Eurasian populations approximately 78 thousand years ago. In the extreme desert specialists, Rueppell's fox (V. rueppellii) and fennec (V. zerda), we identified repeated signatures of selection in genes affecting renal water homeostasis supported by gene expression and physiological differences. Our study provides insights into the mechanisms and genetic underpinnings of a natural experiment of repeated adaptation to extreme conditions.

The extinct Sicilian wolf shows a complex history of isolation and admixture with ancient dogs.

The Sicilian wolf remained isolated in Sicily from the end of the Pleistocene until its extermination in the 1930s-1960s. Given its long-term isolation on the island and distinctive morphology, the genetic origin of the Sicilian wolf remains debated. We sequenced four nuclear genomes and five mitogenomes from the seven existing museum specimens to investigate the Sicilian wolf ancestry, relationships with extant and extinct wolves and dogs, and diversity. Our results show that the Sicilian wolf is most closely related to the Italian wolf but carries ancestry from a lineage related to European Eneolithic and Bronze Age dogs. The average nucleotide diversity of the Sicilian wolf was half of the Italian wolf, with 37-50% of its genome contained in runs of homozygosity. Overall, we show that, by the time it went extinct, the Sicilian wolf had high inbreeding and low-genetic diversity, consistent with a population in an insular environment.

Genetic variability in mitochondrial and nuclear genes of Larus dominicanus (Charadriiformes, Laridae) from the Brazilian coast.

Several phylogeographic studies of seabirds have documented low genetic diversity that has been attributed to bottleneck events or individual capacity for dispersal. Few studies have been done in seabirds on the Brazilian coast and all have shown low genetic differentiation on a wide geographic scale. The Kelp Gull is a common species with a wide distribution in the Southern Hemisphere. In this study, we used mitochondrial and nuclear markers to examine the genetic variability of Kelp Gull populations on the Brazilian coast and compared this variability with that of sub-Antarctic island populations of this species. Kelp Gulls showed extremely low genetic variability for mitochondrial markers (cytb and ATPase) and high diversity for a nuclear locus (intron 7 of the ß-fibrinogen). The intraspecific evolutionary history of Kelp Gulls showed that the variability found in intron 7 of the ß-fibrinogen gene was compatible with the variability expected under neutral evolution but suggested an increase in population size during the last 10,000 years. However, none of the markers revealed evidence of a bottleneck population. These findings indicate that the recent origin of Kelp Gulls is the main explanation for their nuclear diversity, although selective pressure on the mtDNA of this species cannot be discarded.

On taming the effect of transcript level intra-condition count variation during differential expression analysis: A story of dogs, foxes and wolves.

The evolution of RNA-seq technologies has yielded datasets of scientific value that are often generated as condition associated biological replicates within expression studies. With expanding data archives opportunity arises to augment replicate numbers when conditions of interest overlap. Despite correction procedures for estimating transcript abundance, a source of ambiguity is transcript level intra-condition count variation; as indicated by disjointed results between analysis tools. We present TVscript, a tool that removes reference-based transcripts associated with intra-condition count variation above specified thresholds and we explore the effects of such variation on differential expression analysis. Initially iterative differential expression analysis involving simulated counts, where levels of intra-condition variation and sets of over represented transcripts are explicitly specified, was performed. Then counts derived from inter- and intra-study data representing brain samples of dogs, wolves and foxes (wolves vs. dogs and aggressive vs. tame foxes) were used. For simulations, the sensitivity in detecting differentially expressed transcripts increased after removing hyper-variable transcripts, although at levels of intra-condition variation above 5% detection became unreliable. For real data, prior to applying TVscript, ≈20% of the transcripts identified as being differentially expressed were associated with high levels of intra-condition variation, an over representation relative to the reference set. As transcripts harbouring such variation were removed pre-analysis, a discordance from 26 to 40% in the lists of differentially expressed transcripts is observed when compared to those obtained using the non-filtered reference. The removal of transcripts possessing intra-condition variation values within (and above) the 97th and 95th percentiles, for wolves vs. dogs and aggressive vs. tame foxes, maximized the sensitivity in detecting differentially expressed transcripts as a result of alterations within gene-wise dispersion estimates. Through analysis of our real data the support for seven genes with potential for being involved with selection for tameness is provided. TVscript is available at: https://sourceforge.net/projects/tvscript/.

Molecular tracking and prevalence of the red colour morph restricted to a harvested leopard population in South Africa.

The red leopard (Panthera pardus) colour morph is a colour variant that occurs only in South Africa, where it is confined to the Central Bushveld bioregion. Red leopards have been spreading over the past 40 years, which raises the speculation that the prevalence of this phenotype is related to low dispersal of young individuals owing to high off-take in the region. Intensive selective hunting tends to remove large resident male leopards from the breeding population, which gives young male leopards the chance to mate with resident female leopards that are more likely to be their relatives, eventually increasing the frequency of rare genetic variants. To investigate the genetic mechanisms underlying the red coat colour morph in leopards, and whether its prevalence in South Africa relates to an increase in genetic relatedness in the population, we sequenced exons of six coat colour-associated genes and 20 microsatellite loci in twenty Wild-type and four red leopards. The results were combined with demographic data available from our study sites. We found that red leopards own a haplotype in homozygosity identified by two SNPs and a 1 bp deletion that causes a frameshift in the tyrosinase-related protein 1 (TYRP1), a gene known to be involved in the biosynthesis of melanin. Microsatellite analyses indicate clear signs of a population bottleneck and a relatedness of 0.11 among all pairwise relationships, eventually supporting our hypothesis that a rare colour morph in the wild has increased its local frequency due to low natal dispersal, while subject to high human-induced mortality rate.

The Cantabrian capercaillie: A population on the edge.

The capercaillie Tetrao urogallus - the world's largest grouse- is a circumboreal forest species, which only two remaining populations in Spain: one in the Cantabrian mountains in the west and the other in the Pyrenees further east. Both have shown severe declines, especially in the Cantabrian population, which has recently been classified as "Critically Endangered". To develop management plans, information on demographic parameters is necessary to understand and forecast population dynamics. We used spatial capture-recapture (SCR) modeling and non-invasive DNA samples to estimate the current population size in the whole Cantabrian mountain range. In addition, for the assessment of population status, we analyzed the population trajectory over the last 42 years (1978-2019) at 196 leks on the Southern slope of the range, using an integrated population model with a Dail-Madsen model at its core, combined with a multistate capture-recapture model for survival and a Poisson regression for productivity. For 2019, we estimate the size of the entire population at 191 individuals (95% BCI 165-222) for an estimated 60 (48-78) females and 131 (109-157) males. Since the 1970s, our study estimates a shrinkage of the population range by 83%. The population at the studied leks in 2019 was at about 10% of the size estimated for 1978. Apparent annual survival was estimated at 0.707 (0.677-0.735), and per-capita recruitment at 0.233 (0.207-0.262), and insufficient to maintain a stable population. We suggest work to improve the recruitment (and survival) and manage these mountain forests for capercaillie conservation. Also, in the future, management should assess the genetic viability of this population.