Polymorphism of mtDNA regions in brown hare (Lepus europaeus) populations from Vojvodina (Serbia and Montenegro)

It is well known that there is heterogeneity of mitochondrial DNA (mtDNA) within and among natural populations of same species. The polymorphism level of particular regions of mtDNA gives valuable results in detection of population genetic structure. The aim of this paper was to detect polymorphism of three mtDNA regions: cytochrome oxidase I (COI), Control region, and 12S/16S rRNA, by the mtDNA RFLP-PCR method, in three Lepus europaeus populations from Vojvodina province (Serbia and Montenegro). Polymorphism was detected within the two regions, COI and Control region, while 12S/16S rRNA region was monomorphic in all 77 individuals. Eight haplotypes were detected in the brown hare population in Vojvodina, and three were unique for the Srem brown hare population.     

Intraspecific mitochondrial DNA variation in the Colorado potato beetle,Leptinotarsa decemlineata (Coleoptera: Chrysomelidae)

Restriction endonuclease analyses of mitochondrial DNA (mtDNA) were used to examine genetic variability and population structure inLeptinotarsa decemlineata (Say). A group of three enzymes, EcoRI, HpaI, and PstI, was used to reveal polymorphism both within and among some of the 10 populations tested, yielding 16 haplotypes in combination. The frequencies of these 16 haplotypes differed significantly across geographic regions, indicating some partitioning of mtDNA haplotypes. Estimates of mtDNA sequence divergence (δ) between haplotypes ranged from 0.016 to 0.135%, suggesting local differentiation of mtDNA in some populations. Analysis of these data suggests that Texas was colonized by more than one mtDNA lineage, most likely originating in Mexico. We hypothesize that a larger founder size for the initial introductions or high levels of variability in the parent population at the edge of the CPB expanding range led to the initial partitioning of haplotypes observed in samples from Texas.     

The impact of stocking on the genetic structure of European grayling <Emphasis Type="Italic">Thymallus thymallus</Emphasis>, Salmonidae) in two alpine rivers

The European grayling (Thymallus thymallus) is considered to be threatened in several European regions. In recent decades fishery managers have increasingly turned to stocking programs as one way to combat the negative effect that human influence has had on population densities. The present study surveyed the genetic structure of two Danube drainage populations at the Inn and Drau rivers, in Austria, on the basis of mtDNA sequences. Data were placed in the context of the phylogeographic structure of European Thymallus and thus could reveal unexpected geographical mixing due to stocking with allochthonous individuals. Our analyses revealed that regular stocking of fish not originating from their natal Rivers has left genetic traces in both systems surveyed. These traces may be classified as marginal for the Inn river and its tributaries in which 97% of the graylings investigated carried haplotypes belonging to the northern alpine lineage, corresponding to the region through which the Inn flows. In contrast, the genetic composition of the Drau population, situated in the southern Alps, has been seriously altered through the stocking of fish belonging to the northern alpine mtDNA lineage as only 62% of the fishes sampled carried haplotypes representing the native southern alpine lineage.     

DIFFERING LEVELS OF AMONG-POPULATION DIVERGENCE IN THE MITOCHONDRIAL DNA OF PERIODICAL CICADAS RELATED TO HISTORICAL BIOGEOGRAPHY

Mitochondrial DNA (mtDNA) haplotypes were determined for 118 individuals of 13-and 17-year periodical cicadas (genus Magicicada) collected from 16 localities throughout the Midwest and eastern United States. Two distinct mtDNA lineages, identified as A and B, differ by 2.5% based on analysis of fragment patterns and restriction maps. Observed levels of mtDNA diversity within each lineage are low compared to estimates for other taxa. The two lineages are regionally segregated, with the boundary line occurring at a latitude of approximately 33° North. The levels of mtDNA diversity and population genetic structure differ within the two lineages. There is a remarkably low level of mean mtDNA divergence and no genetic structure in lineage A, whereas lineage B exhibits an order of magnitude higher level of mtDNA diversity and significant genetic structure among sampled populations. The low level of mtDNA diversity in cicadas may be attributed to (1) a population bottleneck that most likely occurred during the Pleistocene, (2) recent colonization following the retreat of the glaciers and the expansion of deciduous forests, and/or (3) high among-family reproductive variance (as a consequence of large population size, high fecundity, aggregative behavior of adults, and clumping of eggs). The difference in mtDNA diversity and population genetic structure between the lineages suggests that they experienced different biogeographic histories; we relate this to Pleistocene changes.     

Genetic diversity of Garra rufa Heckel, 1843 (Teleostei: Cyprinidae) in Anatolia

In the present study, mitochondrial DNA polymerase chain reaction-restriction-fragment length polymorphism (PCR-RFLP) and nuclear DNA inter-simple sequence repeat (ISSR) assays were used to assess the phylogenetic and phylogeographic relationships among Garra rufa samples from Anatolia. The complete mtDNA NADH 3/4 dehydrogenase (ND-3/4) gene amplified by PCR was digested with eight restriction enzymes. These enzymes produced 20 composite haplotypes for G. rufa populations. All the mtDNA haplotypes detected were highly diverged from each other and each lineage had a unique genetic profile. The evaluation of mtDNA PCR-RFLP data coupled with geological history of Anatolia indicated a deep genetic divergence among the mtDNA haplotypes of G. rufa populations from drainages of the Mediterranean and Persian Gulf, suggesting an early isolation of Tigris-Euphrates with Orontes river and other rivers draining into the Mediterranean Sea. In general, data from both mtDNA and nDNA were congruent.     

Phylogeographical analysis of Ligia oceanica (Crustacea: Isopoda) reveals two deeply divergent mitochondrial lineages

Isopods of the species Ligia oceanica are typical inhabitants of the rocky intertidal of the northern European coastline. The aim of this study was to assess the genetic structure of this species using mitochondrial and nuclear sequence data. We analysed partial mitochondrial cytochrome c oxidase subunit I (CO1) and 16S rRNA gene sequence data of 161 specimens collected from ten sites ranging from Spain to Norway. For selected specimens, we also sequenced the hypervariable V7 expansion segment of the nuclear 18S rRNA gene as a supplementary marker. Furthermore, we studied the infection rate of all analysed specimens by the alphaproteobacterium Wolbachia. Our analyses revealed two deeply divergent mitochondrial lineages for Ligia oceanica that probably diverged in the late Pliocene to mid Pleistocene. One lineage comprised specimens from northern populations (‘lineage N’) and one primarily those from France and Spain (‘lineage S’). Distribution patterns of the haplotypes and the genetic distances between both lineages revealed two populations that diverged before the Last Glacial Maximum. Given that we found no homogenization of mitochondrial haplotypes, our present results also reject any influence of Wolbachia on the observed mtDNA variability. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 112 , 16–30.     

Patterns of variability at the major histocompatibility class I and class II loci in populations of the endangered cyprinid <Emphasis Type="Italic">Ladigesocypris ghigii</Emphasis>

The patterns of MHC diversity were studied at UAA and DAB1 loci and the two domains involved in the recognition of antigenic peptides (α2 and β1, respectively) in eight Ladigesocypris ghigii populations inhabiting streams and a concrete reservoir, in order to understand the significance of these genes in bottlenecked populations of an endemic species and develop conservation rationale. In agreement with previous study employing RAPD and mtDNA markers (Mamuris et al., Freshw Biol 50:1441–1453, 2005), both loci exhibited a very low level of polymorphism with only two and four alleles detected for UAA and DAB1, respectively. The functional MHC diversity was even lower since UAA alleles were distinguished by a single synonymous substitution. The type of habitat did not affect the level of polymorphism. Our data suggest that DAB1 polymorphism might be the outcome of the positive selection, imposed by the temporal and spatial variation of pathogen load, and the genetic drift as a result of successive habitat shrinkage and deterioration by water abstraction year after year. The populations studied were significantly less diverged at MHC loci than expected based on nuclear and mtDNA markers, suggesting that common parasites might act as causative factors to homogenize selection. Sufficient epidemiological data are required for the interpretation of the results and decision-making on suitable conservation actions.     

Retrospective Estimation of Genetic Diversity of an Extinct Oriental White Stork (Ciconia boyciana) Population in Japan Using Mounted Specimens and Implications for Reintroduction Programs

The Japanese regional population of the Oriental white stork (Ciconia boyciana) became extinct in 1986. Mitochondrial DNA (mtDNA) D-loop region from 20 mounted specimens preserved at public facilities in Toyooka City, Hyogo Prefecture, Japan and its vicinity (n = 17), the area inhabited by the last of the Japanese population, and samples originating from China (n = 3) which were kept at a zoo was analyzed. After extracting DNA from small pieces of skin from mounted specimens, a 1210-bp region of the mtDNA D-loop region was analyzed. The haplotypes among 11 specimens of storks captured or found dead at Toyooka City just before the population became extinct were completely identical. Four haplotypes observed among the mounted specimens preserved in the vicinity of Toyooka City were differentiated from those of captive storks originating from China and Russia in a previous study. Therefore, the last Japanese population might have been a genetically unique group. However, phylogenetic analysis using the maximum likelihood method showed that haplotypes found in the Japanese regional population were closely related to the Chinese and Russian lineages (sequence difference = 2.1%). One mounted specimen collected in 1935 at Izushi village, in the vicinity of Toyooka City, showed the same haplotype as the captive storks from China, suggesting that genetic flow may have historically occurred between the populations of Japan and the continent. Recently, reintroduction for the Oriental white stork has been planned in Toyooka City. The planning for the recovery of extinct populations should not only involve translocation of species to the range from which it disappeared, but also reconstruction of regional populations while considering the genetic lineage between the extinct and introduced populations.     

Phylogeography of northern Dolly Varden Salvelinus malma malma based on analysis of mitochondrial DNA

The northern Dolly Varden, Salvelinus malma malma, is a typical representative of arctic fauna distributed in northeastern Asia and northwestern North America. Because its spawning habitats were affected by Pleistocene glacial advances over most of its natural range, S. m. malma is among the most interesting objects of phylogeographic and microevolutionary studies. We reconstructed the genealogy of mtDNA haplotypes from 27 Alaskan and Asian populations to study the influence of glacial and geological vicariance events on the contemporary population genetic structure, phylogeographic subdivision and distribution of the northern Dolly Varden. Analysis of restriction site states in three PCR‐amplified mtDNA regions (ND1/ND2, ND5/ND6, Cytb/D‐loop; 47% of the mitochondrial genome) resolved 75 haplotypes in 436 fish. Similar patterns of subspecific variation apparent from hierarchical diversity and nested clade analyses of mtDNA haplotypes identify weak spatial differentiation and low levels of divergence. Our results suggest that (1) demographic history has been influenced by historical range expansions and recent isolation by distance, (2) present populations from Asia and North America were colonized from one main Beringian Refugium, and (3) that this taxon's ancestral population probably experienced a bottleneck in the Beringian Refugium during the late Pleistocene (Wisconsin) glacial period.     

Genetic diversity of the Chinese goat in the littoral zone of the Yangtze River as assessed by microsatellite and mtDNA

The objective of this study was to assess the genetic diversity and population structure of goats in the Yangtze River region using microsatellite and mtDNA to better understand the current status of those goat genetic diversity and the effects of natural landscape in fashion of domestic animal genetic diversity. The genetic variability of 16 goat populations in the littoral zone of the Yangtze River was estimated using 21 autosomal microsatellites, which revealed high diversity and genetic population clustering with a dispersed geographical distribution. A phylogenetic analysis of the mitochondrial D‐loop region (482 bp) was conducted in 494 goats from the Yangtze River region. In total, 117 SNPs were reconstructed, and 173 haplotypes were identified, 94.5% of which belonged to lineages A and B. Lineages C, D, and G had lower frequencies (5.2%), and lineage F haplotypes were undetected. Several high‐frequency haplotypes were shared by different ecogeographically distributed populations, and the close phylogenetic relationships among certain low‐frequency haplotypes indicated the historical exchange of genetic material among these populations. In particular, the lineage G haplotype suggests that some west Asian goat genetic material may have been transferred to China via Muslim migration.     

MITOCHONDRIAL DNA VARIATION IN THE FIELD VOLE (MICROTUS AGRESTIS): REGIONAL POPULATION STRUCTURE AND COLONIZATION HISTORY

Restriction fragment length polymorphism analysis of mitochondrial DNA (mtDNA) was used to examine the genetic structure among field voles (Microtus agrestis) from southern and central Sweden. A total of 57 haplotypes was identified in 158 voles from 60 localities. Overall mtDNA diversity was high, but both haplotype and nucleotide diversity exhibited pronounced geographic heterogeneity. Phylogenetic analyses revealed a shallow tree with seven primary mtDNA lineages separated by sequence divergences ranging from 0.6% to 1.0%. The geographic structure of mtDNA diversity and lineage distribution was complex but strongly structured and deviated significantly from an equilibrium situation. The extensive mtDNA diversity observed and the recent biogeographic history of the region suggests that the shallow mtDNA structure in the field vole cannot be explained solely by stochastic lineage sorting in situ or isolation by distance. Instead, the data suggest that the genetic imprints of historical demographic conditions and vicariant geographic events have been preserved and to a large extent determine the contemporary geographic distribution of mtDNA variation. A plausible historical scenario involves differentiation of mtDNA lineages in local populations in glacial refugia, a moving postglacial population structure, and bottlenecks and expansions of mtDNA lineages during the postglacial recolonization of Sweden. By combining the mtDNA data with an analysis of Y-chromosome variation, a specific population unit was identified in southwestern Sweden. This population, defined by a unique mtDNA lineage and fixation of a Y-chromosome variant, probably originated in a population bottleneck in southern Sweden about 12,000 to 13,000 calendar years ago.     

Within‐population genetic effects of mtDNA on metabolic rate in Drosophila subobscura

A growing body of research supports the view that within‐species sequence variation in the mitochondrial genome (mtDNA) is functional, in the sense that it has important phenotypic effects. However, most of this empirical foundation is based on comparisons across populations, and few studies have addressed the functional significance of mtDNA polymorphism within populations. Here, using mitonuclear introgression lines, we assess differences in whole‐organism metabolic rate of adult Drosophila subobscura fruit flies carrying either of three different sympatric mtDNA haplotypes. We document sizeable, up to 20%, differences in metabolic rate across these mtDNA haplotypes. Further, these mtDNA effects are to some extent sex specific. We found no significant nuclear or mitonuclear genetic effects on metabolic rate, consistent with a low degree of linkage disequilibrium between mitochondrial and nuclear genes within populations. The fact that mtDNA haplotype variation within a natural population affects metabolic rate, which is a key physiological trait with important effects on life‐history traits, adds weight to the emergent view that mtDNA haplotype variation is under natural selection and it revitalizes the question as to what processes act to maintain functional mtDNA polymorphism within populations.     

Conservation implications of the evolutionary history and genetic diversity hotspots of the snowshoe hare

With climate warming, the ranges of many boreal species are expected to shift northward and to fragment in southern peripheral ranges. To understand the conservation implications of losing southern populations, we examined range‐wide genetic diversity of the snowshoe hare (Lepus americanus), an important prey species that drives boreal ecosystem dynamics. We analysed microsatellite (8 loci) and mitochondrial DNA sequence (cytochrome b and control region) variation in almost 1000 snowshoe hares. A hierarchical structure analysis of the microsatellite data suggests initial subdivision in two groups, Boreal and southwestern. The southwestern group further splits into Greater Pacific Northwest and U.S. Rockies. The genealogical information retrieved from mtDNA is congruent with the three highly differentiated and divergent groups of snowshoe hares. These groups can correspond with evolutionarily significant units that might have evolved in separate refugia south and east of the Pleistocene ice sheets. Genetic diversity was highest at mid‐latitudes of the species' range, and genetic uniqueness was greatest in southern populations, consistent with substructuring inferred from both mtDNA and microsatellite analyses at finer levels of analysis. Surprisingly, snowshoe hares in the Greater Pacific Northwest mtDNA lineage were more closely related to black‐tailed jackrabbits (Lepus californicus) than to other snowshoe hares, which may result from secondary introgression or shared ancestral polymorphism. Given the genetic distinctiveness of southern populations and minimal gene flow with their northern neighbours, fragmentation and loss of southern boreal habitats could mean loss of many unique alleles and reduced evolutionary potential.     

Mitochondrial DNA diversity, population structure, and gender association in the gynodioecious plant Silene vulgaris

A highly variable mitochondrial DNA (mtDNA) restriction fragment length polymorphism (RFLP) locus is used to assess the population structure of mitochondrial genomes in the gynodioecious plant Silene vulgaris at two spatial scales. Thirteen mtDNA haplotypes were identified within 250 individuals from 18 populations in a 20-km diameter region of western Virginia. The population structure of these mtDNA haplotypes was estimated as thetaST = 0.574 (+/- 0.066 SE) and, surprisingly, genetic differentiation among populations was negatively correlated with geographic distance (Mantel r = -0.246, P < 0.002). Additionally, mtDNA haplotypes were spatially clumped at the scale of meters within one population. Gender in S. vulgaris is determined by an interaction between autosomal male fertility restorers and cytoplasmic male sterility (CMS) factors, and seed fitness is affected by an interaction between gender and population sex ratio; thus, selection acting on gender could influence the distribution of mtDNA RFLP haplotypes. The sex ratio (females:hermaphrodites) varied among mtDNA haplotypes across the entire metapopulation, possibly because the haplotypes were in linkage disequilibrium with different CMS factors. The gender associated with some of the most common haplotypes varied among populations, suggesting that there is also population structure in male fertility restorer genes. In comparison with reports of mtDNA variation from other published studies, we found that S. vulgaris exhibits a large number of mtDNA haplotypes relative to that observed in other species.     

Evaluation of genetic introgression from domesticated pigs into the Ryukyu wild boar population on Iriomote Island in Japan

We evaluated genetic introgression from domesticated pigs into the Ryukyu wild boar (RWB) population on Iriomote Island based on their genetic structure and diversity. We used a combination of mitochondrial DNA D‐loop region (596 bp) polymorphisms and 23 microsatellite markers. RWBs (n = 130) were collected from 18 locations on Iriomote Island and compared with 66 reference samples of European and Asian domestic pigs. We identified six distinct haplotypes, involving 22 single nucleotide polymorphisms (including one insertion) in the RWB population. The phylogenetic tree had two branches: the RWB group and domestic lineage. Fourteen of 130 RWBs (10.8%) belonged to the European domestic lineage, including 11 RWBs from the Panari Islands, northwest of Iriomote Main Island (IMI). The heterozygosity values, total number of alleles, number of effective alleles and polymorphism information content of the RWB groups were lower than those of the European domestic groups. The RWB population on IMI had a lower heterozygous deficiency index (F_IS = 0.059) than did the other populations, which indicates that this population was more inbred. There was a large genetic distance (F_ST = 0.560) between RWBs on IMI and the Meishan populations. Structure analysis using the 23 microsatellite markers revealed that 16 RWBs had an admixture pattern between RWB and domesticated pig breeds. These results suggest that gene flow may have occurred from domestic pigs to RWBs and demonstrate that there was low genetic variation in the IMI population.     

Geographic differentiation of <Emphasis Type="Italic">Cobitis shikokuensis</Emphasis> inferred from mtDNA RFLP analysis

Mitochondrial DNA divergence among populations of the Japanese spinous loach Cobitis shikokuensis, endemic to Shikoku Island, was investigated by restricted fragment length polymorphism analysis. A total of 68 restriction sites on DNA fragments from the cytochrome b to D-loop regions and from the 12S rRNA to 16S rRNA regions, amplified by PCR, were analyzed. A total of 12 haplotypes (plus 6 in outgroups) were detected in 268 specimens collected from 19 localities in seven rivers (and 41 specimens from four localities in three rivers in outgroups). Three of the seven river populations of C. shikokuensis were shown to have unique haplotypes, and four of the seven river populations were monomorphic. The nested structure of the haplotype network for populations of C. shikokuensis exhibited two large clades corresponding to (1) populations from the Shimanto River and its neighbors and (2) two genetically divergent populations in the Shigenobu and Iwamatsu Rivers. The population from the Shimanto River, the largest river inhabited by C. shikokuensis, maintains great haplotype diversity as well as the allozyme diversity previously reported. On the other hand, populations from the Hiji River, the second largest river, which exhibited the highest allozyme diversity, were monomorphic in their mtDNA. The nested clade analysis (NCA) revealed that past fragmentation between the above two clades could occur in the initial distribution process of C. shikokuensis. The large genetic divergence of two river populations from the Shigenobu and Iwamatsu Rivers was inferred to be caused by a process of long distance colonization and fragmentation. MtDNA introgression into the Hiji River population from southern river populations was suggested. Taking genetic divergence into consideration, each river population of C. shikokuensis should be conserved separately as like a distinct species, and conservation programs for the small populations showing less genetic variability should be invoked as soon as possible.     

Comparison between mitochondrial and chloroplast DNA variation in the native range of Silene vulgaris

A detailed survey of mitochondrial and chloroplast diversity in eight populations of Silene vulgaris from Central Europe was conducted for comparison with previously published data on diversity from S. vulgaris populations in the introduced range. Mitochondrial DNA (mtDNA) variation around the coxI gene was assessed with Southern blotting/restriction fragment length polymorphism methods. Chloroplast variation was assessed by sequencing the intergenic spacer separating the trnH and psbA genes. Thirty mtDNA haplotypes and 24 chloroplast DNA (cpDNA) haplotypes were found within 86 individuals. The overall genetic diversity h (0.941 for mitochondrial, and 0.893 for chloroplast markers) and within-population diversity were higher than reported in previous population studies of S. vulgaris in the USA and Europe. The frequency of private alleles was surprisingly high - more than 90% for both kinds of markers. Most of our populations were large and located in relatively undisturbed meadows, whereas surveys in Virginia consisted of smaller roadside populations. The slow rate of population turnover in European populations is discussed as a factor responsible for the relatively high diversity of S. vulgaris in undisturbed areas of its native range. Association between mtDNA and cpDNA haplotypes was also demonstrated. Finally, gender and mtDNA haplotype were associated in the Alps populations, where females were very rare.     

Nucleotide diversity of a ND5 fragment confirms that population expansion is the most suitable explanation for the mtDNA haplotype polymorphism of <Emphasis Type="Italic">Drosophila subobscura</Emphasis>

Results from mitochondria (mt) DNA restriction site analyses (RSAs) have revealed that wild populations of Drosophila subobscura are formed by two common (I and II) and some rare, often endemic, low-frequency haplotypes. In the study reported here, we analysed nucleotide diversity in a 942-bp fragment of the mtDNA ND5 gene in 48 D. subobscura individuals captured from three populations that showed haplotypes I, II or the less common ones, as well as in one additional individual belonging to D. guanche that was taken as an outgroup. RSAs and sequencing results were compared. The two approaches yielded similar nucleotide variability parameters, suggesting a consistency in the results obtained from mtDNA dynamics in natural populations of D. subobscura. Patterns of polymorphism at ND5 are most consistent with the hypothesis of population expansion after a bottleneck that may have occurred since the last glaciation or which may occur seasonally after the summer and winter. However, we cannot rule out that selection has a role in maintaining the two major haplotypes at intermediate frequencies in worldwide populations of D. subobscura.     

The Amazon River system as an ecological barrier driving genetic differentiation of the pink dolphin (Inia geoffrensis)

The pink dolphin (Inia geoffrensis) is widely distributed along the Amazon and Orinoco basins, covering an area of approximately 7 million km². Previous morphological and genetic studies have proposed the existence of at least two evolutionary significant units: one distributed across the Orinoco and Amazon basins and another confined to the Bolivian Amazon. The presence of barriers in the riverine environment has been suggested to play a significant role in shaping present‐day patterns of ecological and genetic structure for this species. In the present study, we examined the phylogeographic structure, lineage divergence time and historical demography using mitochondrial (mt)DNA sequences in different pink dolphin populations distributed in large and small spatial scales, including two neighbouring Brazilian Amazon populations. mtDNA control region (CR) analysis revealed that the Brazilian haplotypes occupy an intermediate position compared to three previously studied geographic locations: the Colombian Amazon, the Colombian Orinoco, and the Bolivian Amazon. On a local scale, we have identified a pattern of maternal isolation between two neighbouring populations from Brazil. Six mtDNA CR haplotypes were identified in Brazil with no sharing between the two populations, as well as specific cytochrome b (cyt b) haplotypes identified in each locality. In addition, we analyzed autosomal microsatellites to investigate male‐mediated gene flow and demographic changes within the study area in Brazil. Data analysis of 14 microsatellite loci failed to detect significant population subdivision, suggesting that male‐mediated gene flow may maintain homogeneity between these two locations. Moreover, both mtDNA and microsatellite data indicate a major demographic collapse within Brazil in the late Pleistocene. Bayesian skyline plots (BSP) of mtDNA data revealed a stable population for Colombian and Brazilian Amazon lineages through time, whereas a population decline was demonstrated in the Colombian Orinoco lineage. Moreover, BSP and Tajima's D and Fu's F_s tests revealed a recent population expansion exclusively in the Bolivian sample. Finally, we estimated that the diversification of the Inia sp. lineage began in the Late Pliocene (approximately 3.1 Mya) and continued throughout the Pleistocene. Published 2011. This article is a US Government work and is in the public domain in the USA. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 812–827.     

Genetic variation of European grayling (Thymallus thymallus) populations in the Western Balkans

In order to elucidate genetic composition of European grayling (Thymallus thymallus) populations in the Western Balkans, the partial mitochondrial DNA (mtDNA) control region was sequenced and 12 microsatellite loci genotyped in 14 populations originating from tributaries of the Adriatic and Danube drainages. Eleven mtDNA haplotypes were found, one confined to the Adriatic clade, one to the Alpine group and the rest to the ‘Balkan’ grayling phylogenetic clade. Haplotypes from the Balkan clade were confined to the Danube drainage and constituted two groups: northern group with haplotypes found in the Slovenian part of the Danube drainage, and southern group, consisting from Bosnia–Herzegovina and Montenegro. Substantial genetic distance between northern and southern groups of haplotypes (0.75–1.8%) and well supported divisions within the northern group indicate very structured grayling population within the studied Danube basin that most probably did not evolve due to vicariance but rather as a consequence of multiple colonization waves that might have occurred during the Pleistocene. Furthermore, genetic distance of ~4% between Adriatic and Danube populations’ haplotypes, suggest that their separation occurred in mid-Pliocene. These findings imply a complex colonization pattern of the Western Balkans drainages. Microsatellite data also confirm high genetic diversity in Western Balkans populations of grayling (on average 7.5 alleles per microsatellite locus and H _exp 0.58). Limited stocking activities were detected based on microsatellites and mtDNA data. Regarding current knowledge of grayling phylogeography appropriate management strategies were proposed to preserve unique, autochthonous grayling populations in Western Balkan.