european huchen ( hucho hucho ) , also called danube salmon , is the largest salmonid species and one of the most endangered members of the salmonidae family ( holcik 1995 ) .
it is endemic to the danube drainage in central europe , occurring in cool montane and submontane reaches of large streams and swift rivers ( kottelat and freyhof 2007 ) .
european huchen , together with its sister species taimen ( hucho taimen ) are exclusively freshwater residents .
furthermore , this fish species can live for more than 20  years , and adults are large piscivorous predators of up to 1.65  m standard length ( sl ) and with a maximum mass of 60  kg ( kottelat and freyhof 2007 ) .
european huchen is a good indicator for overall stream connectivity and ecosystem health due to the species stringent water and substratum quality requirements as well as its facultative migratory spawning behavior .
european huchen is listed in annex ii and annex v of the european fauna flora habitat ( ffh ) directive 92/43/ewg and is considered as an endangered species with respect to its global distribution ( iucn 2013 ) , including poland ( prdboa 2002 ) . in recent years
the size of the european huchen populations has strongly declined through the effects of anthropogenic habitat alterations , such as : river regulation by constructing of dams and weirs , siltation of spawning ground , hydropower development as well as both industrial and agricultural pollution ( holcik 1990 ; witkowski et al .
2013a ) . in austria , european huchen inhabit 10  % of their previous distribution , considering self - sustaining populations ( schmutz et al .
distinctive characters of the species make european huchen a popular target species for anglers and a flagship species for conservation of running waters ( geist et al .
european huchen are bred in artificial conditions for a fishery supplementation of the local populations . according to some authors , due to low natural reproduction success
most of the wild populations of the european huchen depend on the stocking activities ( holcik 1995 ; witkowski et al .
european huchen occurs in dunajec , poprad , and san rivers , where it is protected by a fishery supplementation and enacted rigorous fishing restrictions . at present , there is only one fish farm , situated in lopuszna ( southern poland ) producing between 800 thousand to 1 million individuals of hatch annually , which is the only available source for european huchen stocking material in poland ( witkowski et al .
furthermore , available data suggests that hatchery stock of european huchen in lopuszna ( poland ) was established by small number of captured spawners . in 1955
only six ( three females and three males ) and in 1963 another few fish were caught from czarna orava river .
additionally , in 1985 , this stock was only once supplemented by 30 females originated from slovakia ( witkowski et al .
genetic data on the european huchen are still sparse and mainly limited to phylogenetic analysis of a few individuals in higher order systematic studies ( e.g. , phillips et al . 1995 ; crespi and fulton 2004 ) or research focused on the related taxa , such as : taimen ( froufe et al . 2005 ; guangxiang et al .
; liu et al . 2011 ) , lenok ( brachymystax lenok ) ( xia et al .
recently , a conservation genetic study on geographically limited european huchen populations from austria , slovenia , bosnia - herzogevenia , montenegro , germany , slovakia , and ukraine has been performed , which provided some fragmentary information about their genetic structure and variability ( geist et al . 2009 ; weiss et al . 2011 ) .
although information on the genetic diversity and differentiation is urgently needed for the conservation management of european huchen these data are still unavailable for the polish population of this fish species .
additionally , there are many studies that demonstrate the benefits of conservation management based on genetic data ( hansen 2002 ; gum et al .
thus , studies on the genetic structure of stock or population are important during conservation management for endangered fish species being protected by a rehabilitation measure , such as european huchen
. the major objectives of the present study were assessment of the current genetic diversity of european huchen broodstock in poland by means of microsatellite dna analysis and comparison of it with samples from german , slovak and ukrainian broodstocks .
the results of the research will provide baseline data for the improvement of existing conservation program management of the european huchen broodstock in poland .
fin clips from a total of 135 european huchen specimens were non - invasively sampled from four european huchen broodstocks and utilized for the genetic analysis .
fish tissues originated from fish farms localized in : poland ( restocking centre and trout hatchery lopuszna ) , germany ( fish farm lindbergmuehle , bavaria ) , slovakia ( fish farm pribovce , martin province ) and ukraine ( fish farm
baniliv , chemivtsi province ) during the years 2011 - 2013 ( fig .  1 ) .
small ( < 1  cm ) pelvic or pectoral fin clips were placed in eppendorf tubes and kept in 96  % ethanol at a temperature of 4  c until dna extraction .
dna was isolated from collected fin clips using standard chelex 100 procedure ( walsh et al . 1991).fig .
1sampling locations of european huchen : 1 . restocking centre and trout hatchery lopuszna , poland
a total of 29 microsatellite fragments , previously established for salmonids : bletet-2 , bletet-9 , bletri-2 ( froufe et al . 2004 ) , ssa-171 , ssa-197 , ssa-85 ( oreilly et al . 1996 ) , str-543 , str-85 ( presa and guymard 1996 ) , str-12 ( poteaux et al . 1999 ) , sfo-292 , sfo-262 ( perry et al . 2005 ) , str-15 , str-60 ( sonstebo et al . 2007 ) , str-73 ( estoup et al . 1993 ) , ogo-2 ( olsen et al . 1998 ) , sfo-18 ( geist et al . 2009 ) , hljz-003 , hljz-023 ( guangxiang et al .
2006 ) , omm-1016 , omm-1032 , omm-1077 , omm-1088 , omm-1097 , omm-1125 , omm-5000 , omm-5017 , hljz-031 , hljz-056 as well as hljz-069 ( you - yi et al .
reaction mixtures were prepared in a total volume of 25  l with a 0.8  l dna template ( 4.6    0.5  g / ml ) , 5.0  l of 10 pcr reaction buffer ( 100  mm tris - hcl ph  9.0 ,
500  mm kcl , triton x-100 ) , 0.5  l of each primer , 0.5  l ( 500  m ) of each deoxynucleotide triphosphate ( dntp ) , 0.8 - 1.0  l mgcl2 ( 25  mm/l ) and 0.6 unit gotaq flexi dna polymerase ( lucigen , middleton , wi , usa ) .
re - distilled water was used to bring the reaction mixture to the desired final volume .
amplification was performed with a mastercycler gradient thermocycler ( eppendorf , germany ) under the following conditions : an initial denaturation at 94  c for 3  min , followed by 34  cycles at 94  c for 30  s , annealing at 5361  c ( table  1 ) for 45  s , elongation at 72  c for 45  s and a final elongation step at 72  c for 7 to 10  min.table 1characterization of ten microsatellite loci applied in study of european huchen : locus designation , primer sequences , optimal annealing temperature ( t 
 a ) , mgcl2 concentration , number of observed alleles ( a 
 o ) , allele size range , fluorescence dye ( dye ) used for detection on the abi 3130 and source reference .
tetrasomic loci were underlinedlocusprimer sequence ( 5    3 ) 
 t 
 a ( c)mgcl2 ( mm ) 
 a 
 o 
 allele size range ( bp)dyereference 
 bletet-9 
 f : actggatagaaagacctgtgg530.815281422nedfroufe et al .
2004 
 r : ccacagctcagggcagggagt 
 hljz-056 
 f : ctctgtctcatctcgctt570.84224236petyou - yi et al .
1998 
 r : gtttcttcgactgtttcctctgtgttgag 
 omm-1032 
 f : gcgaggaagagaaagtagtag570.87203227nedyou - yi et al .
2009 
 r : cccatcttctctctgattatg 
 omm-1077 
 f : ggctgaccagagaaagactagttc601.026325467petyou - yi et al .
2009 
 r : tgttacggtgtctgacatgc 
 omm-1088 
 f : ctacaggccaacactacaatc610.86112159petyou - yi et al .
1996 
 r : tggcagggatttgacataac 
  characterization of ten microsatellite loci applied in study of european huchen : locus designation , primer sequences , optimal annealing temperature ( t 
 a ) , mgcl2 concentration , number of observed alleles ( a 
 o ) , allele size range , fluorescence dye ( dye ) used for detection on the abi 3130 and source reference .
tetrasomic loci were underlined genotyping of microsatellite dna fragments was conducted using an applied biosystems 3130 genetic analyzer . in order to enable genotyping of pcr products forward primers were labeled with different fluorescent reporter dyes ( pet - red , vic - green , 6-fam - blue and ned - yellow ) ( table  1 ) . the genescan 600 liz size standard was utilized as a reference for determining the length of examined dna fragments .
individual microsatellite loci amplified using primers with different attached fluorescent dyes were arranged into sets and analyzed in multiplex mode . in order to visualize the results , software provided by manufacturer genemaper v4.1 software and data collection software v3.0 ( applied biosystems , california , usa ) were used .
the genetic profiles containing the list of alleles detected within the studied loci were prepared for each fish .
the micro - checker software ( version 2.2.3 ) was used to check microsatellite null alleles , scoring errors due to stuttering and large allele drop - out in samples ( van oosterhout et al .
the observed number of alleles per locus , allele frequency , number of private alleles , allelic range and allelic richness  ( a 
 r ) were computed by genepop software ( version 4.2.1 ) ( rousset 2008 ) .
the observed ( ho ) and expected heterozygosity ( he ) , the exact hardy - weinberg ( h - w ) equilibrium test , as well as linkage disequilibrium ( ld ) were calculated using arlequin software ( version 3.5 ) ( excoffier and lischer 2010 ) .
the genetic differences between the sampled broodstocks and fixation index ( fis ) were calculated using fstat software ( version 2.9.3 ) ( goudet 2001 ) .
the polymorphism information content ( pic value ) was also calculated employing powermarker software ( version 3.25 ) ( liu and muse 2005 ) . genetic divergence between studied broodstocks of european huchen was analyzed using two different parameters : genetic differentiation index ( fst ) and nei s genetic distance .
shannon s index ( i ) for each loci within tested broodstocks was calculated using popgene software ( version 1.3.2 ) ( yeh and boylet 1997 ) .
the upgma dendrogram based on nei s genetic distance was constructed by mega6 ( version 6.0.5 ) ( tamura et al . 2013 ) .
an analysis of molecular variance ( amova ) was done with the arlequin package 3.5 for measuring variance within and between broodstocks .
the likely occurrence of bottleneck or the founder effect and their influence on within - broodstock genetic variability was based on the garza - williamson index  ( m ) , which was computed applying arlequin software 3.5 .
a test for bottleneck assessment was also conducted using the bottleneck software ( version 1.9 ) ( piry et al .
1999 ) , which tests for departure from mutation drift equilibrium based on heterozygosity excess or deficiency .
recent broodstock bottlenecks assuming a stepwise mutation model ( smm ) and infinite allele model ( iam ) for four broodstocks of european huchen were tested .
this method is based on the assumption that in non - bottlenecked broodstock ( close to mutation drift equilibrium ) the value of expected heterozygosity ( he ) is equal to heq ( heterozygosity expected in a mutation - drift equilibrium ) .
the excess of he over heq is the evidence of severe reduction in broodstock effective size that may occur because of a bottleneck event . in order to accommodate the obtained genotypic data to the requirements of employed software , every tetrasomic locus
was examined as two disomic loci and as result the mean values was considered for estimation of genetic parameters .
fin clips from a total of 135 european huchen specimens were non - invasively sampled from four european huchen broodstocks and utilized for the genetic analysis .
fish tissues originated from fish farms localized in : poland ( restocking centre and trout hatchery lopuszna ) , germany ( fish farm lindbergmuehle , bavaria ) , slovakia ( fish farm pribovce , martin province ) and ukraine ( fish farm
baniliv , chemivtsi province ) during the years 2011 - 2013 ( fig .  1 ) .
small ( < 1  cm ) pelvic or pectoral fin clips were placed in eppendorf tubes and kept in 96  % ethanol at a temperature of 4  c until dna extraction .
dna was isolated from collected fin clips using standard chelex 100 procedure ( walsh et al . 1991).fig .
1sampling locations of european huchen : 1 . restocking centre and trout hatchery lopuszna , poland
a total of 29 microsatellite fragments , previously established for salmonids : bletet-2 , bletet-9 , bletri-2 ( froufe et al . 2004 ) , ssa-171 , ssa-197 , ssa-85 ( oreilly et al . 1996 ) , str-543 , str-85 ( presa and guymard 1996 ) , str-12 ( poteaux et al . 1999 ) , sfo-292 , sfo-262 ( perry et al . 2005 ) , str-15 , str-60 ( sonstebo et al . 2007 ) , str-73 ( estoup et al . 1993 ) , ogo-2 ( olsen et al . 1998 ) , sfo-18 ( geist et al . 2009 ) , hljz-003 , hljz-023 ( guangxiang et al .
2006 ) , omm-1016 , omm-1032 , omm-1077 , omm-1088 , omm-1097 , omm-1125 , omm-5000 , omm-5017 , hljz-031 , hljz-056 as well as hljz-069 ( you - yi et al . 2009 ) were tested for cross - amplification in four broodstocks of european huchen .
reaction mixtures were prepared in a total volume of 25  l with a 0.8  l dna template ( 4.6    0.5  g / ml ) , 5.0  l of 10 pcr reaction buffer ( 100  mm tris - hcl ph  9.0 ,
500  mm kcl , triton x-100 ) , 0.5  l of each primer , 0.5  l ( 500  m ) of each deoxynucleotide triphosphate ( dntp ) , 0.8 - 1.0  l mgcl2 ( 25  mm/l ) and 0.6 unit gotaq flexi dna polymerase ( lucigen , middleton , wi , usa ) .
re - distilled water was used to bring the reaction mixture to the desired final volume .
amplification was performed with a mastercycler gradient thermocycler ( eppendorf , germany ) under the following conditions : an initial denaturation at 94  c for 3  min , followed by 34  cycles at 94  c for 30  s , annealing at 5361  c ( table  1 ) for 45  s , elongation at 72  c for 45  s and a final elongation step at 72  c for 7 to 10  min.table 1characterization of ten microsatellite loci applied in study of european huchen : locus designation , primer sequences , optimal annealing temperature ( t 
 a ) , mgcl2 concentration , number of observed alleles ( a 
 o ) , allele size range , fluorescence dye ( dye ) used for detection on the abi 3130 and source reference .
tetrasomic loci were underlinedlocusprimer sequence ( 5    3 ) 
 t 
 a ( c)mgcl2 ( mm ) 
 a 
 o 
 allele size range ( bp)dyereference 
 bletet-9 
 f : actggatagaaagacctgtgg530.815281422nedfroufe et al .
2004 
 r : ccacagctcagggcagggagt 
 hljz-056 
 f : ctctgtctcatctcgctt570.84224236petyou - yi et al .
1998 
 r : gtttcttcgactgtttcctctgtgttgag 
 omm-1032 
 f : gcgaggaagagaaagtagtag570.87203227nedyou - yi et al .
2009 
 r : cccatcttctctctgattatg 
 omm-1077 
 f : ggctgaccagagaaagactagttc601.026325467petyou - yi et al .
2009 
 r : tgttacggtgtctgacatgc 
 omm-1088 
 f : ctacaggccaacactacaatc610.86112159petyou - yi et al .
1996 
 r : tggcagggatttgacataac 
  characterization of ten microsatellite loci applied in study of european huchen : locus designation , primer sequences , optimal annealing temperature ( t 
 a ) , mgcl2 concentration , number of observed alleles ( a 
 o ) , allele size range , fluorescence dye ( dye ) used for detection on the abi 3130 and source reference .
genotyping of microsatellite dna fragments was conducted using an applied biosystems 3130 genetic analyzer . in order to enable genotyping of pcr products forward
primers were labeled with different fluorescent reporter dyes ( pet - red , vic - green , 6-fam - blue and ned - yellow ) ( table  1 ) .
the genescan 600 liz size standard was utilized as a reference for determining the length of examined dna fragments .
individual microsatellite loci amplified using primers with different attached fluorescent dyes were arranged into sets and analyzed in multiplex mode . in order to visualize the results , software provided by manufacturer genemaper v4.1 software and data collection software v3.0 ( applied biosystems , california , usa )
the genetic profiles containing the list of alleles detected within the studied loci were prepared for each fish .
the micro - checker software ( version 2.2.3 ) was used to check microsatellite null alleles , scoring errors due to stuttering and large allele drop - out in samples ( van oosterhout et al .
the observed number of alleles per locus , allele frequency , number of private alleles , allelic range and allelic richness  ( a 
 r ) were computed by genepop software ( version 4.2.1 ) ( rousset 2008 ) .
the observed ( ho ) and expected heterozygosity ( he ) , the exact hardy - weinberg ( h - w ) equilibrium test , as well as linkage disequilibrium ( ld ) were calculated using arlequin software ( version 3.5 ) ( excoffier and lischer 2010 ) .
the genetic differences between the sampled broodstocks and fixation index ( fis ) were calculated using fstat software ( version 2.9.3 ) ( goudet 2001 ) .
the polymorphism information content ( pic value ) was also calculated employing powermarker software ( version 3.25 ) ( liu and muse 2005 ) . genetic divergence between studied broodstocks of european huchen
was analyzed using two different parameters : genetic differentiation index ( fst ) and nei s genetic distance .
shannon s index ( i ) for each loci within tested broodstocks was calculated using popgene software ( version 1.3.2 ) ( yeh and boylet 1997 ) .
the upgma dendrogram based on nei s genetic distance was constructed by mega6 ( version 6.0.5 ) ( tamura et al . 2013 ) .
an analysis of molecular variance ( amova ) was done with the arlequin package 3.5 for measuring variance within and between broodstocks .
the likely occurrence of bottleneck or the founder effect and their influence on within - broodstock genetic variability was based on the garza - williamson index  ( m ) , which was computed applying arlequin software 3.5 . a test for bottleneck assessment
was also conducted using the bottleneck software ( version 1.9 ) ( piry et al .
1999 ) , which tests for departure from mutation drift equilibrium based on heterozygosity excess or deficiency .
recent broodstock bottlenecks assuming a stepwise mutation model ( smm ) and infinite allele model ( iam ) for four broodstocks of european huchen were tested .
this method is based on the assumption that in non - bottlenecked broodstock ( close to mutation drift equilibrium ) the value of expected heterozygosity ( he ) is equal to heq ( heterozygosity expected in a mutation - drift equilibrium ) .
the excess of he over heq is the evidence of severe reduction in broodstock effective size that may occur because of a bottleneck event . in order to accommodate the obtained genotypic data to the requirements of employed software , every tetrasomic locus
was examined as two disomic loci and as result the mean values was considered for estimation of genetic parameters .
in the present study the 29 microsatellite dna fragments were applied to evaluate the genetic diversity in four broodstocks of european huchen . among these 29 microsatellites
, 11 ( str-12 , str-15 , str-60 , str-73 , str-543 , ssa-85 , ssa-171 , bletet-2 , sfo-292 , omm-1016 , omm-5017 ) were monomorphic in all broodstocks , one ( str-85 ) did not produce any pcr product and seven loci ( hljz-003 , hljz-023 , hljz-031 , hljz-069 , omm-1097 , omm-1125 , omm-5000 ) produced non specified bands ( stutter bands ) .
ten microsatellite loci ( bletet-9 , bletri-2 , hljz-056 , ogo-2 , omm-1032 , omm-1077 , omm-1088 , sfo-18 , sfo-262 , ssa-197 ) were selected for further analyses because they were characterized by good quality products and were polymorphic ( table  1 ) . among these ,
loci : bletet-9 , bletri-2 , hljz-056 , ogo-2 , omm-1032 , omm-1088 , sfo-18 , ssa-197 were considered as disomic and omm-1077 , sfo-262 were tetrasomic .
the size of the alleles in an individual locus varied between 112 base pairs ( bp ) and 467  bp .
the number of amplified alleles per locus ranged from one ( sfo-262 ) to 26 ( omm-1077 ) with an average of 8.1 alleles per locus ( table  1 ) .
the examined broodstocks of european huchen differed in the number of alleles detected in locus as well as in the overall number of alleles identified across all investigated loci .
moreover , the allelic frequency distribution for a number of loci was quite different among the four studied broodstocks .
the null alleles were detected in german and slovak broodstocks , which appeared in omm-1077 loci at 0.086 and 0.063 frequencies , respectively ( table  2).table 2allele frequencies at microsatellite loci in studied european huchen broodstocks .
null : null allelelocusallelesallele frequencies by locuspolandgermanyslovakiaukraine 
 sfo-262 
 1861.0000.5160.6931.0002200.0000.4450.2160.0002280.0000.0390.0910.000 
 ssa-197 
 1320.0000.0000.0910.0001360.0000.1090.0000.0001400.0000.2340.0110.0001440.0000.0000.0910.0001480.9330.5000.2950.9481560.0170.0940.0910.0001600.0000.0000.1360.0001800.0500.0630.2840.052 
 sfo-18 
 1630.0000.0780.0000.0001680.4330.6090.2730.4661940.3330.0780.6480.2412070.2330.2340.0800.293 
 ogo-2 
 2480.0000.0310.0000.0002560.0000.2030.1480.0002580.4000.4840.2610.2412660.6000.2810.5910.7242700.0000.0000.0000.034 
 bletri-2 
 1180.0000.2340.0800.0341240.0500.0000.0000.0171270.9500.7660.9200.948 
 bletet-9 
 2810.0000.0940.0000.0002920.0000.1720.0000.0003040.1830.0000.0000.5173080.2000.1880.0680.0523120.0330.0310.0000.0003370.0000.0160.1700.0003410.0000.0000.0110.0003530.4170.3750.0910.2073570.0830.0470.0000.1903650.0330.0780.1360.0003900.0000.0000.1360.0003940.0330.0000.0000.0004100.0170.0000.0000.0344180.0000.0000.3750.0004220.0000.0000.0110.000 
 hljz-056 
 2240.0000.0470.0000.0002300.1170.0630.1140.1382330.1830.3440.0000.1382360.7000.5470.8860.724 
 omm-1032 
 2030.2830.0000.2270.2762050.0000.2340.0800.0002090.0670.0780.0000.2592110.4830.4060.6930.4482150.0000.0630.0000.0002250.1670.0000.0000.0172270.0000.2190.0000.000 
 omm-1088 
 1120.1330.2970.4090.0691160.4170.1090.1360.5001200.2670.4380.2610.1381510.0000.0000.0000.0171550.1830.0310.1820.2761590.0000.1250.0110.000 
 omm-1077 

null : null allele one microsatellite locus , sfo-262 , was monomorphic in polish and ukrainian broodstocks ( table  2 ) .
eight microsatellite loci , bletet-9 , hljz-056 , ogo-2 , omm-1032 , omm-1077 , omm-1088 , sfo-18 , and ssa-197 , were highly polymorphic in the groups of fish studied ( table  2 ) .
the high degree of polymorphism of studied loci implies that each locus is informative and could be used in population studies .
the genetic diversity parameters ( ho , he , a 
 o , a 
 e , a 
 r , i , and pic ) of the four broodstocks of european huchen are shown in tables  3 and 4 .
the mean values of the polymorphism information content ( pic ) in fish from poland , germany , slovakia , and ukraine were 0.430 , 0.590 , 0.496 , and 0.405 , respectively .
the mean allelic richness varied from 3.944 to 5.200 in studied broodstocks of european huchen ( table  3 ) .
the german broodstock was characterized by the highest number of the private alleles ( 15 ) in the studied microsatellite loci.table 3genetic diversity parameters of four european huchen broodstocks analyzed .
a 
 r : allelic richness , a 
 o : observed alleles , a 
 e : expected alleles , i : shannon s index , pic : polymorphism information content , mono : monomorphic locusbroodstocklocus 
 a 
 r 

 a 
 o 

 a 
 e 

 i 

 pic 
 private allelespoland 
 sfo-262 
 1.0001mono0.00000.000 
 ssa-197 
 2.96731.14430.28240.121 
 sfo-18 
 3.00032.83021.06810.572 

3.00031.86140.81130.416 
 omm-1032 
 4.00042.88921.18790.596 
 omm-1088 
 4.00043.37711.29690.653 
 omm-1077 
 11.349145.34451.90500.777346 , 359 , 418mean4.1284.42.70600.91180.430germany 
 sfo-262 
 2.95332.14580.82030.426 
 ssa-197 
 5.00053.03411.32390.627136 
 sfo-18 
 4.00042.28061.04020.510163 
 ogo-2 
 3.99342.80931.14000.580248 
 bletri-2 
 2.00021.55980.54450.294 
 bletet-9 
 7.89884.47161.72220.748281 , 292 
 hljz-056 
 3.99942.36221.01390.501224 
 omm-1032 
 5.00053.59931.41090.678215 , 227 
 omm-1088 
 4.99353.24561.33250.642 
 omm-1077 
 12.167177.61712.19390.845340 , 349 , 353 , 355 , 402 , 416 , 467mean5.2005.63.31251.25420.590slovakia 
 sfo-262 
 2.99831.86020.79370.407 
 ssa-197 
 6.65974.72771.69430.756132 , 144 , 160 
 sfo-18 
 3.00031.99900.83700.431 
 ogo-2 
 3.00032.27630.94410.495 
 bletri-2 
 2.00021.17160.27770.136 
 bletet-9 
 7.31784.54461.71560.753341 , 390 , 418 , 422 
 hljz-056 
 2.00021.25230.35410.181 
 omm-1032 
 3.00031.85710.79210.405 
 omm-1088 
 4.65953.47891.34890.664 
 omm-1077 
 7.589115.70251.71310.733378 , 447mean4.2224.72.88701.04700.496ukraine 
 sfo-262 
 1.0001mono0.00000.000 
 ssa-197 
 2.00021.10880.20360.093 
 sfo-18 

3.00032.77101.05870.567 
 ogo-2 
 3.00031.71280.69290.354270 
 bletri-2 
 3.00031.11020.23650.097 
 bletet-9 
 5.00052.85571.25160.602 
 hljz-056 
 3.00031.77800.78020.397 
 omm-1032 
 4.00042.90501.13470.588 
 omm-1088 
 5.00052.85571.22950.596151 
 omm-1077 
 7.500105.09371.72510.755mean3.9444.222.46560.8310.405 
 table 4comparison of observed ( ho ) and expected ( he ) heterozygosity , expected heterozygosity ( heq ) in an infinite allele model ( iam ) and stepwise mutation model ( smm ) , garza - williamson index ( m ) and fixation index ( fis ) in studied european huchen boodstocks : p - level of significance , mono : monomorphic locus .
deviations statistically significant at p  <  0.05broodstock / locus 
 ho 

 he 

 p 

 iam 

 smm 

 fis 

 m 

 heq 

 p 

 heq 

 p 
 poland 
 sfo-262 
 mono0.0001.000 
 ssa-197 
 0.1330.1281.0000.3440.1640.4760.0180.0400.091 
 sfo-18 
 0.8000.6580.3240.3640.0240.4730.0410.2210.075 
 ogo-2 
 0.4670.4881.0000.2050.0800.2500.1270.0450.222 
 bletri-2 
 0.1000.0971.0000.2020.4270.2560.2870.0360.500 
 bletet-9 
 0.6670.7550.0450.7090.4040.8130.0820.1190.075 
 hljz-056 
 0.4670.4710.8830.3580.3530.4760.3810.0090.429 
 omm-1032 
 0.7670.6650.1970.4680.0830.6080.3290.1560.174 
 omm-1088 
 0.8000.7160.6620.4620.0210.6040.0770.1200.091 
 omm-1077 
 1.0000.8160.0740.8120.4330.8790.0200.2410.135  mean0.5200.4780.5760.4360.2210.5370.1510.0710.279germany 
 sfo-262 
 0.9690.5420.0000.3560.1670.4710.3440.8140.070 
 ssa-197 
 0.6870.6810.7570.5390.1800.6820.4090.0100.111 
 sfo-18 
 0.5940.5701.0000.4590.3100.5980.3380.0420.089 
 ogo-2 
 0.7500.6540.3510.4630.1080.6020.3800.1490.210 
 bletri-2 
 0.4060.3650.6530.2160.2760.2460.3290.1160.200 
 bletet-9 
 0.5940.7890.0070.7010.2170.8100.2480.2500.094 
 hljz-056 
 0.6250.5860.7450.4580.2640.5950.3720.0680.308 
 omm-1032 
 0.7500.7340.8840.5500.0650.6840.3040.0230.217 
 omm-1088 
 0.5620.7030.0340.5390.1430.6770.4700.2020.104 
 omm-1077 
 0.9220.8760.0020.8190.1480.8860.2210.0550.058  mean0.6860.6500.4430.5100.1880.6250.3420.0820.146slovakia 
 sfo-262 
 0.6140.4680.0360.3350.3100.4640.4070.3170.070 
 ssa-197 
 0.9540.7970.0000.6360.0500.7710.3630.2000.143 
 sfo-18 
 0.5230.5050.3200.3300.2110.4580.4420.0350.075 
 ogo-2 
 0.6820.5670.2640.3300.0980.4600.2290.2050.273 
 bletri-2 
 0.1590.1481.0000.1880.4910.2330.4250.0750.200 
 bletet-9 
 0.6590.7890.0030.6740.1380.8030.3090.1660.070 
 hljz-056 
 0.1820.2040.4370.1930.3990.2420.4970.1090.286 
 omm-1032 
 0.2500.4670.0000.3370.3140.4640.4120.4670.333 
 omm-1088 
 0.8180.7210.0860.5240.0640.6680.2920.1370.104 
 omm-1077 
 0.9660.7850.0000.6450.0880.7640.2560.2460.033  mean0.5810.5450.2150.4190.2160.5330.3630.0470.159ukraine 
 sfo-262 
 mono0.0001.000 
 ssa-197 
 0.1030.1001.0000.2100.4240.2460.3140.0370.061 
 sfo-18 
 0.6900.6500.3830.3610.0240.4790.0600.0620.075 
 ogo-2 
 0.2410.4230.0300.3540.4190.4730.3100.4340.231 
 bletri-2 
 0.1030.1011.0000.3540.1240.4730.0130.0240.300 
 bletet-9 
 0.5520.6610.0340.5570.2980.6850.2950.1680.047 
 hljz-056 
 0.4830.4450.5330.3630.4040.4790.3230.0860.429 
 omm-1032 
 0.6550.6670.5210.4650.1010.6050.2900.0180.174 
 omm-1088 
 0.5520.6610.0860.5510.2740.6860.3110.1680.114 
 omm-1077 
 0.9650.8170.0060.7110.0910.8100.3940.1000.116 
a 
 r : allelic richness , a 
 o : observed alleles , a 
 e : expected alleles , i : shannon s index , pic : polymorphism information content , mono : monomorphic locus comparison of observed ( ho ) and expected ( he ) heterozygosity , expected heterozygosity ( heq ) in an infinite allele model ( iam ) and stepwise mutation model ( smm ) , garza - williamson index ( m ) and fixation index ( fis ) in studied european huchen boodstocks : p - level of significance , mono : monomorphic locus .
deviations statistically significant at p  <  0.05 the mean observed heterozygosity ( ho ) in the investigated broodstocks ranged from 0.434 ( ukraine ) to 0.686 ( germany ) and were close to the mean values expected under h - w equilibrium ( he ) ( table  4 ) .
slovak broodstock had the highest level of deviation ( five loci ) and three of these loci exhibited heterozygote excess ( table  4 ) .
table  4 also shows expected heterozygosity in two models of a mutation - drift equilibrium ( heq ) . in all broodstocks , under the infinite allele model ( iam ) and stepwise mutation ( smm ) models ,
heterozygosity excess were detected in most of analyzed loci ; however , observed he > heq differences were significant in only polish ( two out of ten loci were analyzed ) and ukrainian ( one out of ten loci were analyzed ) broodstocks . all investigated loci differed in terms of the garza - williamson index ( m )
m values were the lowest for the most polymorphic loci ( such as bletet-9 and omm-1077 ) and the highest for the least polymorphic loci ( bletri-2 , hljz-056 ) .
the mean observed m values in the examined broodstocks varied from 0.146 to 0.279 . the mean fixation index ( fis ) in fish from poland , germany , slovakia , and ukraine were 0.071 , 0.082 , 0.047 , and 0.053 , respectively ( table  4 ) .
the genetic differentiation ( fst ) of the four broodstocks of european huchen was 0.1139 .
analysis of the genetic structure of the studied broodstocks with amova method revealed that 13.20  % of the genetic diversity was distributed among broodstocks and 86.80  % occurred among individuals within the broodstocks .
the nei s genetic distances among the four broodstocks and the fst matrix are shown in table  5 .
the highest genetic divergence , calculated as fst values , occurred between ukraine and slovakia ( 0.1731 ) . the lowest genetic divergence and genetic distance  were observed between poland and ukraine ( 0.0234 and 0.0333 , respectively ) , while the highest genetic distances were observed among german , slovak , and ukrainian broodstocks ( fig . 
2).table 5nei s genetic distance ( below diagonal ) and fst ( above diagonal ) of four european huchen broodstocks analyzedbroodstockpolandgermanyslovakiaukrainepoland0.13610.13610.0234germany0.22260.11920.1711slovakia0.19280.24270.1731ukraine0.03330.28570.2483 
 fig .
2upgma dendrogram based on nei s genetic distance illustrating relationships among four european huchen broodstocks under current study 
  nei s genetic distance ( below diagonal ) and fst ( above diagonal ) of four european huchen broodstocks analyzed upgma dendrogram based on nei s genetic distance illustrating relationships among four european huchen broodstocks under current study
such molecular based screening is essential for understanding population genetic differentiation among reared stocks within species , inferring parentage in mixed - family assemblages , maintaining genetic variability in populations , estimating the effective size of populations as well as inferring the effects of selection within raised stocks ( vrijenhoek 1998 ; frankham et al .
microsatellite dna analysis technique successfully applied for the first time in the polish broodstock of european huchen provided new information about genetic structure of this valuable salmonid fish species .
the comparison of the overall number of private alleles , allelic ranges and its frequency clearly demonstrates genetic differences between studied broodstocks of european huchen .
this confirms the hypothesis that the broodstocks studied do not share a common gene pool at the population genetic level
. genetic diversity , which can be evaluated as the allele number by locus and the mean number across loci , differed among the studied groups of fish .
examined broodstocks of european huchen exhibited moderate genetic diversity ( the mean he and pic per broodstock were between 0.40.5 ) with the exception of german broodstock which presented high genetic diversity ( the mean he and pic were close to 0.6 ) .
similar moderate indices of genetic diversity to those observed in the present study were reported for some populations of taimen from heilongjiang river basin ( china ) ( guangxiang et al .
2006 ; you - yi et al . 2009 ; liu et al . 2011 ) and for the huchen populations from europe ( geist et al . 2009 ; weiss et al .
2011 ) , evidencing that our results are comparable to those described in other studies ( liang et al .
all genetic diversity parameters ( ho , he , a 
 o , a 
 e , a 
 r , i , and pic ) consistently ranked the broodstocks in order of decreasing diversity ( germany > slovakia
> poland > ukraine ) . however , there were no significant differences between these parameters among the four studied broodstocks ( p  >  0.05 ) .
however , from 8 to 42  % of the studied loci deviated from h - w equilibrium , suggesting that both the genotype and gene frequencies fluctuated continuously as a response to the different stocking conditions .
the average fis value was positive only in ukrainian broodstock ( 0.053 ) indicating a small overall deficiency of heterozygotes , while the rest of the analyzed broodstocks of european huchen were characterized by negative fis indicators exhibiting small overall excess of heterozygotes against hardy - weinberg expectations .
in contrast , in domesticated stocks of fish the observed heterozygosity ( ho ) sometimes exceeds the expected heterozygosity ( he ) .
for example , this tendency was observed in chinook salmon ( oncorhynchus tshawytscha ) stocks ( kim et al .
2004 ) and paddlefish ( polyodon spathula ) ( kaczmarczyk et al . 2012 ) .
this excess may be a consequence of the use of a non - random subset of the broodstock in the hatchery conditions ( luikart and cornuet 1998 ; 1999 ) .
the bottleneck test applied in the present study showed evidence of a bottleneck under the infinite allele ( iam ) and stepwise mutation models ( smm ) in polish and ukrainian broodstocks .
the heterozygosity excess at observed loci in mentioned broodstocks might be indicative of a small founding population size or bottleneck events as was hypothesized in the case of the hatchery stocks of brook trout ( salvelinus fontinalis ) ( fopp - bayat et al .
2010 ) and paddlefish ( kaczmarczyk et al . 2012 ) . in the present study
the average values of the garza - williamson index for every examined broodstock were surprisingly low ( 0.1460.279 ) suggesting that all investigated broodstocks suffered from bottleneck or founder effects in the past ( garza and williamson 2001 ; tzika et al .
the historical information together with genetic data on polish broodstock of european huchen confirm the occurrence of the bottleneck effect in the past .
our results showed that some genetic diversity among the broodstocks occurred ( 13.20  % ) ; however , the majority of diversity were observed between individuals within a broodstocks ( 86.80  % ) .
similar genetic diversity ( 13.10  % ) was observed among two hatchery stocks of barfin flounder ( verasper moseri ) from china ( hongyu et al .
lower genetic diversity were reported between four taimen populations from heilongjiang river basin ( 6.12  % ) ( you - yi et al .
these presented data suggest that there is substantial genetic differentiation between tested broodstocks of european huchen . a low level ( 0.023 ) of genetic differentiation
contrastingly , the genetic distances observed between german , slovak , and ukrainian broodstocks were moderate because the rate of genetic distance of these groups of fish was 0.1731 ( balloux and lugon - moulin 2002 ) .
numerous examples evidencing that mixed - source reintroductions by genetically distant populations may result in outbreeding depression ( gharrett et al . 1999 ; alacs et al . 2007 ; huff et al .
thus , any enrichment to the genetic pool of a conserved population should be done with genetically similar material . according to the present results ,
if a supplement of the genetic pool of polish broodstock of european huchen will be necessary in the future , the ukrainian broodstock seems to be the most suitable material to maintain the genetic diversity of this stock .
a continuous genetic monitoring of hatchery stocks is an important tool for development of sustainable conservation management programs .
management based on supportive and captive breeding carries the risk of adverse modifications to a gene pool , resulting in loss of genetic diversity , decrease of heterozygosity , or inbreeding depression .
the consequence of inbreeding depression in hatchery stocks can be the negative effects on : hatching rate , fry survival , growth rate , and spawning performance , as well as inefficiency of feed conversion , occurrence of deformations and presence of short - lived albino fish ( guo - sheng et al .
; wang et al . 2002 ; ala - honkola et al . 2009 ) .
supportive breeding systems are complex and the long - term genetic consequences are unclear and difficult to predict ( duchesne and bernatchez 2002 ) . for example
, the common practice of mixed - milt fertilization , systematic selection of specified phenotypes and feeding of salmonid fish with trout pellets under hatchery conditions may lead to selection of genotypes that have lower viability in the wild ( flagg et al .
the efficacy of artificial stocking in order to maintain natural populations is one of the most controversial topics in fishery management .
it is considered that the post - release survival rate of reared fish is essentially lower than that of wild fish ( hoekstra et al . 2007 ; araki et al . 2008 ) .
the long - term assessment of results of historic stocking of salmonids in europe , including european huchen are very poorly known .
thus , the effects of hatchery fish on wild populations remain an open question and a topic of major concern . in the case of european huchen ,
stocking activities have not been well documented and no official information is therefore available regarding the purpose of stocking and the origin or the numbers of stocked fish .
the potential for evolution and adaptation to new environments is limited by the genetic diversity of the population .
natural selection favors individuals that are better adapted to live in natural environment than in captivity , thus research programs aimed at retaining the evolutionary genetic adaptive potential of european huchen stocks are very important .
therefore the optimum degree of genetic differentiation and a high genetic diversity should be permanently controlled . moreover
, european huchen habitat continues to be threatened by further anthropogenic factors influence and particularly unfavorable conditions may adversely impact on the populations characterized by low genetic variation level .
the lack of information on the genetic structure of european huchen broodstock in poland essentially limits the sustainable conservation of this species . therefore , baseline genetic data are crucial to guide future population specific conservation programs and research efforts on european huchen in poland .
 four broodstocks of european huchen ( hucho hucho ) from : poland , germany , slovakia , and ukraine were investigated using ten microsatellite dna loci . 
 microsatellite dna analysis was successfully applied for the first time in the polish broodstock of this fish species . 
 the genetic variation and genetic distance between these broodstocks were evaluated . 
 in addition , we examined the potential effects of a genetic bottleneck on the genetic variation of the broodstocks . 
 the european huchen broodstocks exhibited moderate genetic diversity ( pic  =  0.4050.496 and i  =  0.8311.047 ) with the exception of german broodstock which presented higher genetic diversity ( pic  =  0.590 and i  =  1.254 ) . observed ( ho ) and expected ( he ) heterozygosity across the investigated loci in all broodstocks ranged from 0.434 to 0.686 and from 0.452 to 0.650 , respectively . 
 overall , the studied broodstocks were in hardy - weinberg equilibrium ( hwe ) ; however , from 8 to 42  % of the loci deviated from hwe in each stock . the garza - williamson index ( m  =  0.1460.279 ) and values of the heterozygosity excess revealed a reduction of genetic variation in all studied broodstocks because of the founder or bottleneck effect 
 . the analysis of genetic differentiation ( fst ) and nei s genetic distance between pairs of broodstocks revealed that polish and ukrainian broodstocks of european huchen were characterized by the closest genetic distance . 
 in contrast , the highest genetic divergence parameters ( fst and nei s distance ) were observed among german , slovak , and ukrainian broodstocks . 
Introduction
Materials and methods
Sample collection and DNA extraction
Polymerase chain reaction amplification (PCR)
Genotyping
Data analysis
Results
Discussion