Mitochondrial DNA and Y-chromosomal diversity in ancient populations of domestic sheep (Ovis aries) in Finland: comparison with contemporary sheep breeds

Background

Several molecular and population genetic studies have focused on the native sheep breeds of Finland. In this work, we investigated their ancestral sheep populations from Iron Age, Medieval and Post-Medieval periods by sequencing a partial mitochondrial DNA D-loop and the 5’-promoter region of the SRY gene. We compared the maternal (mitochondrial DNA haplotypes) and paternal (SNP oY1) genetic diversity of ancient sheep in Finland with modern domestic sheep populations in Europe and Asia to study temporal changes in genetic variation and affinities between ancient and modern populations.

Results

A 523-bp mitochondrial DNA sequence was successfully amplified for 26 of 36 sheep ancient samples i.e. five, seven and 14 samples representative of Iron Age, Medieval and Post-Medieval sheep, respectively. Genetic diversity was analyzed within the cohorts. This ancient dataset was compared with present-day data consisting of 94 animals from 10 contemporary European breeds and with GenBank DNA sequence data to carry out a haplotype sharing analysis. Among the 18 ancient mitochondrial DNA haplotypes identified, 14 were present in the modern breeds. Ancient haplotypes were assigned to the highly divergent ovine haplogroups A and B, haplogroup B being the major lineage within the cohorts. Only two haplotypes were detected in the Iron Age samples, while the genetic diversity of the Medieval and Post-Medieval cohorts was higher. For three of the ancient DNA samples, Y-chromosome SRY gene sequences were amplified indicating that they originated from rams. The SRY gene of these three ancient ram samples contained SNP G-oY1, which is frequent in modern north-European sheep breeds.

Conclusions

Our study did not reveal any sign of major population replacement of native sheep in Finland since the Iron Age. Variations in the availability of archaeological remains may explain differences in genetic diversity estimates and patterns within the cohorts rather than demographic events that occurred in the past. Our ancient DNA results fit well with the genetic context of domestic sheep as determined by analyses of modern north-European sheep breeds.     

Mitochondrial DNA sequence variation and phylogenetic relationships among Iberian pigs and other domestic and wild pig populations

Nucleotide sequences of mitochondrial DNA (mtDNA) cytochrome B gene (1140 bp) and control region (707 bp) were used to determine the phylogenetic relationships among 51 pig samples representing ancient and current varieties of Iberian pigs (26), Spanish wild boars (seven) and other domestic pigs (18) of cosmopolitan (Duroc, Large White, Landrace, Pietrain and Meishan) and local (Spotted Black Jabugo, Basque and Mangalitza) breeds. A neighbour-joining tree constructed from pairwise distances provide evidence of the European origin of both Iberian pigs and Spanish wild boars. The introgression of Asian mtDNA haplotypes in the genetic pool of the Iberian breed seems unlikely. Four estimates of sequence divergence between European and Asian clades were calculated from the two main domains of the D-loop region and the synonymous and nonsynonymous nucleotide substitutions in the cytochrome B gene. The time since the divergence of pig ancestors was estimated at about 600,000 years before present.     

Genome sequencing and analysis of Mangalica,a fatty local pig of Hungary

Background

Mangalicas are fatty type local/rare pig breeds with an increasing presence in the niche pork market in Hungary and in other countries. To explore their genetic resources, we have analysed data from next-generation sequencing of an individual male from each of three Mangalica breeds along with a local male Duroc pig. Structural variations, such as SNPs, INDELs and CNVs, were identified and particular genes with SNP variations were analysed with special emphasis on functions related to fat metabolism in pigs.

Results

More than 60 Gb of sequence data were generated for each of the sequenced individuals, resulting in 11× to 19× autosomal median coverage. After stringent filtering, around six million SNPs, of which approximately 10% are novel compared to the dbSNP138 database, were identified in each animal. Several hundred thousands of INDELs and about 1,000 CNV gains were also identified. The functional annotation of genes with exonic, non-synonymous SNPs, which are common in all three Mangalicas but are absent in either the reference genome or the sequenced Duroc of this study, highlighted 52 genes in lipid metabolism processes. Further analysis revealed that 41 of these genes are associated with lipid metabolic or regulatory pathways, 49 are in fat-metabolism and fatness-phenotype QTLs and, with the exception of ACACA, ANKRD23, GM2A, KIT, MOGAT2, MTTP, FASN, SGMS1, SLC27A6 and RETSAT, have not previously been associated with fat-related phenotypes.

Conclusions

Genome analysis of Mangalica breeds revealed that local/rare breeds could be a rich source of sequence variations not present in cosmopolitan/industrial breeds. The identified Mangalica variations may, therefore, be a very useful resource for future studies of agronomically important traits in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-761) contains supplementary material, which is available to authorized users.     

Ancient DNA sheds light on the ancestry of pre-hispanic Canarian pigs

Background

Canarian Black (CB) pigs belong to an autochthonous and endangered breed, which is spread throughout the Canarian archipelago. It is commonly accepted that they represent a relic of the pig populations that were bred by the Berbers in North Africa over millennia. It is important to note that the geographic isolation of the Canary Islands has preserved this genetic legacy intact from foreign introgressions until the Spanish conquest of the archipelago in the 15^th century. Ten years ago, it was demonstrated that, in CB pigs, the frequency of the Asian A2 cytochrome-b haplogroup reached 73%. The current work aimed at investigating whether this observation is explained by either a recent or an ancient introgression of CB pigs with Far Eastern pigs.

Results

Genetic analyses of 23 ancient samples from pre-hispanic Canarian pigs (420 to 2500 years before present) showed that Near Eastern and Far Eastern genetic signatures were totally absent in the primitive Canarian pre-hispanic pigs. Indeed, the haplotypes detected in these pigs were closely related to those of North African and European wild boars.

Conclusions

Our results demonstrate that the high frequency of the Far Eastern mitochondrial cytochrome B A2 haplotype in modern Canarian Black pigs probably corresponds to a relatively recent introgression with British breeds.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0115-7) contains supplementary material, which is available to authorized users.     

TYR基因外显子1的序列变异

为了分析家猪与野猪的遗传多样性及起源,测定了来自１２个中国地方家猪品种、３个欧洲引进猪品种以及8个中国野猪和2个越南野猪共36个个体的酪氨酸酶基因（TYR）外显子1的序列,共检出６个单核苷酸多态性位点（SNPs）,且这6个位点的变异均为同义突变,根据这些变异可将酪氨酸酶基因DNA序列归结为４种单倍型。结合已发表的数据,构建了简约中介网络图。 在网络图中,单倍型TYR*2主要为欧洲家猪与欧洲野猪和三条亚洲家猪染色体。大部分亚洲家猪和野猪共享单倍型TYR*1,表明这是一个亚洲类型的单倍型;同时也有部分欧洲家猪与野猪携带这一单倍型。 而单倍型TYR*3和TYR*4为本研究检测到的稀有单倍型,这两种单倍型主要由中国家猪与亚洲野猪组成。这种网络图结构支持家猪的欧洲和亚洲独立起源学说,同时也表明相当部分的欧洲家猪品种受到亚洲猪的基因渗透,而少量中国家猪和日本野猪也受到了欧洲猪的基因渗透。     

A signature inferred from Drosophila mitotic genes predicts survival of breast cancer patients

Introduction

The classification of breast cancer patients into risk groups provides a powerful tool for the identification of patients who will benefit from aggressive systemic therapy. The analysis of microarray data has generated several gene expression signatures that improve diagnosis and allow risk assessment. There is also evidence that cell proliferation-related genes have a high predictive power within these signatures.

Methods

We thus constructed a gene expression signature (the DM signature) using the human orthologues of 108 Drosophila melanogaster genes required for either the maintenance of chromosome integrity (36 genes) or mitotic division (72 genes).

Results

The DM signature has minimal overlap with the extant signatures and is highly predictive of survival in 5 large breast cancer datasets. In addition, we show that the DM signature outperforms many widely used breast cancer signatures in predictive power, and performs comparably to other proliferation-based signatures. For most genes of the DM signature, an increased expression is negatively correlated with patient survival. The genes that provide the highest contribution to the predictive power of the DM signature are those involved in cytokinesis.

Conclusion

This finding highlights cytokinesis as an important marker in breast cancer prognosis and as a possible target for antimitotic therapies.     

Genotyping of Capreolus pygargus fossil DNA from Denisova cave reveals phylogenetic relationships between ancient and modern populations

Background

The extant roe deer (Capreolus Gray, 1821) includes two species: the European roe deer (C. capreolus) and the Siberian roe deer (C. pygargus) that are distinguished by morphological and karyotypical differences. The Siberian roe deer occupies a vast area of Asia and is considerably less studied than the European roe deer. Modern systematics of the Siberian roe deer remain controversial with 4 morphological subspecies. Roe deer fossilized bones are quite abundant in Denisova cave (Altai Mountains, South Siberia), where dozens of both extant and extinct mammalian species from modern Holocene to Middle Pleistocene have been retrieved.

Methodology/Principal Findings

We analyzed a 629 bp fragment of the mitochondrial control region from ancient bones of 10 Holocene and four Pleistocene Siberian roe deer from Denisova cave as well as 37 modern specimen belonging to populations from Altai, Tian Shan (Kyrgyzstan), Yakutia, Novosibirsk region and the Russian Far East. Genealogical reconstructions indicated that most Holocene haplotypes were probably ancestral for modern roe deer populations of Western Siberia and Tian Shan. One of the Pleistocene haplotypes was possibly ancestral for modern Yakutian populations, and two extinct Pleistocene haplotypes were close to modern roe deer from Tian Shan and Yakutia. Most modern geographical populations (except for West Siberian Plains) are heterogeneous and there is some tentative evidence for structure. However, we did not find any distinct phylogenetic signal characterizing particular subspecies in either modern or ancient samples.

Conclusion/Significance

Analysis of mitochondrial DNA from both ancient and modern samples of Siberian roe deer shed new light on understanding the evolutionary history of roe deer. Our data indicate that during the last 50,000 years multiple replacements of populations of the Siberian roe deer took place in the Altai Mountains correlating with climatic changes. The Siberian roe deer represent a complex and heterogeneous species with high migration rates and without evident subspecies structure. Low genetic diversity of the West Siberian Plain population indicates a recent bottleneck or founder effect.     

Genetic Diversity of Moose (Alces alces L.) in Eurasia

Polymorphism of nucleotide sequence of D-loop fragment of the mitochondrial DNA was studied in 20 moose from several local populations on the territory of Eurasia. Three main haplotype variants of D-loop were detected by molecular phylogenetic method, which formed three clusters named European, Asian, and American. Intraspecies variation in the length of HVSI of D-loop of the mitochondrial DNA of moose was revealed. In the Far Eastern and Yakutian moose, haplotypes with a 75-bp deletion were found, which were most similar with haplotypes (also with the deletion), earlier observed in North American moose [1]. The highest diversity of the haplotypes of mitochondrial DNA is characteristic of Yakutia and the Far East (where three haplotype variants were found), which demonstrates the probable role of the region as the center of the species origin or as the region of ancient population mixture. The geographic region might be considered as a probable source of ancient moose migrations from Asia to America, basing on the data of distribution of mitochondrial haplotypes of D-loop and alleles of MhcAlal-DRB1. Divergence of nucleotide sequences of haplotypes with the 75-bp deletion (forming the American cluster on the phylogenetic tree) was the lowest (0.4%), which evidences respectively recent origin of the group of haplotypes. In Europe, only haplotypes of mitochondrial DNA referred to European variant were observed. Basing on analysis of variation of nucleotide sequences of D-loop, exon 4 of -Casein and exon 2 of MhcAlal-DRB1, we demonstrated that Eurasian moose studied belong to the unique species, which has probably passed through a bottle neck. The time of the origin of modern diversity of D-loop haplotypes of the species was estimated as 0.075–0.15 Myr ago.     

Mitochondrial DNA Sequence and Phylogenetic Evaluation of Geographically Disparate Sus scrofa Breeds

Next generation sequencing of mitochondrial DNA (mtDNA) facilitates studies into the metabolic characteristics of production animals and their relation to production traits. Sequence analysis of mtDNA from pure-bred swine with highly disparate production characteristics (Mangalica Blonde, Mangalica Swallow-bellied, Meishan, Turopolje, and Yorkshire) was initiated to evaluate the influence of mtDNA polymorphisms on mitochondrial function. Herein, we report the complete mtDNA sequences of five Sus scrofa breeds and evaluate their position within the phylogeny of domestic swine. Phenotypic traits of Yorkshire, Mangalica Blonde, and Swallow-belly swine are presented to demonstrate their metabolic characteristics. Our data support the division of European and Asian breeds noted previously and confirm European ancestry of Mangalica and Turopolje breeds. Furthermore, mtDNA differences between breeds suggest function-altering changes in proteins involved in oxidative phosphorylation such as ATP synthase 6 (MT-ATP6), cytochrome oxidase I (MT-CO1), cytochrome oxidase III (MT-CO3), and cytochrome b (MT-CYB), supporting the hypothesis that mtDNA polymorphisms contribute to differences in metabolic traits between swine breeds. Our sequence data form the basis for future research into the roles of mtDNA in determining production traits in domestic animals. Additionally, such studies should provide insight into how mtDNA haplotype influences the extreme adiposity observed in Mangalica breeds.     

Y-SNPs do not indicate hybridisation between European aurochs and domestic cattle

Background

Previous genetic studies of modern and ancient mitochondrial DNA have confirmed the Near Eastern origin of early European domestic cattle. However, these studies were not able to test whether hybridisation with male aurochs occurred post-domestication. To address this issue, Götherström and colleagues (2005) investigated the frequencies of two Y-chromosomal haplotypes in extant bulls. They found a significant influence of wild aurochs males on domestic populations thus challenging the common view on early domestication and Neolithic stock-rearing. To test their hypothesis, we applied these Y-markers on Neolithic bone specimens from various European archaeological sites.

Methods and Findings

Here, we have analysed the ancient DNA of 59 Neolithic skeletal samples. After initial molecular sexing, two segregating Y-SNPs were identified in 13 bulls. Strikingly, our results do not support the hypothesis that these markers distinguish European aurochs from domesticated cattle.

Conclusions

The model of a rapid introduction of domestic cattle into Central Europe without significant crossbreeding with local wild cattle remains unchallenged.     

Genetic diversity of moose (Alces alces L.) in Eurasia

Polymorphism of nucleotide sequence of D-loop fragment of the mitochondrial DNA was studied in 20 moose from several local populations on the territory of Eurasia. Three main haplotype variants of D-loop were detected by molecular phylogenetic method, which formed three clusters named European, Asian and American. Intraspecies variation in the length of HVSI of D-loop of the mitochondrial DNA of moose was revealed. In the Far Eastern and Yakutian moose, haplotypes with a 75-bp deletion were found, which were most similar with haplotypes (also with the deletion), earlier observed in North American moose [1]. The highest diversity of the haplotypes of mitochondrial DNA is characteristic of Yakutia and the Far East (where three haplotype variants were found), which demonstrates the probable role of the region as the center of the species or as the region of ancient population mixture. The geographic region might be considered as a probable source of ancient moose migrations from Asia to America, basing on the data of distribution of mitochondrial haplotypes of D-loop and alleles of MhcAlal-DRB1. Divergence of nucleotide sequences of haplotypes with the 75-bp deletion (forming the American cluster on the phylogenetic tree) was the lowest (0.4%), which evidences respectively recent origin of the group of haplotypes. In Europe, only haplotypes of mitochondrial DNA referred to European variant were observed. Basing on analysis of variation of nucleotide sequences of D-loop, exon 4 of kappa-Casein and exon 2 of MhcALal-DRB1, we demonstrated that Eurasian moose studied belong to the unique species, which has probably passed through a bottle neck. The time of the origin of modern diversity of D-loop haplotypes of the species was estimated as 0.075-0.15 Myr ago.     

Origin of Hungarian indigenous chicken breeds inferred from mitochondrial DNA D-loop sequences

In this study, we assessed the maternal origin of six Hungarian indigenous chicken breeds using mitochondrial DNA information. Sequences of Hungarian chickens were compared with the D-loop chicken sequences annotated in the GenBank and to nine previously described reference haplotypes representing the main haplogroups of chicken. The first 530 bases of the D-loop region were sequenced in 74 chickens of nine populations. Eleven haplotypes (HIC1-HIC11) were observed from 17 variable sites. Three sequences (HIC3, HIC8 and HIC9) of our chickens were found as unique to Hungary when searched against the NCBI GenBank database. Hungarian domestic chicken mtDNA sequences could be assigned into three clades and probably two maternal lineages. Results indicated that 86% of the Hungarian haplotypes are related to the reference sequence that likely originated from the Indian subcontinent, while the minor part of our sequences presumably derive from South East Asia, China and Japan.     

Are wild boars roaming Ireland once more?

Wild boars (Sus scrofa) have been increasingly sighted in the wild in Ireland during the last few years, likely due to illegal releases and/or escapees. The species has since been designated an invasive species in Ireland, which is seen as controversial by some because of uncertainties about the historic status of the species in Ireland. However, just as pertinent to the argument is the genetic purity of these individuals currently found in Ireland: are these pure wild boars? We carried out a genetic assessment of 15 wild boars shot in Ireland between 2009 and 2012 using 14 microsatellites and mitochondrial DNA (mtDNA). These were compared to European wild boar, domestic pig breeds and a hybrid population of ‘wild boar’ from England. Microsatellite analysis revealed that almost all the Irish individuals belonged to the ‘domestic pig’ genetic cluster, with only three individuals being classified as hybrids. All but two individuals carried Asian mtDNA haplotypes, indicating a domestic pig origin. It is clear from this study that the individuals currently found in Ireland are not pure wild boars and this result has to be factored into any management/eradication plans.     

Multiple maternal origins of native modern and ancient horse populations in China

To obtain more knowledge of the origin and genetic diversity of domestic horses in China, this study provides a comprehensive analysis of mitochondrial DNA (mtDNA) D-loop sequence diversity from nine horse breeds in China in conjunction with ancient DNA data and evidence from archaeological and historical records. A 247-bp mitochondrial D-loop sequence from 182 modern samples revealed a total of 70 haplotypes with a high level of genetic diversity. Seven major mtDNA haplogroups (A–G) and 16 clusters were identified for the 182 Chinese modern horses. In the present study, nine 247-bp mitochondrial D-loop sequences of ancient remains of Bronze Age horse from the Chifeng region of Inner Mongolia in China ( c. 4000–2000a bp ) were used to explore the origin and diversity of Chinese modern horses and the phylogenetic relationship between ancient and modern horses. The nine ancient horses carried seven haplotypes with rich genetic diversity, which were clustered together with modern individuals among haplogroups A, E and F. Modern domestic horse and ancient horse data support the multiple origins of domestic horses in China. This study supports the argument that multiple successful events of horse domestication, including separate introductions of wild mares into the domestic herds, may have occurred in antiquity, and that China cannot be excluded from these events. Indeed, the association of Far Eastern mtDNA types to haplogroup F was highly significant using Fisher's exact test of independence ( P  = 0.00002), lending support for Chinese domestication of this haplogroup. High diversity and all seven mtDNA haplogroups (A–G) with 16 clusters also suggest that further work is necessary to shed more light on horse domestication in China.     

Mitochondrial D-loop sequence variation among Italian horse breeds

The genetic variability of the mitochondrial D-loop DNA sequence in seven horse breeds bred in Italy (Giara, Haflinger, Italian trotter, Lipizzan, Maremmano, Thoroughbred and Sarcidano) was analysed. Five unrelated horses were chosen in each breed and twenty-two haplotypes were identified. The sequences obtained were aligned and compared with a reference sequence and with 27 mtDNA D-loop sequences selected in the GenBank database, representing Spanish, Portuguese, North African, wild horses and an Equus asinus sequence as the outgroup. Kimura two-parameter distances were calculated and a cluster analysis using the Neighbour-joining method was performed to obtain phylogenetic trees among breeds bred in Italy and among Italian and foreign breeds. The cluster analysis indicates that all the breeds but Giara are divided in the two trees, and no clear relationships were revealed between Italian populations and the other breeds. These results could be interpreted as showing the mixed origin of breeds bred in Italy and probably indicate the presence of many ancient maternal lineages with high diversity in mtDNA sequences.     

Genetic Footprints of Iberian Cattle in America 500 Years after the Arrival of Columbus

Background

American Creole cattle presumably descend from animals imported from the Iberian Peninsula during the period of colonization and settlement, through different migration routes, and may have also suffered the influence of cattle directly imported from Africa. The introduction of European cattle, which began in the 18th century, and later of Zebu from India, has threatened the survival of Creole populations, some of which have nearly disappeared or were admixed with exotic breeds. Assessment of the genetic status of Creole cattle is essential for the establishment of conservation programs of these historical resources.

Methodology/Principal Findings

We sampled 27 Creole populations, 39 Iberian, 9 European and 6 Zebu breeds. We used microsatellite markers to assess the origins of Creole cattle, and to investigate the influence of different breeds on their genetic make-up. The major ancestral contributions are from breeds of southern Spain and Portugal, in agreement with the historical ports of departure of ships sailing towards the Western Hemisphere. This Iberian contribution to Creoles may also include some African influence, given the influential role that African cattle have had in the development of Iberian breeds, but the possibility of a direct influence on Creoles of African cattle imported to America can not be discarded. In addition to the Iberian influence, the admixture with other European breeds was minor. The Creoles from tropical areas, especially those from the Caribbean, show clear signs of admixture with Zebu.

Conclusions/Significance

Nearly five centuries since cattle were first brought to the Americas, Creoles still show a strong and predominant signature of their Iberian ancestors. Creole breeds differ widely from each other, both in genetic structure and influences from other breeds. Efforts are needed to avoid their extinction or further genetic erosion, which would compromise centuries of selective adaptation to a wide range of environmental conditions.     

Assessing signatures of selection through variation in linkage disequilibrium between taurine and indicine cattle

Background

Signatures of selection are regions in the genome that have been preferentially increased in frequency and fixed in a population because of their functional importance in specific processes. These regions can be detected because of their lower genetic variability and specific regional linkage disequilibrium (LD) patterns.

Methods

By comparing the differences in regional LD variation between dairy and beef cattle types, and between indicine and taurine subspecies, we aim at finding signatures of selection for production and adaptation in cattle breeds. The VarLD method was applied to compare the LD variation in the autosomal genome between breeds, including Angus and Brown Swiss, representing taurine breeds, and Nelore and Gir, representing indicine breeds. Genomic regions containing the top 0.01 and 0.1 percentile of signals were characterized using the UMD3.1 Bos taurus genome assembly to identify genes in those regions and compared with previously reported selection signatures and regions with copy number variation.

Results

For all comparisons, the top 0.01 and 0.1 percentile included 26 and 165 signals and 17 and 125 genes, respectively, including TECRL, BT.23182 or FPPS, CAST, MYOM1, UVRAG and DNAJA1.

Conclusions

The VarLD method is a powerful tool to identify differences in linkage disequilibrium between cattle populations and putative signatures of selection with potential adaptive and productive importance.     

Genetic differentiation of European grayling (Thymallus thymallus) populations in Serbia,based on mitochondrial and nuclear DNA analyses

Background

The structure and diversity of grayling (Thymallus thymallus) populations have been well studied in most of its native habitat; however the southernmost populations of the Balkan Peninsula remain largely unexplored. The purpose of this study was to assess the genetic diversity of Serbian grayling populations, detect the impact of stocking and provide guidelines for conservation and management.

Methods

Eighty grayling individuals were collected from four rivers (Ibar, Lim, Drina and Rzav). The mitochondrial DNA control region (CR; 595 bp of the 3''end and 74 bp of flanking tRNA) and the ATP6 gene (630 bp fragment) were sequenced for 20 individuals (five from each locality). In addition, all individuals were genotyped with 12 microsatellite loci. The diversity and structure of the populations as well as the recent and ancient population declines were studied using specialized software.

Results

We detected three new haplotypes in the mtDNA CR and four haplotypes in the ATP6 gene of which three had not been described before. Previously, one CR haplotype and two ATP6 gene haplotypes had been identified as allochthonous, originating from Slovenia. Reconstruction of phylogenetic relations placed the remaining two CR haplotypes from the River Danube drainage of Serbia into a new clade, which is related to the previously described sister Slovenian clade. These two clades form a new Balkan clade. Microsatellite marker analysis showed that all four populations are genetically distinct from each other without any sign of intra-population structure, although stocking of the most diverse population (Drina River) was confirmed by mtDNA analysis. Recent and historical population declines of Serbian grayling do not differ from those of other European populations.

Conclusions

Our study shows that (1) the Ibar, Lim and Drina Rivers grayling populations are genetically distinct from populations outside of Serbia and thus should be managed as native populations in spite of some introgression in the Drina River population and (2) the Rzav River population is not appropriate for further stocking activities since it originates from stocked Slovenian grayling. However, the Rzav River population does not represent an immediate danger to other populations because it is physically isolated from these.