Weak large‐scale population genetic structure in a philopatric seabird,the European Shag Phalacrocorax aristotelis

Quantifying population genetic structure is fundamental to testing hypotheses regarding gene flow, population divergence and dynamics across large spatial scales. In species with highly mobile life‐history stages, where it is unclear whether such movements translate into effective dispersal among discrete philopatric breeding populations, this approach can be particularly effective. We used seven nuclear microsatellite loci and mitochondrial DNA (ND2) markers to quantify population genetic structure and variation across 20 populations (447 individuals) of one such species, the European Shag, spanning a large geographical range. Despite high breeding philopatry, rare cross‐sea movements and recognized subspecies, population genetic structure was weak across both microsatellites and mitochondrial markers. Furthermore, although isolation‐by‐distance was detected, microsatellite variation provided no evidence that open sea formed a complete barrier to effective dispersal. These data suggest that occasional long‐distance, cross‐sea movements translate into gene flow across a large spatial scale. Historical factors may also have shaped contemporary genetic structure: cluster analyses of microsatellite data identified three groups, comprising colonies at southern, mid‐ and northern latitudes, and similar structure was observed at mitochondrial loci. Only one private mitochondrial haplotype was found among subspecies, suggesting that this current taxonomic subdivision may not be mirrored by genetic isolation.     

Phylogeography of the pine processionary moth Thaumetopoea wilkinsoni in the Near East

Phylogeographic structure of the eastern pine processionary moth Thaumetopoea wilkinsoni was explored in this study by means of nested clade phylogeographic analyses of COI and COII sequences of mitochondrial DNA and Bayesian estimates of divergence times. Intraspecific relationships were inferred and hypotheses tested to understand historical spread patterns and spatial distribution of genetic variation. Analyses revealed that all T. wilkinsoni sequences were structured in three clades, which were associated with two major biogeographic events, the colonization of the island of Cyprus and the separation of southwestern and southeastern Anatolia during the Pleistocene. Genetic variation in populations of T. wilkinsoni was also investigated using amplified fragment length polymorphisms and four microsatellite loci. Contrasting nuclear with mitochondrial data revealed recurrent gene flow between Cyprus and the mainland, related to the long-distance male dispersal. In addition, a reduction in genetic variability was observed at both mitochondrial and nuclear markers at the expanding boundary of the range, consistent with a recent origin of these populations, founded by few individuals expanding from nearby localities. In contrast, several populations fixed for one single mitochondrial haplotype showed no reduction in nuclear variability, a pattern that can be explained by recurrent male gene flow or selective sweeps at the mitochondrial level. The use of both mitochondrial and nuclear markers was essential in understanding the spread patterns and the population genetic structure of T. wilkinsoni, and is recommended to study colonizing species characterized by sex-biased dispersal.     

Contrasting patterns of mitochondrial and microsatellite genetic structure among Western European populations of tawny owls (Strix aluco)

A recent study of mitochondrial phylogeography of tawny owls (Strix aluco) in western Europe suggested that this species survived the Pleistocene glaciations in three allopatric refugia located in Iberia, Italy, and the Balkans, and the latter was likely the predominant source of postglacial colonization of northern Europe. New data from seven microsatellite loci from 184 individual owls distributed among 14 populations were used to assess the genetic congruence between nuclear and mitochondrial DNA (mtDNA) markers. Microsatellites corroborated the major phylogeographical conclusions reached on the basis of the mtDNA sequences, but also showed important differences leading to novel inferences. Microsatellites corroborated the three major refugia and supported the Balkan origin of northern populations. When corrected for differences in effective population size, microsatellites and mtDNA yielded generally congruent overall estimates of population structure (N*ST=0.12 vs. RST=0.16); however, there was substantial heterogeneity in the RST among the seven nuclear loci that was not correlated with heterozygosity. Populations representing the Balkans postglacial expansion interact with populations from the other two refugia forming two clines near the Alps and the Pyrenees. In both cases, the apparent position of the contact zones differed substantially between markers due to the genetic composition of populations sampled in northern Italy and Madrid. Microsatellite data did not corroborate the lower genetic diversity of northern, recently populated regions as was found with mtDNA; this discrepancy was taken as evidence for a recent bottleneck recovery. Finally, this study suggests that congruence among genetic markers should be more likely in cases of range expansion into new areas than when populations interact across contact zones.     

Molecular and genetic characterization of the allelic variants of Du215, Du281, Du323, and Du47G microsatellite loci in parthenogenetic lizard <Emphasis Type="Italic">Darevskia armeniaca</Emphasis> (Lacertidae)

A key issue in the study of unisexual (parthenogenetic) vertebrate species is the determination of their genetic and clonal diversity. In pursuing this aim, various markers of nuclear and mitochondrial genomes can be used. The most effective genetic markers include microsatellite DNA, characterized by high variability. The development and characterization of such markers is a necessary step in the genetic studies of parthenogenetic species. In the present study, using locus-specific PCR, for the first time, an analysis of allelic polymorphism of four microsatellite loci is performed in the populations of parthenogenetic species Darevskia armeniaca. In the studied populations, allelic variants of each locus are identified, and the nucleotide sequences of each allele are determined. It is demonstrated that allele differences are associated with the variation in the structure of microsatellite clusters and single nucleotide substitutions at fixed distances in flanking DNA regions. Structural allele variations form haplotype markers that are specific to each allele and are inherited from their parental bisexual species. It is established which of the parental alleles of each locus were inherited by the parthenogenetic species. The characteristics of the distribution and frequency of the alleles of microsatellite loci in the populations of D. armeniaca determining specific features of each population are obtained. The observed heterozygosity of the populations at the studied loci and the mutation rates in genome regions, as well as Nei’s genetic distances between the studied populations, are determined, and the phylogenetic relationships between them are established.     

Microgeographic population genetic structure in the northern water snake, Nerodia sipedon sipedon detected using microsatellite DNA loci

We describe the isolation and genetic characterization of eight microsatellite DNA loci from the northern water snake, Nerodia sipedon sipedon and use these loci to analyse levels of genetic differentiation between local (< 2 km apart) populations of these snakes in Ontario. These loci are variable, with expected heterozygosities ranging from 0.28 to 0.91, and can correctly exclude nonsires in parentage analyses with a high probability (0.998). Population analyses reveal significant deviation from expected heterozygosity levels for one population, probably a result of a null allele(s) at a single locus and small but significant levels of genetic differentiation among all three populations. This demonstrates that microgeographic genetic structure exists in this species, possibly due to limited dispersal.     

Subterranean termite phylogeography reveals multiple postglacial colonization events in southwestern Europe

A long‐standing goal of evolutionary biology is to understand how paleoclimatic and geological events shape the geographical distribution and genetic structure within and among species. Using a diverse set of markers (cuticular hydrocarbons, mitochondrial and nuclear gene sequences, microsatellite loci), we studied Reticulitermes grassei and R. banyulensis, two closely related termite species in southwestern Europe. We sought to clarify the current genetic structure of populations that formed following postglacial dispersal from refugia in southern Spain and characterize the gene flow between the two lineages over the last several million years. Each marker type separately provided a fragmented picture of the evolutionary history at different timescales. Chemical analyses of cuticular hydrocarbons and phylogenetic analyses of mitochondrial and nuclear genes showed clear separation between the species, suggesting they diverged following vicariance events in the Late Miocene. However, the presence of intermediate chemical profiles and mtDNA introgression in some Spanish colonies suggests ongoing gene flow. The current genetic structure of Iberian populations is consistent with alternating isolation and dispersal events during Quaternary glacial periods. Analyses of population genetic structure revealed postglacial colonization routes from southern Spain to France, where populations underwent strong genetic bottlenecks after traversing the Pyrenees resulting in parapatric speciation.     

Contemporary genetic structure of Xantus's Hummingbird (Basilinna xantusii) in the Baja California peninsula

Genetic structure and phylogeographic patterns of natural populations are of great importance in assessing the conservation status of species. These population properties can be estimated using molecular markers of either mitochondrial DNA (mtDNA) or nuclear DNA (nDNA) to understand the historical, ecological and dispersal patterns that influence genetic exchange within and between populations. Basilinna xantusii is a sexually dimorphic hummingbird endemic to the Baja California Peninsula (BCP). It comprises three ancestral mitochondrial lineages linked to vicariant events, late Pleistocene climate changes and the geographical distribution of oases. This study aimed to determine and understand the current population genetic structure of this hummingbird. The genotypes of 16 microsatellite loci from 100 individuals collected across the geographical range of this species were compared with mtDNA sequences previously published. Cluster analyses identified five populations, two with almost no genetic admixture in the northern part of the BCP and three others with varying levels of admixed ancestry across the BCP. In San José de Magdalena, at the northernmost end of the range of Xantus's Hummingbird, 40% of individuals collected belong to one genetic cluster and the remaining 60% to another, both genetic clusters showing very little admixed ancestry. We hypothesize that, despite being in sympatry, these individuals do not interbreed, unlike the other populations where individuals showed ancestry coefficients of the other genetic groups. The philopatric behaviour of males and the long-range dispersal capacity of females probably determine the observed genetic differentiation pattern. The mito-nuclear discordance detected could be due to the molecular markers used and to female-biased dispersal. Gene flow is asymmetric in this species, being greater from north to south than vice versa, which is probably related to differences in the seasonality of precipitation across the BCP and to urbanization of the oases.     

Genetic variation of microsatellite loci in a bottlenecked species: the northern hairy-nosed wombat Lasiorhinus krefftii

We investigate the utility of hypervariable microsatellite loci to measure genetic variability remaining in the northern hairy-nosed wombat, one of Australia's rarest mammals. This species suffered a dramatic range and population reduction over the past 120 years and now exists as a single colony of about 70 individuals at Epping Forest National Park, central Queensland. Because our preliminary research on mitochondrial DNA and multilocus DNA fingerprints did not reveal informative variation in this population, we chose to examine variation in microsatellite repeats, a class of loci known to be highly polymorphic in mammals. To assess the suitability of various wombat populations as a reference for comparisons of genetic variability and subdivision we further analysed mitochondrial DNA cytochrome b sequence, using phylogenetic methods. Our results show that appreciable levels of variation still exist in the Epping Forest colony although it has only 41% of the heterozygosity shown in a population of a closely-related species. From museum specimens collected in 1884, we also assessed microsatellite variation in an extinct population of the northern hairy-nosed wombat, from Deniliquin, New South Wales, 2000 km to the south of the extant population. The apparent loss of variation in the Epping Forest colony is consistent with an extremely small effective population size throughout its 120-year decline.     

Comparative allozyme and microsatellite population structure in a narrow endemic plant species, Centaurea corymbosa Pourret (Asteraceae)

Centaurea corymbosa Pourret (Asteraceae) is a narrow endemic species known only from six populations located in a 3-km2 area in the south of France. Earlier field experiments have suggested that pollen and seed dispersal were highly restricted within and among populations. Consistent with the field results, populations were highly differentiated for five allozyme loci and among-population variation fitted an isolation-by-distance model. In the present study, we investigated the genetic structure of C. corymbosa using six microsatellite loci. As with allozymes, microsatellites revealed no within-population structure and a large differentiation among populations. However, allozyme loci were less powerful than microsatellites in detecting the extent of gene flow assessed by assignment tests. The patterns of structuration greatly varied among loci for both types of marker; we suggest that differences in single-locus pattern could mainly be an effect of stochastic variation for allozymes and an effect of variation in mutation rate for microsatellites. In contrast to the multilocus results, the two most polymorphic microsatellite loci did not show any isolation-by-distance pattern. Our results suggest that highly variable loci might not always be the best suited markers to quantify levels of gene flow among populations.     

Isolation and characterization of microsatellite markers from the limpet Cellana grata

This note describes the development of nine polymorphic microsatellite loci in the limpet Cellana grata to investigate population structure and cohort variation in this species. The number of alleles ranged from seven to 22 and observed heterozygosity ranged from 0.62 to 0.95. Deviation from the Hardy-Weinberg equilibrium was detected in two loci, both as a result of heterozygote deficiency. Null alleles were detected in one of these loci. These genetic markers will be used to investigate the genetic structure of C. grata populations, as well as variation among cohorts of this common intertidal species.     

Genetic variation of Avicennia marina (Forsk.) Vierh. (Avicenniaceae) in Vietnam revealed by microsatellite and AFLP markers

Genetic variation of Avicennia marina in the costal area of Vietnam was examined using microsatellite and AFLP markers. By using five microsatellite loci a total of 21 alleles were detected. The average number of alleles per locus per population ranged from 1.667 to 3.000. The observed heterozygosity varied from 0.180 to 0.263, with an average of 0.210 indicating relatively low level of genetic variation comparing to the previous studies on A. marina in the worldwide range. The expected heterozygosity was larger than the observed heterozygosity leading to positive inbreeding coefficients in all the six populations. Highly significant departures from Hardy-Weinberg Equilibrium were detected in four populations. AFLP analysis revealed a total of 386 loci, of which 232 (60.1%) were polymorphic. In congruent with microsatellite markers relatively low levels of genetic variation were detected at both gene and nucleotide levels (H = 0.086; pi = 0.0054). Reduced level of genetic variation was found in the central population, and in the southern populations. Both microsatellite and AFLP markers revealed large genetic differentiation (F(ST) = 0.262 and 0.338, respectively) indicating strong genetic structure among regional populations. Pairwise genetic distance by AFLP showed two populations in the north and the other two in the south are closely related each other.     

Population structure and genetic variation in the endangered Giant Kangaroo Rat (<Emphasis Type="Italic">Dipodomys ingens</Emphasis>)

Populations of the endangered giant kangaroo rat, Dipodomys ingens (Heteromyidae), have suffered increasing fragmentation and isolation over the recent past, and the distribution of this unique rodent has become restricted to 3% of its historical range. Such changes in population structure can significantly affect effective population size and dispersal, and ultimately increase the risk of extinction for endangered species. To assess the fine-scale population structure, gene flow, and genetic diversity of remnant populations of Dipodomys ingens, we examined variation at six microsatellite DNA loci in 95 animals from six populations. Genetic subdivision was significant for both the northern and southern part of the kangaroo rat’s range although there was considerable gene flow among southern populations. While regional gene diversity was relatively high for this endangered species, hierarchical F-statistics of northern populations in Fresno and San Benito counties suggested non-random mating and genetic drift within subpopulations. We conclude that effective dispersal, and therefore genetic distances between populations, is better predicted by ecological conditions and topography of the environment than linear geographic distance between populations. Our results are consistent with and complimentary to previous findings based on mtDNA variation of giant kangaroo rats. We suggest that management plans for this endangered rodent focus on protection of suitable habitat, maintenance of connectivity, and enhancement of effective dispersal between populations either through suitable dispersal corridors or translocations.     

Limited genetic differentiation in Labeo rohita (Hamilton 1822) populations as revealed by microsatellite markers

Labeo rohita, popularly known as rohu is a widely cultured species in the whole Indian subcontinent. Knowledge of the genetic diversity of this species is important to support management and conservation programs which will subsequently help in sustainable production of this species. DNA markers, mostly microsatellite markers are excellent tool to evaluate genetic variation of populations. Genetic variation of three wild and one farm population was assessed using eleven microsatellite loci. In analyzing 192 samples, the number of alleles ranged from 4 to 23; observed heterozygosity 0.500 to 0.870 and expected heterozygosity from 0.389 to 0.878. Exact test for Hardy Weinberg disequilibrium revealed that each riverine sample had at least one locus not in equilibrium except one river. Negative inbreeding coefficients (F_IS) were observed across populations indicating very high level of genetic diversity but little genetic differentiation among populations.     

Phylogeography and conservation of the bull shark (Carcharhinus leucas) inferred from mitochondrial and microsatellite DNA

The bull shark (Carcharhinus leucas) is a widely distributed, large coastal shark species known to travel long distances. These characteristics, coupled with the species?? long life span and late age of maturity, would lead one to predict significant global genetic exchange among bull shark populations. By contrast, data show localized depletion in some areas of large coastal shark fisheries, indicating some geographic isolation may exist. We examined genetic variation in the control region of mitochondrial DNA and at five nuclear microsatellite loci in bull sharks sampled from the western Atlantic to investigate the degree of population subdivision. The average per sample haplotype and nucleotide diversity in the mtDNA (0.51 ± 0.26 and 0.12% ± 0.12, respectively) and expected heterozygosity (0.84) in the microsatellite loci contrast sharply in having lower and higher values (respectively) relative to many other shark species. Significant structure exists between the Brazilian and all northern populations at the mtDNA control region (pairwise ??_ST > 0.8, P < 0.001), but not at the nuclear microsatellite loci. Adjacent northern populations show weak to no genetic differentiation for both markers. These results are congruent with restricted maternal gene flow between populations caused by female site fidelity to nursery areas. We estimate the current effective population size to be around 160,000 and 221,000 individuals for the southern and northern Atlantic populations, respectively. The philopatric habits and the relatively low levels of mtDNA genetic diversity observed in bull sharks must be considered in the conservation of this species. Our results indicate that effective bull shark management strategies will require local, regional, and international attention and cooperation.     

Assessment of genetic diversity in the critically endangered Australian corroboree frogs, Pseudophryne corroboree and Pseudophryne pengilleyi, identifies four evolutionarily significant units for conservation

The iconic and brightly coloured Australian northern corroboree frog, Pseudophryne pengilleyi , and the southern corroboree frog, Pseudophryne corroboree are critically endangered and may be extinct in the wild within 3 years. We have assembled samples that cover the current range of both species and applied hypervariable microsatellite markers and mitochondrial DNA sequences to assess the levels and patterns of genetic variation. The four loci used in the study were highly variable, the total number of alleles observed ranged from 13 to 30 and the average number of alleles per locus was 19. Expected heterozygosity of the four microsatellite loci across all populations was high and varied between 0.830 and 0.935. Bayesian clustering analyses in structure strongly supported four genetically distinct populations, which correspond exactly to the four main allopatric geographical regions in which the frogs are currently found. Individual analyses performed on the separate regions showed that breeding sites within these four regions could not be separated into distinct populations. Twelve mtND2 haplotypes were identified from 66 individuals from throughout the four geographical regions. A statistical parsimony network of mtDNA haplotypes shows two distinct groups, which correspond to the two species of corroboree frog, but with most of the haplotype diversity distributed in P. pengilleyi . These results demonstrate an unexpectedly high level of genetic diversity in both species. Our data have important implications for how the genetic diversity is managed in the future. The four evolutionarily significant units must be protected and maintained in captive breeding programmes for as long as it is possible to do.     

The population genetic structure of a large temperate pollinator species, Bombus pascuorum (Scopoli) (Hymenoptera: Apidae)

The genetic population structure of the bumble bee Bombus pascuorum was studied using six microsatellite loci and a partial sequence of the mitochondrial gene cytochrome b . Eighteen populations from central and northern Europe were included in the analysis. Observed levels of genetic variability and heterozygosity were high. Estimates of population differentiation based on F - and Φ-statistics revealed significant genetic differentiation among B. pascuorum populations and suggest that two partially isolated gene pools, separated by the Alps, do exist. The distribution of mtDNA haplotypes supports this view and presents direct evidence for gene flow across the Alps. Estimates of the number of migrants exchanged among populations north of the Alps suggest that historical events may have left a strong imprint on population structure.     

Reduced genetic variation and strong genetic population structure in the freshwater killifish Valencia letourneuxi (Valenciidae) based on nuclear and mitochondrial markers

The genetic variation of the critically endangered Corfu killifish (Valencia letourneuxi), an endemic freshwater fish species of the western Balkans, was assessed for nine populations sampled in eight water systems in western continental Greece, the Peloponnese and the Ionian Island of Corfu, using mitochondrial and microsatellite markers. The analyses were based on data from three mtDNA regions (D‐loop, COI and 16S rRNA sequences) and 14 microsatellite loci. Samples from the congeneric species Valencia hispanica and the phylogenetically closely related species Aphanius fasciatus were also used in the study as outgroups. Both the mitochondrial and the microsatellite analyses revealed three distinct population groupings associated with the geographical distribution of the populations: one southern group occupying rivers draining to the Patraikos Gulf, the second one including the populations flowing into the Amvrakikos Gulf and the third, more northern group, including the other populations from rivers in Corfu Island and Epirus flowing into the Ionian Sea. Within these groupings there is limited genetic differentiation between populations; in addition, there is reduced intrapopulation genetic variation, evidenced by low heterozygosity values, number of alleles and haplotype diversity. In terms of taxonomic implications and appropriate management actions for conservation, our data suggest that the major population groups should be regarded at least as three distinct conservation units (CUs), with translocation and restocking actions to take place only within the geographical range of the CU concerned. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 334–349.     

Microsatellite DNA markers for Australian geckos

Nine highly polymorphic microsatellite markers were developed for Gehyra variegata from an enriched library screened for tetranucleotide repeats. An average of 19.44 alleles per locus and an average observed heterozygosity of 0.86 were found in 81 samples from fragmented populations in the Western Australian wheatbelt. Different gecko species exhibit varying rates of extinction in fragmented landscapes and the genetic markers were developed to test whether these difference are due to different rates of dispersal. The loci were tested for their ability to cross amplify in 12 species of Australian Gekkonidae and Pygopodidae, in an attempt to identify species which might be amenable to population genetic analysis using these microsatellite markers.     

Development of microsatellite markers in black locust (Robinia pseudoacacia) using a dual‐supression‐PCR technique

Black locust (Robinia pseudoacacia) is an economically and ecologically important tree species in the world. We isolated seven polymorphic microsatellite loci from R. pseudoacacia using a dual‐suppression‐PCR technique. These loci provided microsatellite markers with high polymorphism ranging from three to 12 alleles per locus and expected heterozygosity between 0.538 and 0.944. The markers are now available for more detailed investigation of population genetic structure and pollen and seed dispersal.     

Limited ecological population connectivity suggests low demands on self-recruitment in a tropical inshore marine fish (Eleutheronema tetradactylum: Polynemidae)

The diversity of geographic scales at which marine organisms display genetic variation mirrors the biophysical and ecological complexity of dispersal by pelagic larvae. Yet little is known about the effect of larval ecology on genetic population patterns, partly because detailed data of larval ecology do not yet exist for most taxa. One species for which this data is available is Eleutheronema tetradactylum, a tropical Indo-West Pacific shorefish. Here, we use a partial sequence mitochondrial cytochrome oxidase subunit 1 (COI) marker and five microsatellite loci to survey the genetic structure of E. tetradactylum across northern Australia. Structure was found throughout the range and isolation by distance was strong, explaining approximately 87 and 64% of the genetic variation in microsatellites and mtDNA, respectively. Populations separated by as little as 15 km also showed significant genetic structure, implying that local populations are mainly insular and self-seeding on an ecological time frame. Because the larvae of E. tetradactylum have lower swimming performance and poor orientation compared with other tropical fishes, even modest larval abilities may permit self-recruitment rather than passive dispersal.