Genome‐wide association study identifies genomic regions of association for cruciate ligament rupture in Newfoundland dogs

Cranial cruciate ligament rupture (CCLR) is the most common cause of pelvic limb lameness in dogs. To investigate the genetic basis of canine CCLR, we conducted a genome‐wide association study using a canine SNP array in Newfoundland pedigree dogs with and without CCLR (n = 96). We identified three main chromosomal regions of CCLR association (on chromosomes 1, 3 and 33). Each of these regions was confirmed by Sequenom genotyping in a further cohort of Newfoundlands (n = 271). The results, particularly SNPs identified in the SORCS2 and SEMA5B genes, suggest that there may be neurological pathways involved in susceptibility to canine CCLR.     

Sequence analysis of three canine adipokine genes revealed an association between TNF polymorphisms and obesity in Labrador dogs

Obesity is an emerging health problem in purebred dogs. Due to their crucial role in energy homeostasis control, genes encoding adipokines are considered candidate genes, and their variants may be associated with predisposition to obesity. Searching for polymorphism was carried out in three adipokine genes (TNF, RETN and IL6). The study was performed on 260 dogs, including lean (n = 109), overweight (n = 88) and obese (n = 63) dogs. The largest cohort was represented by Labrador Retrievers (n = 136). Altogether, 24 novel polymorphisms were identified: 12 in TNF (including one missense SNP), eight in RETN (including one missense SNP) and four in IL6. Distributions of five common SNPs (two in TNF, two in RETN and one in IL6) were further analyzed with regard to body condition score. Two SNPs in the non‐coding parts of TNF (c.‐40A>C and c.233+14G>A) were associated with obesity in Labrador dogs. The obtained results showed that the studied adipokine genes are highly polymorphic and two polymorphisms in the TNF gene may be considered as markers predisposing Labrador dogs to obesity.     

Development and evaluation of a genome‐wide Coffee 8.5K SNP array and its application for high‐density genetic mapping and for investigating the origin of Coffea arabica L.

Coffee species such as Coffea canephora P. (Robusta) and C. arabica L. (Arabica) are important cash crops in tropical regions around the world. C. arabica is an allotetraploid (2n = 4x = 44) originating from a hybridization event of the two diploid species C. canephora and C. eugenioides (2n = 2x = 22). Interestingly, these progenitor species harbour a greater level of genetic variability and are an important source of genes to broaden the narrow Arabica genetic base. Here, we describe the development, evaluation and use of a single‐nucleotide polymorphism (SNP) array for coffee trees. A total of 8580 unique and informative SNPs were selected from C. canephora and C. arabica sequencing data, with 40% of the SNP located in annotated genes. In particular, this array contains 227 markers associated to 149 genes and traits of agronomic importance. Among these, 7065 SNPs (~82.3%) were scorable and evenly distributed over the genome with a mean distance of 54.4 Kb between markers. With this array, we improved the Robusta high‐density genetic map by adding 1307 SNP markers, whereas 945 SNPs were found segregating in the Arabica mapping progeny. A panel of C. canephora accessions was successfully discriminated and over 70% of the SNP markers were transferable across the three species. Furthermore, the canephora‐derived subgenome of C. arabica was shown to be more closely related to C. canephora accessions from northern Uganda than to other current populations. These validated SNP markers and high‐density genetic maps will be useful to molecular genetics and for innovative approaches in coffee breeding.     

SNP‐based heritability and genetic architecture of tarsal osteochondrosis in North American Standardbred horses

Osteochondrosis is a common developmental orthopedic disease characterized by a failure of endochondral ossification. Standardbred horses are recognized as being predisposed to tarsal osteochondrosis. Prior heritability estimates for tarsal osteochondrosis in European Standardbreds and related trotting breeds have been based on pedigree data and range from 17–29%. Here, we report on genetic architecture and heritability based on high‐density genotyping data in a cohort of North American Standardbreds (n = 479) stringently phenotyped for tarsal osteochondrosis. Whole‐genome array genotyping data were imputed to ~2 million single nucleotide polymorphisms (SNPs). SNP‐based heritability of osteochondrosis in this population was explained by 2326 SNPs. The majority of these SNPs (86.6%) had small effects, whereas fewer SNPs had moderate or large effects (10% and 2.9% respectively), which is consistent with a polygenic/complex disease. Heritability was estimated at 0.24 ± 0.16 using two methods of restricted maximum likelihood analysis, as implemented in gcta (with and without a weighted relatedness matrix) and ldak software. Estimates were validated using bootstrapping. Heritability estimates were within the range previously reported and suggest that osteochondrosis is moderately heritable but that a significant portion of disease risk is due to environmental factors and/or genotype × environment interactions. Future identification of the genes/variants that have the most impact on disease risk may allow early recognition of high‐risk individuals.     

Evaluation of a SNP map of 6q24–27 confirms diabetic nephropathy loci and identifies novel associations in type 2 diabetes patients with nephropathy from an African-American population

Previously, we performed a genome scan for type 2 diabetes (T2DM) using 638 African-American (AA) affected sibling pairs from 247 families; non-parametric linkage analysis suggested evidence of linkage at 6q24–27 (LOD 2.26). To comprehensively evaluate this region, we performed a two-stage association study by first constructing a SNP map of 754 SNPs selected from HapMap on the basis of linkage disequilibrium (LD) in 300 AAT2DM end-stage renal disease (ESRD) subjects, 311 AA controls, 43 European American controls and 45 Yoruba Nigerian samples (Set 1). Replication analyses were conducted in an independent population of 283 AA T2DM-ESRD subjects and 282 AA controls (Set 2). In addition, we adjusted for the impact of admixture on association results by using ancestry informative markers (AIMs). In Stage 1, 137 (18.2%) SNPs showed nominal evidence of association (P < 0.05) in one or more of tests of association: allelic (n = 33), dominant (n = 36), additive (n = 29), or recessive (n = 34) genotypic models, and 2- (n = 47) and 3-SNP (n = 43) haplotypic analyses. These SNPs were selected for follow-up genotyping. Stage 2 analyses confirmed association with a predicted 2-SNP “risk” haplotype in the PARK2 gene. Also, two intergenic SNPs showed consistent genotypic association with T2DM-ESRD: rs12197043 and rs4897081. Combined analysis of all subjects from both stages revealed nominal associations with 17 SNPs within genes, including suggestive associations in ESR1 and PARK2. This study confirms known diabetic nephropathy loci and identifies potentially novel susceptibility variants located within 6q24–27 in AA. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Atopic dermatitis in West Highland white terriers is associated with a 1.3-Mb region on CFA 17

Canine atopic dermatitis (AD) is an allergic inflammatory skin disease that shares similarities with AD in humans. Canine AD is likely to be an inherited disease in dogs and is common in West Highland white terriers (WHWTs). We performed a genome-wide association study using the Affymetrix Canine SNP V2 array consisting of over 42,800 single nucleotide polymorphisms, on 35 atopic and 25 non-atopic WHWTs. A gene-dropping simulation method, using SIB-PAIR, identified a projected 1.3 Mb area of association (genome-wide P = 6 × 10⁻⁵ to P = 7 × 10⁻⁴) on CFA 17. Nineteen genes on CFA 17, including 1 potential candidate gene (PTPN22), were located less than 0.5 Mb from the interval of association identified on the genome-wide association analysis. Four haplotypes within this locus were differently distributed between cases and controls in this population of dogs. These findings suggest that a major locus for canine AD in WHWTs may be located on, or in close proximity to an area on CFA 17.     

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.     

Expression of TGFβ1 in pulmonary vein stenosis after radiofrequency ablation in chronic atrial fibrillation of dogs

The development of pulmonary vein stenosis has recently been described after radiofrequency ablation (RF) to treat atrial fibrillation (AF). The purpose of this study was to examine expression of TGFβ1 in pulmonary vein stenosis after radiofrequency ablation in chronic atrial fibrillation of dogs. About 28 mongrel dogs were randomly assigned to the sham-operated group (n = 7), the AF group (n = 7), AF + RF group (n = 7), and RF group (n = 7). In AF or AF + RF groups, dogs underwent chronic pulmonary vein (PV) pacing to induce sustained AF. RF application was applied around the PVs until electrical activity was eliminated. Histological assessment of pulmonary veins was performed using hematoxylin and eosin staining; TGFβ1 gene expression in pulmonary veins was examined by RT-PCR analysis; expression of TGFβ1 protein in pulmonary veins was assessed by Western blot analysis. Rapid pacing from the left superior pulmonary vein (LSPV) induced sustained AF in AF group and AF + RF group. Pulmonary vein ablation terminated the chronic atrial fibrillation in dogs. Histological examination revealed necrotic tissues in various stages of collagen replacement, intimal thickening, and cartilaginous metaplasia with chondroblasts and chondroclasts. Compared with sham-operated and AF group, TGFβ1 gene and protein expressions was increased in AF + RF or RF groups. It was concluded that TGFβ1 might be associated with pulmonary vein stenosis after radiofrequency ablation in chronic atrial fibrillation of dogs. Shufeng Li and Hongli Li contributed equally to the work.     

Serological Evidence of Infection by Paracoccidioides brasiliensis in Dogs with Leishmaniasis

Paracoccidioidomycosis is a systemic mycosis prevalent in Latin American countries, caused by the dimorphic fungi Paracoccidioides brasiliensis and P. lutzii. The habitat of these fungi in nature remains undefined, although it is believed that infection occurs by inhalation of infective propagules present in soil. Sentinel animals, such as dogs, can be valuable epidemiological markers of paracoccidioidomycosis. Taking into account that paracoccidioidomycosis and visceral leishmaniasis may occur in the same area, the objective of this study was to evaluate the occurrence of P. brasiliensis infection in dogs positive for Leishmania sp. Serum samples of dogs positive (n = 199) and negative (n = 101) for Leishmania sp. were analyzed by the immunodiffusion test using P. brasiliensis exoantigen, and 22 samples (7.3%) were positive. The serum samples positive in the immunodiffusion test were also analyzed by Western blotting using the P. brasiliensis gp43 recombinant protein, and 86% of the samples were positive. A high positive correlation (r = 0.96) between positivity for Leishmania sp. and P. brasiliensis was observed. These data suggest an association between leishmaniasis and paracoccidioidomycosis in dogs.

     

Viral gut metagenomics of sympatric wild and domestic canids,and monitoring of viruses: Insights from an endangered wolf population

Animal host–microbe interactions are a relevant concern for wildlife conservation, particularly regarding generalist pathogens, where domestic host species can play a role in the transmission of infectious agents, such as viruses, to wild animals. Knowledge on viral circulation in wild host species is still scarce and can be improved by the recent advent of modern molecular approaches. We aimed to characterize the fecal virome and identify viruses of potential conservation relevance of diarrheic free‐ranging wolves and sympatric domestic dogs from Central Portugal, where a small and threatened wolf population persists in a highly anthropogenically modified landscape. Using viral metagenomics, we screened diarrheic stools collected from wolves (n = 8), feral dogs (n = 4), and pet dogs (n = 6), all collected within wolf range. We detected novel highly divergent viruses as well as known viral pathogens with established effects on population dynamics, including canine distemper virus, a novel bocavirus, and canine minute virus. Furthermore, we performed a 4‐year survey for the six wolf packs comprising this endangered wolf population, screening 93 fecal samples from 36 genetically identified wolves for canine distemper virus and the novel bocavirus, previously identified using our metagenomics approach. Our novel approach using metagenomics for viral screening in noninvasive samples of wolves and dogs has profound implications on the knowledge of both virology and wildlife diseases, establishing a complementary tool to traditional screening methods for the conservation of threatened species.     

Dietary Enrichment with Medium Chain Triglycerides (AC-1203) Elevates Polyunsaturated Fatty Acids in the Parietal Cortex of Aged Dogs: Implications for Treating Age-Related Cognitive Decline

Dogs demonstrate an age-related cognitive decline, which may be related to a decrease in the concentration of omega-3 polyunsaturated fatty acids (n-3 PUFA) in the brain. Medium chain triglycerides (MCT) increase fatty acid oxidation, and it has been suggested that this may raise brain n-3 PUFA levels by increasing mobilization of n-3 PUFA from adipose tissue to the brain. The goal of the present study was to determine whether dietary MCT would raise n-3 PUFA concentrations in the brains of aged dogs. Eight Beagle dogs were randomized to a control diet (n = 4) or an MCT (AC-1203) enriched diet (n = 4) for 2 months. The animals were then euthanized and the parietal cortex was removed for phospholipid, cholesterol and fatty acid determinations by gas-chromatography. Dietary enrichment with MCT (AC-1203) resulted in a significant increase in brain phospholipid and total lipid concentrations (P < 0.05). In particular, n-3 PUFA within the phospholipid, unesterified fatty acid, and total lipid fractions were elevated in AC-1203 treated subjects as compared to controls (P < 0.05). Brain cholesterol concentrations did not differ significantly between the groups (P > 0.05). These results indicate that dietary enrichment with MCT, raises n-3 PUFA concentrations in the parietal cortex of aged dogs.     

Genomic regions affecting backfat thickness and cannon bone circumference identified by genome‐wide association study in a Duroc pig population

We performed a genome‐wide association study using the porcine 60K SNP array to detect QTL regions for nine traits in a three‐generational Duroc samples (n = 651), viz. generations 1, 2 and 3 from a population selected over five generations using a closed nucleus breeding scheme. We applied a linear mixed model for association mapping to detect SNP effects, adjusting for fixed effects (sex and season) and random polygenic effects (reflecting genetic relatedness), and derived a likelihood ratio statistic for each SNP using the efficient mixed‐model association method. We detected a region on SSC6 for backfat thickness (BFT) and on SSC7 for cannon bone circumference (CANNON), with a genome‐wide significance of P < 0.01 after Bonferroni correction. These regions had been detected previously in other pig populations. Six genes are located in the BFT‐associated region, while the CANNON‐associated region includes 66 genes. In the future, significantly associated SNPs, derived by sequencing the coding regions of the six genes in the BFT region, can be used in marker‐assisted selection of BFT, whereas haplotypes constructed from the SSC7 region with strong LD can be used to select for the CANNON trait in our resource family.     

Development of an integrated 200K SNP genotyping array and application for genetic mapping,genome assembly improvement and genome wide association studies in pear (Pyrus)

Pear (Pyrus; 2n = 34), the third most important temperate fruit crop, has great nutritional and economic value. Despite the availability of many genomic resources in pear, it is challenging to genotype novel germplasm resources and breeding progeny in a timely and cost‐effective manner. Genotyping arrays can provide fast, efficient and high‐throughput genetic characterization of diverse germplasm, genetic mapping and breeding populations. We present here 200K AXIOM^® PyrSNP, a large‐scale single nucleotide polymorphism (SNP) genotyping array to facilitate genotyping of Pyrus species. A diverse panel of 113 re‐sequenced pear genotypes was used to discover SNPs to promote increased adoption of the array. A set of 188 diverse accessions and an F₁ population of 98 individuals from ‘Cuiguan’ × ‘Starkrimson’ was genotyped with the array to assess its effectiveness. A large majority of SNPs (166 335 or 83%) are of high quality. The high density and uniform distribution of the array SNPs facilitated prediction of centromeric regions on 17 pear chromosomes, and significantly improved the genome assembly from 75.5% to 81.4% based on genetic mapping. Identification of a gene associated with flowering time and candidate genes linked to size of fruit core via genome wide association studies showed the usefulness of the array in pear genetic research. The newly developed high‐density SNP array presents an important tool for rapid and high‐throughput genotyping in pear for genetic map construction, QTL identification and genomic selection.     

A 34K SNP genotyping array for Populus trichocarpa: Design,application to the study of natural populations and transferability to other Populus species

Genetic mapping of quantitative traits requires genotypic data for large numbers of markers in many individuals. For such studies, the use of large single nucleotide polymorphism (SNP) genotyping arrays still offers the most cost‐effective solution. Herein we report on the design and performance of a SNP genotyping array for Populus trichocarpa (black cottonwood). This genotyping array was designed with SNPs pre‐ascertained in 34 wild accessions covering most of the species latitudinal range. We adopted a candidate gene approach to the array design that resulted in the selection of 34 131 SNPs, the majority of which are located in, or within 2 kb of, 3543 candidate genes. A subset of the SNPs on the array (539) was selected based on patterns of variation among the SNP discovery accessions. We show that more than 95% of the loci produce high quality genotypes and that the genotyping error rate for these is likely below 2%. We demonstrate that even among small numbers of samples (n = 10) from local populations over 84% of loci are polymorphic. We also tested the applicability of the array to other species in the genus and found that the number of polymorphic loci decreases rapidly with genetic distance, with the largest numbers detected in other species in section Tacamahaca. Finally, we provide evidence for the utility of the array to address evolutionary questions such as intraspecific studies of genetic differentiation, species assignment and the detection of natural hybrids.     

Identification of SNP markers for common CNV regions and association analysis of risk of subarachnoid aneurysmal hemorrhage in Japanese population

Copy number variation (CNV) is emerging as a new tool for understanding human genomic variation, but its relationship with human disease is not yet fully understood. The data for a total of 317,503 genotypes were collected for a genome-wide association study of subarachnoid aneurismal hemorrhage (SAH) in a Japanese population (cases and controls, n = 497) using Illumina HumanHap300 BeadChip^®. To identify multi-allelic CNV markers, we visually inspected all genotype clusters of 317,503 SNP markers covering the whole genome using Illumina’s BeadStudio 3.0^® software. As a result, we identified 597 multi-allelic CNV markers for common (copy loss frequency > 0.05) CNV regions in a Japanese population (n = 497). The identified CNV markers shared the following characteristics: enrichment of Hardy–Weinberg disequilibria, Mendelian inconsistency among families, and high missing genotype rate. All annotated information for those markers is summarized in our database (http://www.snp-genetics.com/user/srch.htm). In addition, we performed case-control association analyses of identified multi-allelic CNV markers with the risk of subarachnoid aneurysmal hemorrhage. One SNP marker (rs1242541) within a CNV region neighboring the Sel-1 suppressor of lin-12-like protein (SEL1L) was significantly associated with a risk of SAH (P = 0.0006). We also validated the CNV around rs1242541 using real-time quantitative polymerase chain reaction (PCR). Information and methods used in this study would be helpful for accurate genotyping of SNPs on CNV regions, which could be used for association analysis of SNP markers within CNV regions.     

Multi-drug resistant,biofilm-producing high-risk clonal lineage of Klebsiella in companion and household animals

Antimicrobial resistance is a global emergency which needs one health approach to address. The present study was conducted to detect the prevalence of beta-lactamase and biofilm-producing Klebsiella strains in rectal swabs (n = 624) collected from healthy dogs, cats, sheep and goats reared as companion or household animals in India. The dogs and cats were frequently exposed to third- or fourth-generation cephalosporins for therapy. The sheep and goats were occasionally exposed to antibiotics and had environmental exposure. Phenotypical ESBL (n = 93) and ACBL (n = 88)-producing Klebsiella were isolated significantly more (P < 0·05) from companion animals than household animals. Majority of the Klebsiella possessed bla_CTX-M-15. The sequences bla_CTX-M-15.2, bla_CTX-M-197 and bla_CTX-M-225 are reported first time from the companion animals. All ACBL-producing isolates possessed bla_AmpC. The present study detected 65·8% of Klebsiella strains as biofilm producers possessing the studied biofilm associated genes. The isolates showed phenotypical resistance against chloramphenicol, tetracycline, doxycycline, co-trimoxazole, ampicillin, cefotaxime/clavulanic acid. The present study showed that companion and household animals (dogs, cats, sheep, goats) may act as a carrier of ESBL/biofilm-producing, multi-drug resistant, high-risk clonal lineage of Klebsiella.     

Spatial dynamics and mixing of bluefin tuna in the Atlantic Ocean and Mediterranean Sea revealed using next‐generation sequencing

The Atlantic bluefin tuna is a highly migratory species emblematic of the challenges associated with shared fisheries management. In an effort to resolve the species’ stock dynamics, a genomewide search for spatially informative single nucleotide polymorphisms (SNPs) was undertaken, by way of sequencing reduced representation libraries. An allele frequency approach to SNP discovery was used, combining the data of 555 larvae and young‐of‐the‐year (LYOY) into pools representing major geographical areas and mapping against a newly assembled genomic reference. From a set of 184,895 candidate loci, 384 were selected for validation using 167 LYOY. A highly discriminatory genotyping panel of 95 SNPs was ultimately developed by selecting loci with the most pronounced differences between western Atlantic and Mediterranean Sea LYOY. The panel was evaluated by genotyping a different set of LYOY (n = 326), and from these, 77.8% and 82.1% were correctly assigned to western Atlantic and Mediterranean Sea origins, respectively. The panel revealed temporally persistent differentiation among LYOY from the western Atlantic and Mediterranean Sea (F_ST = 0.008, p = .034). The composition of six mixed feeding aggregations in the Atlantic Ocean and Mediterranean Sea was characterized using genotypes from medium (n = 184) and large (n = 48) adults, applying population assignment and mixture analyses. The results provide evidence of persistent population structuring across broad geographic areas and extensive mixing in the Atlantic Ocean, particularly in the mid‐Atlantic Bight and Gulf of St. Lawrence. The genomic reference and genotyping tools presented here constitute novel resources useful for future research and conservation efforts.     

SNPs in the bovine IL-10 receptor are associated with somatic cell score in Canadian dairy bulls

Altering the balance between pro- and anti-inflammatory responses can influence an animal’s susceptibility to acute or chronic inflammatory disease; bovine mastitis is no exception. Genetic variation in the form of single nucleotide polymorphisms (SNPs) may alter the function and expression of genes that regulate inflammation, making them important candidates for defining an animal’s risk of developing acute or chronic mastitis. The objective of the present study was to identify SNPs in genes that regulate anti-inflammatory responses and test their association with estimated breeding values (EBVs) for somatic cell score (SCS), a trait highly correlated with the incidence of mastitis. These genes included bovine interleukin-10 (IL-10) and its receptor (IL-10R), and transforming growth factor β1 (TGF-β1) and its receptor (TGF-βR). Sequencing-pooled DNA allowed for the identification of SNPs in IL-10 (n = 2), IL-10Rα (n = 6) and β (n = 2), and TGF-β1 (n = 1). These SNPs were subsequently genotyped in a cohort of Holstein (n = 500), Jersey (n = 83), and Guernsey (n = 50) bulls. Linear regression analysis identified significant SNP effects for IL-10Rα 1185C>T with SCS. Haplotype IL-10Rα AAT showed a significant effect on increasing SCS compared to the most common haplotype. The results presented here indicate that SNPs in IL-10Rα may contribute to variation in the SCS of dairy cattle. Although functional studies are necessary to ascertain whether these SNPs are causal polymorphisms or merely in linkage with the true causal SNP(s), a selection program incorporating these markers could have a beneficial influence on the average SCS and productivity of a dairy herd.     

Genome‐wide analyses suggest parallel selection for universal traits may eclipse local environmental selection in a highly mobile carnivore

Ecological and environmental heterogeneity can produce genetic differentiation in highly mobile species. Accordingly, local adaptation may be expected across comparatively short distances in the presence of marked environmental gradients. Within the European continent, wolves (Canis lupus) exhibit distinct north–south population differentiation. We investigated more than 67‐K single nucleotide polymorphism (SNP) loci for signatures of local adaptation in 59 unrelated wolves from four previously identified population clusters (northcentral Europe n = 32, Carpathian Mountains n = 7, Dinaric‐Balkan n = 9, Ukrainian Steppe n = 11). Our analyses combined identification of outlier loci with findings from genome‐wide association study of individual genomic profiles and 12 environmental variables. We identified 353 candidate SNP loci. We examined the SNP position and neighboring megabase (1 Mb, one million bases) regions in the dog (C. lupus familiaris) genome for genes potentially under selection, including homologue genes in other vertebrates. These regions included functional genes for, for example, temperature regulation that may indicate local adaptation and genes controlling for functions universally important for wolves, including olfaction, hearing, vision, and cognitive functions. We also observed strong outliers not associated with any of the investigated variables, which could suggest selective pressures associated with other unmeasured environmental variables and/or demographic factors. These patterns are further supported by the examination of spatial distributions of the SNPs associated with universally important traits, which typically show marked differences in allele frequencies among population clusters. Accordingly, parallel selection for features important to all wolves may eclipse local environmental selection and implies long‐term separation among population clusters.