Evidence for transoceanic migrations by loggerhead sea turtles in the southern Pacific Ocean

Post-hatchling loggerhead turtles (Caretta caretta) in the northern Pacific and northern Atlantic Oceans undertake transoceanic developmental migrations. Similar migratory behaviour is hypothesized in the South Pacific Ocean as post-hatchling loggerhead turtles are observed in Peruvian fisheries, yet no loggerhead rookeries occur along the coast of South America. This hypothesis was supported by analyses of the size-class distribution of 123 post-hatchling turtles in the South Pacific and genetic analysis of mtDNA haplotypes of 103 nesting females in the southwest Pacific, 19 post-hatchlings stranded on the southeastern Australian beaches and 22 post-hatchlings caught by Peruvian longline fisheries. Only two haplotypes (CCP1 93% and CCP5 7%) were observed across all samples, and there were no significant differences in haplotype frequencies between the southwest Pacific rookeries and the post-hatchlings. By contrast, the predominant CCP1 haplotype is rarely observed in North Pacific rookeries and haplotype frequencies were strongly differentiated between the two regions (F_st=0.82; p=<0.00001). These results suggest that post-hatchling loggerhead turtles emerging from the southwest Pacific rookeries are undertaking transoceanic migrations to the southeastern Pacific Ocean, thus emphasizing the need for a broader focus on juvenile mortality throughout the South Pacific to develop effective conservation strategies.     

The genetic structure of the loggerhead sea turtle (Caretta caretta) in the Mediterranean as revealed by nuclear and mitochondrial DNA and its conservation implications

The population genetic structure of the loggerhead sea turtle (Caretta caretta) nesting in the eastern Mediterranean was assessed by sequencing a fragment of the control region of the mitochondrial DNA (n = 190) and seven microsatellites (n = 112). The two types of markers revealed genetic structuring (mtDNA: γ_st = 0.212, P < 0.001; nDNA F _st = 0.006, P < 0.001), thus indicating that both females and males are philopatric and that gene flow between populations is restricted. Mitochondrial DNA data indicate that the female populations nesting on the islands of Crete and Cyprus have suffered a recent bottleneck or colonization event. However, no bottleneck or founder effect was revealed by nuclear markers, thus indicating male-mediated gene flow from other populations that would increase nuclear genetic variability. Crete, and to a lower extent Cyprus, are thought to play a central role in such male-mediated gene flow that may reduce the negative effect of genetic drift or inbreeding on the small populations of Lebanon and Israel. This population structure indicates that assessing population relevance only on the basis of genetic variability and size would be misleading, as some populations not fulfilling those requirements may play a relevant role in genetic exchange and hence contribute to the overall genetic variability.     

Genetic structure of Octopus vulgaris (Cephalopoda,Octopodidae) in the central Mediterranean Sea inferred from the mitochondrial COIII gene

The polymorphism of the mitochondrial gene cytochrome oxidase III was studied in the Mediterranean octopus, Octopus vulgaris Cuvier, 1797. A total of 202 specimens from seven sampling sites were analysed with the aim of elucidating patterns of genetic structure in the central Mediterranean Sea and to give an insight into the phylogeny of the Octopus genus. Phylogenetic analyses showed that individuals from the central Mediterranean belong to the O. vulgaris species whose limits should nevertheless be clarified. Concerning genetic structure, two high-frequency haplotypes were present in all locations. The overall genetic divergence (Φ_ST = 0.05, P < 0.05) indicated a significant genetic structuring in the study area and an AMOVA highlighted a significant break between western and eastern Mediterranean basins (Φ_CT = 0.094, P < 0.05). Possible explanations for the observed patterns of genetic structuring are discussed with reference to their relevance for fisheries management.     

Mitochondrial DNA Diversity of Orchid Bee Euglossa fimbriata (Hymenoptera: Apidae) Populations Assessed by PCR-RFLP

Euglossa fimbriata is a euglossine species widely distributed in Brazil and occurring primarily in Atlantic Forest remnants. In this study, the genetic mitochondrial structure of E. fimbriata from six Atlantic Forest fragments was studied by RFLP analysis of three PCR-amplified mtDNA gene segments (16S, COI-COII, and cyt b). Ten composite haplotypes were identified, six of which were exclusive and represented singleton mitotypes. Low haplotype diversity (0.085–0.289) and nucleotide diversity (0.000–0.002) were detected within samples. AMOVA partitioned 91.13% of the overall genetic variation within samples and 8.87% (?_st = 0.089; P < 0.05) among samples. Pairwise comparisons indicated high levels of differentiation among some pairs of samples (?_st = 0.161–0.218; P < 0.05). These high levels indicate that these populations of E. fimbriata, despite their highly fragmented landscape, apparently have not suffered loss of genetic variation, suggesting that this particular population is not currently endangered.     

Extensive hybridization in hawksbill turtles (Eretmochelys imbricata) nesting in Brazil revealed by mtDNA analyses

Bahia state hosts over 90% of hawksbill (Eretmochelys imbricata) nests registered in the main nesting sites monitored by Projeto Tamar-IBAMA in Brazil. The genetic diversity of this hawksbill population (n=119) was assayed through the analyses of 752 bp of the mitochondrial DNA control region in nesting females. Seven distinct haplotypes, defined by 125 polymorphic sites, were found. Most of the individuals (n=67) display four typical hawksbill haplotypes, 50 individuals display two haplotypes characteristic of the loggerhead turtle (Caretta caretta) and two individuals had a haplotype affiliated with the olive ridley (Lepidochelys olivacea). These results demonstrate hybridization between the hawksbills and two species that nest along the Bahia coast. Of special interest is the high occurrence of loggerhead × hawksbill hybrids (42%), which display loggerhead mtDNA haplotypes but are characterized morphologically as hawksbills. The true hawksbill haplotypes present only three variable sites and low genetic diversity values (h=0.358±0.069; π=0.0005±0.0001). The occurrence of several nesting individuals with identical mtDNA from another species may also suggest a long history of introgression between species producing likely F2 or further generation hybrids. Marine turtle hybrids have been previously reported, but the high frequency observed in Bahia is unprecedented. Such introgression may influence evolutionary pathways for all three species, or may introduce novel morphotypes that develop apart from the parental species. The presence of a unique hybrid swarm has profound conservation implications and will significantly influence the development and implementation of appropriate management strategies for these species.     

GLOBAL PHYLOGEOGRAPHY OF THE LOGGERHEAD TURTLE (CARETTA CARETTA) AS INDICATED BY MITOCHONDRIAL DNA HAPLOTYPES

Restriction-site analyses of mitochondrial DNA (mtDNA) from the loggerhead sea turtle (Caretta caretta) reveal substantial phylogeographic structure among major nesting populations in the Atlantic, Indian, and Pacific oceans and the Mediterranean sea. Based on 176 samples from eight nesting populations, most breeding colonies were distinguished from other assayed nesting locations by diagnostic and often fixed restriction-site differences, indicating a strong propensity for natal homing by nesting females. Phylogenetic analyses revealed two distinctive matrilines in the loggerhead turtle that differ by a mean estimated sequence divergence p = 0.009, a value similar in magnitude to the deepest intraspecific mtDNA node (p = 0.007) reported in a global survey of the green sea turtle Chelonia mydas. In contrast to the green turtle, where a fundamental phylogenetic split distinguished turtles in the Atlantic Ocean and the Mediterranean Sea from those in the Indian and Pacific oceans, genotypes representing the two primary loggerhead mtDNA lineages were observed in both Atlantic–Mediterranean and Indian-Pacific samples. We attribute this aspect of phylogeographic structure in Caretta caretta to recent interoceanic gene flow, probably mediated by the ability of this temperate-adapted species to utilize habitats around southern Africa. These results demonstrate how differences in the ecology and geographic ranges of marine turtle species can influence their comparative global population structures.     

Population stock structure of leatherback turtles (Dermochelys coriacea) in the Atlantic revealed using mtDNA and microsatellite markers

This study presents a comprehensive genetic analysis of stock structure for leatherback turtles (Dermochelys coriacea), combining 17 microsatellite loci and 763 bp of the mtDNA control region. Recently discovered eastern Atlantic nesting populations of this critically endangered species were absent in a previous survey that found little ocean-wide mtDNA variation. We added rookeries in West Africa and Brazil and generated longer sequences for previously analyzed samples. A total of 1,417 individuals were sampled from nine nesting sites in the Atlantic and SW Indian Ocean. We detected additional mtDNA variation with the longer sequences, identifying ten polymorphic sites that resolved a total of ten haplotypes, including three new variants of haplotypes previously described by shorter sequences. Population differentiation was substantial between all but two adjacent rookery pairs, and F _ST values ranged from 0.034 to 0.676 and 0.004 to 0.205 for mtDNA and microsatellite data respectively, suggesting that male-mediated gene flow is not as widespread as previously assumed. We detected weak (F _ST = 0.008 and 0.006) but significant differentiation with microsatellites between the two population pairs that were indistinguishable with mtDNA data. POWSIM analysis showed that our mtDNA marker had very low statistical power to detect weak structure (F _ST < 0.005), while our microsatellite marker array had high power. We conclude that the weak differentiation detected with microsatellites reflects a fine scale level of demographic independence that warrants recognition, and that all nine of the nesting colonies should be considered as demographically independent populations for conservation. Our findings illustrate the importance of evaluating the power of specific genetic markers to detect structure in order to correctly identify the appropriate population units to conserve.     

Large‐scale connectivity,cryptic population structure,and relatedness in Eastern Pacific Olive ridley sea turtles (Lepidochelys olivacea)

Endangered species are grouped into genetically discrete populations to direct conservation efforts. Mitochondrial control region (mtCR) haplotypes are used to elucidate deep divergences between populations, as compared to nuclear microsatellites that can detect recent structuring. When prior populations are unknown, it is useful to subject microsatellite data to clustering and/or ordination population inference. Olive ridley sea turtles (Lepidochelys olivacea) are the most abundant sea turtle, yet few studies have characterized olive ridley population structure. Recently, clustering results of olive ridleys in the Eastern Tropical Pacific Ocean suggested weak structuring (F_ST = 0.02) between Mexico and Central America. We analyzed mtCR haplotypes, new microsatellite genotypes from Costa Rica, and preexisting microsatellite genotypes from olive ridleys across the Eastern Tropical Pacific, to further explore population structuring in this region. We subjected inferred populations to multiple analyses to explore the mechanisms behind their structuring. We found 10 mtCR haplotypes from 60 turtles nesting at three sites in Costa Rica, but did not detect divergence between Costa Rican sites, or between Central America and Mexico. In Costa Rica, clustering suggested one population with no structuring, but ordination suggested four cryptic clusters with moderate structuring (F_ST = 0.08, p < .001). Across the Eastern Tropical Pacific, ordination suggested nine cryptic clusters with moderate structuring (F_ST = 0.103, p < .001) that largely corresponded to Mexican and Central American populations. All ordination clusters displayed significant internal relatedness relative to global relatedness (p < .001) and contained numerous sibling pairs. This suggests that broadly dispersed family lineages have proliferated in Eastern Tropical Pacific olive ridleys and corroborates previous work showing basin‐wide connectivity and shallow population structure in this region. The existence of broadly dispersed kin in Eastern Tropical Pacific olive ridleys has implications for management of olive ridleys in this region, and adds to our understanding of sea turtle ecology and life history, particularly in light of the natal‐homing paradigm.     

Genetic diversity of silver pomfret (Pampus argenteus) populations from the China Sea based on mitochondrial DNA control region sequences

The genetic diversity of silver pomfret (Pampus argenteus) from the Bohai, East China, and South China Seas was investigated using mitochondrial DNA (mtDNA) control region sequence data. We found high levels of variation with 17 haplotypes among the 45 individuals (h = 0.88, π = 0.006). AMOVA analysis detected significant structuring among China Sea silver pomfret (P < 0.05, F_ST = 0.050), which indicates significant genetic differentiation. No significant differentiation between Bohai Sea and East China Sea samples was detected (P > 0.05). Whether silver pomfret from the South China Sea exhibit a different demographic history than those from the Bohai and East China Seas remains to be determined.     

Molecular evidence of male-biased dispersal in loggerhead turtle juveniles

Serum testosterone levels and mtDNA haplotypes were obtained from 65 juvenile loggerhead turtles (Caretta caretta, L.) incidentally caught in the central Mediterranean. The group of specimens carrying a haplotype specific for the northwest Atlantic had higher testosterone levels, and so included more males, than the other one. Since primary sex ratios of northwest Atlantic colonies are strongly skewed towards females, results indicate a male bias among Atlantic turtles entering the Mediterranean. This demonstrates for the first time a sex-biased dispersal of specimens in the pelagic phase, an important factor to be considered in conservation programs.     

Population structure and conservation implications for the loggerhead sea turtle of the Cape Verde Islands

The Cape Verde Islands harbour the second largest nesting aggregation of the globally endangered loggerhead sea turtle in the Atlantic. To characterize the unknown genetic structure, connectivity, and demographic history of this population, we sequenced a segment of the mitochondrial (mt) DNA control region (380 bp, n = 186) and genotyped 12 microsatellite loci (n = 128) in females nesting at three islands of Cape Verde. No genetic differentiation in either haplotype or allele frequencies was found among the islands (mtDNA F _ST = 0.001, P > 0.02; nDNA F _ST = 0.001, P > 0.126). However, population pairwise comparisons of the mtDNA data revealed significant differences between Cape Verde and all previously sequenced Atlantic and Mediterranean rookeries (F _ST = 0.745; P < 0.000). Results of a mixed stock analysis of mtDNA data from 10 published oceanic feeding grounds showed that feeding grounds of the Madeira, Azores, and the Canary Islands, in the Atlantic Ocean, and Gimnesies, Pitiüses, and Andalusia, in the Mediterranean sea, are feeding grounds used by turtles born in Cape Verde, but that about 43% (±19%) of Cape Verde juveniles disperse to unknown areas. In a subset of samples (n = 145) we evaluated the utility of a longer segment (~760 bp) amplified by recently designed mtDNA control region primers for assessing the genetic structure of Atlantic loggerhead turtles. The analysis of the longer fragment revealed more variants overall than in the shorter segments. The genetic data presented here are likely to improve assignment and population genetic analyses, with significant conservation and research applications.     

Evidence of olive ridley mitochondrial genome introgression into loggerhead turtle rookeries of Sergipe, Brazil

The coastline of Sergipe state hosts the main Brazilian nesting sites of Lepidochelys olivacea (Eschscholtz, 1829). The second most abundant species of turtles in Sergipe is Caretta caretta (Linnaeus, 1758). Both sea turtle species, respectively known as olive ridley and loggerhead, are currently listed as endangered by the International Union for the Conservation of Nature and Natural Resources. The genetic diversity of the Sergipe loggerhead population (N = 51) was assayed by analyzing 627 bp from the control region of mitochondrial DNA in nesting females. Three haplotypes were identified: CC-A4, CC-A24 and CC × LO. The last one was recorded for specimens considered hybrids because they represent L. olivacea’s mtDNA, but had the external morphology of C. caretta or of a mixture of both species. Based on the two types of hybrids, it was hypothesized that at least two hybridization events had occurred: a more ancient hybridization event, accompanied by introgression (F2 or later backcrosses), and a recent one (F1), both of which involving the same L. olivacea haplotype. The incidence of L. olivacea mitochondrial genome introgression into the C. caretta rookeries was only observed in Sergipe, which could be related to the large numbers of L. olivacea in this region and an overlap of reproduction periods and distribution areas of both species. This may also be associated to global warming since it might alter the sex ratio of sea turtles, thus facilitating interspecific mating. Awareness of gene flow between these species will significantly influence the development and implementation of adequate management strategies.     

Population Genetic Structure of the Yangtze Finless Porpoise (Neophocaena phocaenoides asiaeorientalis): Implications for Management and Conservation

Understanding the population genetic structure is a prerequisite for conservation of a species. The degree of genetic variability characteristic of the mitochondrial DNA control region has been widely exploited in studies of population genetic structure and can be useful in identifying meaningful population subdivisions. To estimate the genetic profile of the Yangtze finless porpoise (Neophocaena phocaenoides asiaeorientalis), an endangered freshwater population endemic to China, the complete mtDNA control region was examined in 39 individuals belonging to seven different stocks inhabiting the middle and lower reaches of the Yangtze River. Very low genetic diversity was found (nucleotide diversity 0.0011± 0.0002 and haplotypic diversity 0.65± 0.05). The mtDNA genetic pattern of the Yangtze population appears to indicate a founder event in its evolutionary history and to support the marine origin for this population. Analyses by F_st and Φ_st yielded statistically significant population genetic structure (F_st = 0.44, P < 0.05; Φ_st = 0.36, P < 0.05). These results may have significant implications for the management and conservation of the Yangtze finless porpoise in the future.     

Genetic variation and clonal diversity in populations of Nelumbo nucifera (Nelumbonaceae) in central China detected by ISSR markers

To obtain accurate estimates of population structure for purposes of conservation planning for wild lotus (Nelumbo nucifera Gaertn.) in central China, genetic diversity among and within six populations, and clonal diversity within another two populations of the species were analyzed. The genetic diversity was high (percentage of polymorphic bands, PPB = 90.0%; Shannon's information index, I = 0.383 ± 0.234) at the species level, but low within individual study populations (PPB = 35.8%; Shannon's information index I = 0.165 ± 0.241). The mean coefficient of gene differentiation (G_st) was 0.570, indicating that 43.0% of the genetic diversity resided within the population. Analysis of molecular variance (AMOVA) indicated that 50.47% of the genetic diversity among the study populations was attributed to geographical location while 12.3% was attributed to differences in their habitats. An overall value of mean estimated number of gene flow (N_m = 0.377) indicated that there was limited gene flow among the sampled populations. The level of clonal diversity found within the populations was considerably high (Simpson's diversity index, D = 0.985) indicating that clonal diversity contributes to a major extent to the overall genetic variation in the genetic structure of N. nucifera. On the basis of the high G_st and D values detected in this study we recommend that any future conservation plans for this species should be specifically designed to include those representative populations with the highest genetic variation for both in situ conservation and germplasm collection expeditions.     

Genetic divergence between two phenotypically distinct bottlenose dolphin ecotypes suggests separate evolutionary trajectories

Due to their worldwide distribution and occupancy of different types of environments, bottlenose dolphins display considerable morphological variation. Despite limited understanding about the taxonomic identity of such forms and connectivity among them at global scale, coastal (or inshore) and offshore (or oceanic) ecotypes have been widely recognized in several ocean regions. In the Southwest Atlantic Ocean (SWA), however, there are scarce records of bottlenose dolphins differing in external morphology according to habitat preferences that resemble the coastal‐offshore pattern observed elsewhere. The main aim of this study was to analyze the genetic variability, and test for population structure between coastal (n = 127) and offshore (n = 45) bottlenose dolphins sampled in the SWA to assess whether their external morphological distinction is consistent with genetic differentiation. We used a combination of mtDNA control region sequences and microsatellite genotypes to infer population structure and levels of genetic diversity. Our results from both molecular marker types were congruent and revealed strong levels of structuring (microsatellites F_ST = 0.385, p < .001; mtDNA F_ST =  0.183, p < .001; Φ_ST = 0.385, p < .001) and much lower genetic diversity in the coastal than the offshore ecotype, supporting patterns found in previous studies elsewhere. Despite the opportunity for gene flow in potential “contact zones”, we found minimal current and historical connectivity between ecotypes, suggesting they are following discrete evolutionary trajectories. Based on our molecular findings, which seem to be consistent with morphological differentiations recently described for bottlenose dolphins in our study area, we recommend recognizing the offshore bottlenose dolphin ecotype as an additional Evolutionarily Significant Unit (ESU) in the SWA. Implications of these results for the conservation of bottlenose dolphins in SWA are also discussed.     

Genetic diversity and population structure analysis of Qinghai-Tibetan plateau schizothoracine fish (Gymnocypris dobula) based on mtDNA D-loop sequences

Gymnocypris dobula is a commercially important fishery species and mostly distributed in the freshwater of Tibetan plateau. In this study, genetic diversity and intraspecific population differentiation were examined by using mitochondrial DNA D-loop sequences in 97 individuals sampled from three localities (Pali, Lasa and Yanghu) in the Tibetan plateau. Two hundred and fourteen polymorphic sites and 50 haplotypes were defined among the three localities. Genetic diversity analysis showed that the highest genetic diversity level was found in Pali population. Phylogenetic relationships analysis results indicated that closer phylogenetic relationships were found between the Yanghu and Lasa populations. Genetic population differentiation analysis indicated that the majority of variation (84.91%) was attributed to variations among populations and the largest differentiation was found between Pali and Lasa localities (F_st = 0.874, Nm = 0.036). Moreover, the historical demographic events were assessed by implementing the mismatch distribution analysis, Fajima's D test and Fu's Fs test. The results indicated that the Pali population had undergone a demographic expansion, possibly within the last 0.163 MYA (Million Years Ago). Our study firstly identified the population genetic structure of the G. dobula, which could be helpful for artificial breeding, fishery stock identification and resource conservation for this species.     

Assessing the genetic diversity and population structure of Culter alburnus in China based on mitochondrial 16S rRNA and COI gene sequences

The genetic diversity and population genetic structure of Culter alburnus were investigated using mitochondrial cytochrome c oxidase subunit I (COI) and ribosomal 16S subunit (16S rRNA) gene sequences. A total of 89 individuals from four localities were included in the analysis. Overall, 12 polymorphic sites were observed and 10 haplotypes were defined. The C. alburnus populations were characterized by high haplotype diversity (0.587 ± 0.047) and low nucleotide diversity (0.00197 ± 0.00073). Pairwise fixation index (F_ST) analysis indicated significant genetic differentiation among different populations. The hierarchical analysis of molecular variance (AMOVA) analysis also showed significant genetic divergence (φ_ST = 0.3792, P < 0.01) among these populations. The present results suggest that subdivisions exist among four C. alburnus populations, and should be considered as different management unit for effective conservation and management purposes. This study provides new information for genetic assessment and will be crucial for establishing fisheries management and strategies for this species.     

Genetic and morphological analyses of tub gurnard Chelidonichthys lucerna populations in Turkish marine waters

Genetic and morphological structure of tub gurnard Chelidonichthys lucerna populations in Turkish marine waters were investigated with mtDNA sequencing of 16S rRNA and morphological characters. C. lucerna samples were collected from the Black Sea, Marmara, Aegean and northeastern Mediterranean coasts of Turkey. The lowest genetic diversity was found in the northeastern Mediterranean (Iskenderun Bay) population, while the highest was in the Marmara population with overall average value of genetic diversity within populations. A total of 14 haplotypes was found, and the highest haplotype diversity was in the Black Sea whereas the lowest was in the northeastern Mediterranean population (Iskenderun Bay). The Black Sea and Iskenderun Bay populations showed the least genetic divergence (0.001081), while the highest was between the Marmara Sea and northeastern Mediterranean (Antalya Bay) populations (0.002067). Pairwise comparisons of genetic distance revealed statistically significant differences (P < 0.05) between the Marmara and both the Aegean and northeastern Mediterranean (Antalya Bay) samples. Neighbour joining tree analyses clustered the northeastern Mediterranean populations (Antalya Bay and Iskenderun Bay) as genetically more interrelated populations, whereas the Aegean Sea population was clustered as most isolated one. Discriminant function analysis of morphological characters showed that only the Black Sea population is differentiated from the other populations.     

Genetic diversity of natural populations of Larix principis-rupprechtii in Shanxi Province,China

Prince Rupprecht's Larch (Larix principis-rupprechtii Mayr.) is one of dominant components of middle and high elevation forests in North China. Shanxi Province is well known as “the Hometown of Prince Rupprecht's Larch” in China. In this study, six natural populations of this species across Shanxi were selected to investigate the genetic variation of the species using amplified fragment length polymorphism (AFLP) markers. Results showed that in comparison with some other species of Larix, higher genetic diversity was revealed at the species level for L. principis-rupprechtii (percentage of polymorphic loci PPL = 71.9%, Nei's gene diversity H_E = 0.225, Shannon information index I = 0.341). Most of genetic variation existed within populations (80.5%), while the genetic differentiation among populations was significant (p < 0.001) and higher (G_st = 0.194) than most other species of Larix. The differentiation can be attributed to the limited gene flow (N_m = 1.035) among populations, which could be due to the spatial isolation and habitat fragmentation. The six populations can be divided into three groups based on the Nei's genetic distances between populations (from 0.033 to 0.076). There was no significant correlation (r = 0.268, p > 0.05) between genetic distance and geographic distance among populations. The measures for in-situ or ex-situ conservation should be taken to preserve the genetic diversity of this species.     

Genetic structure of expanding wolf (Canis lupus) populations in Italy and Croatia,and the early steps of the recolonization of the Eastern Alps

After centuries of range contraction and demographic declines wolves are now expanding in Europe, colonizing regions from where they have been absent for centuries. Wolf colonizing the western Alps originate by the expansion of the Italian population. Vagrant wolves of Italian and Dinaric-Balkan origins have been recently observed in the Eastern Alps. In this study we compared the genetic structure of wolf populations in Italy and Croatia, aiming to identify the sources of the ongoing recolonization of the Eastern Alps. DNA samples, extracted from 282 Italian and 152 Croatian wolves, were genotyped at 12 autosomal microsatellites (STR), four Y-linked STR and at the hypervariable part of the mitochondrial DNA control-region (mtDNA CR1). Wolves in Croatia and Italy underwent recent demographic bottlenecks, but they differ in genetic diversity and population structure. Wolves in Croatia were more variable at STR loci (N_A = 7.4, H_O = 0.66, H_E = 0.72; n = 152) than wolves in Italy (N_A = 5.3, H_O = 0.57, H_E = 0.58; n = 282). We found four mitochondrial DNA (mtDNA CR1) and 11 Y-STR haplotypes in Croatian wolves, but only one mtDNA CR1 and three Y-STR haplotypes in Italy. Wolves in Croatia were subdivided into three genetically distinct subpopulations (in Dalmatia, Gorski kotar and Lika regions), while Italian wolves were not sub-structured. Assignment testing shows that the eastern and central Alps are recolonized by wolves dispersing from both the Italian and Dinaric populations. The recolonization of the Alps will predictably continue in the future and the new population will be genetically admixed and very variable with greater opportunities for local adaptations and survival.