High variation and very low differentiation in wide ranging plains zebra (Equus quagga): insights from mtDNA and microsatellites

Patterns of genetic differentiation in the plains zebra ( Equus quagga ) were analysed using mitochondrial DNA control region variation and seven microsatellites. The six morphologically defined subspecies of plains zebra lacked the population genetic structure indicative of distinct evolutionary units. Both marker sets showed high levels of genetic variation and very low levels of differentiation. There was no geographical structuring of mitochondrial DNA haplotypes in the phylogenetic tree, and the plains zebra showed the lowest overall differentiation recorded in any African ungulate studied so far. Arid-adapted African ungulates have shown significant regional genetic structuring in support of the Pleistocene refuge theory. This was not the case in the zebra, and the data are discussed in relation to the impact of Pleistocene climate change on a nonbovid member of the savannah ungulate community. The only other species showing a similar absence of genetic structuring is the African buffalo ( Syncerus caffer ), but this taxon lacks the high levels of morphological variation present in the plains zebra.     

Genetic subdivision, glacial refugia and postglacial recolonization in the golden-striped salamander, Chioglossa lusitanica (Amphibia: urodela)

The golden-striped salamander (Chioglossa lusitanica) is an ecologically specialized species, endemic to north-western Iberia. Patterns of genetic variation were assessed at seven polymorphic enzyme loci and one mitochondrial DNA (mtDNA) marker (cytochrome b) in 17 populations across its range. Estimates of enzyme genetic diversity revealed a high degree of genetic subdivision (FST = 0.68), mainly attributable to the existence of two groups of populations. The groups were located, respectively, north and south of the Mondego River, indicating that this river coincided with a major historical barrier to gene flow. A significant decrease in genetic variability from the Mondego northwards was associated with the Douro and Minho rivers. mtDNA sequence variation revealed a congruent pattern of two haplotype groups (d = 2.2%), with a geographical distribution resembling that of allozymes. The pattern and depth of genetic variation is consistent with the following hypotheses: (i) subdivision of an ancestral range of the species prior to the middle Pleistocene; (ii) secondary contact between populations representing historical refugia; (iii) relatively recent range expansion giving rise to the northern part of the species range; and (iv) loss of genetic variation through founder effects during range expansion across major rivers.     

Population genetic structure of North Atlantic, Mediterranean Sea and Sea of Cortez fin whales, Balaenoptera physalus (Linnaeus 1758): analysis of mitochondrial and nuclear loci

Samples were collected from 407 fin whales, Balaenoptera physalus , at four North Atlantic and one Mediterranean Sea summer feeding area as well as the Sea of Cortez in the Pacific Ocean. For each sample, the sex, the sequence of the first 288 nucleotides of the mitochondrial (mt) control region and the genotype at six microsatellite loci were determined. A significant degree of divergence was detected at all nuclear and mt loci between North Atlantic/Mediterranean Sea and the Sea of Cortez. However, the divergence time estimated from the mt sequences was substantially lower than the time elapsed since the rise of the Panama Isthmus, suggesting occasional gene flow between the North Pacific and North Atlantic ocean after the separation of the two oceans. Within the North Atlantic and Mediterranean Sea, significant levels of heterogeneity were observed in the mtDNA between the Mediterranean Sea, the eastern (Spain) and the western (the Gulf of Maine and the Gulf of St Lawrence) North Atlantic. Samples collected off West Greenland and Iceland could not be unequivocally assigned to either of the two areas. The homogeneity tests performed using the nuclear data revealed significant levels of divergence only between the Mediterranean Sea and the Gulf of St Lawrence or West Greenland. In conclusion, our results suggest the existence of several recently diverged populations in the North Atlantic and Mediterranean Sea, possibly with some limited gene flow between adjacent populations, a population structure which is consistent with earlier population models proposed by Kellogg, Ingebrigtsen, and Sergeant.     

Pronounced genetic population structure in a potentially vagile fish species (Pristipomoides multidens, Teleostei; Perciformes; Lutjanidae) from the East Indies triangle

The East Indies triangle, bordered by the Phillipines, Malay Peninsula and New Guinea, has a high level of tropical marine species biodiversity. Pristipomoides multidens is a large, long-lived, fecund snapper species that is distributed throughout the East Indies and Indo-Pacific. Samples were analysed from central and eastern Indonesia and northern Australia to test for genetic discontinuities in population structure. Fish (n = 377) were collected from the Indonesian islands of Bali, Sumbawa, Flores, West Timor, Tanimbar and Tual along with 131 fish from two northern Australian locations (Arafura and Timor Seas) from a previous study. Genetic variation in the control region of the mitochondrial genome was assayed using restriction fragment length polymorphism and direct sequencing. Haplotype diversity was high (0.67-0.82), as was intraspecific sequence divergence (range 0-5.8%). F(ST) between pairs of populations ranged from 0 to 0.2753. Genetic subdivision was apparent on a small spatial scale; F(ST) was 0.16 over 191 km (Bali/Sumbawa) and 0.17 over 491 km (Bali/Flores). Constraints to dispersal that contribute to, and maintain, the observed degree of genetic subdivision are experienced presumably by all life history stages of this tropical marine finfish. The constraints may include (1) little or no movement of eggs or larvae, (2) little or no home range or migratory movement of adults and (3) loss of larval cohorts due to transport of larvae away from suitable habitat by prevailing currents.     

When east meets west: population structure of a high-latitude resident species,the boreal chickadee (Poecile hudsonicus)

The population genetic structure of northern boreal species has been strongly influenced both by the Quaternary glaciations and the presence of contemporary barriers, such as mountain ranges and rivers. We used a combination of mitochondrial DNA (mtDNA), nuclear microsatellites and spatial distribution modelling to study the population genetic structure of the boreal chickadee (Poecile hudsonicus), a resident passerine, and to investigate whether historical or contemporary barriers have influenced this northern species. MtDNA data showed evidence of eastern and western groups, with secondary admixture occurring in central Canada. This suggests that the boreal chickadee probably persisted in multiple glacial refugia, one in Beringia and at least one in the east. Palaeo-distribution modelling identified suitable habitat in Beringia (Alaska), Atlantic Canada and the southern United States, and correspond to divergence dates of 60–96 kya. Pairwise F_ST values for both mtDNA and microsatellites were significant for all comparisons involving Newfoundland, though mtDNA data suggest a more recent separation. Furthermore, unlike mtDNA data, nuclear data support population connectivity among the continental populations, possibly due to male-biased dispersal. Although both are significant, the isolation-by-distance signal is much stronger for mtDNA (r²=0.51) than for microsatellites (r²=0.05), supporting the hypothesis of male-biased dispersal. The population structure of the boreal chickadee was influenced by isolation in multiple refugia and contemporary barriers. In addition to geographical distance, physical barriers such as the Strait of Belle Isle and northern mountains in Alaska are restricting gene flow, whereas the Rocky Mountains in the west are a porous barrier.     

Phylogeography of the brown hare (Lepus europaeus) in Europe: a legacy of south‐eastern Mediterranean refugia?

Aim We analysed the population genetics of the brown hare (Lepus europaeus) in order to test the hypothesis that this species migrated into central Europe from a number of late glacial refugia, including some in Asia Minor. Location Thirty‐three localities in Greece, Bulgaria, Italy, Croatia, Serbia, Poland, Switzerland, Austria, France, Germany, the Netherlands, Spain, the United Kingdom, Turkey and Israel. Methods In total, 926 brown hares were analysed for mitochondrial DNA (mtDNA) variation by restriction fragment length polymorphism (RFLP) performed on polymerase chain reaction‐amplified products spanning cytochrome b (cyt b)/control region (CR), cytochrome oxidase I (COI) and 12S–16S rRNA. In addition, sequence analysis of the mtDNA CR‐I region was performed on 69 individuals, and the data were compared with 137 mtDNA CR‐I sequences retrieved from GenBank. Results The 112 haplotypes detected were partitioned into five phylogeographically well‐defined major haplogroups, namely the ‘south‐eastern European type haplogroup’ (SEEh), ‘Anatolian/Middle Eastern type haplogroup’ (AMh), ‘European type haplogroup, subgroup A’ (EUh‐A), ‘European type haplogroup, subgroup B’ (EUh‐B) and ‘Intermediate haplogroup’ (INTERh). Sequence data retrieved from GenBank were consistent with the haplogroups determined in this study. In Bulgaria and north‐eastern Greece numerous haplotypes of all five haplogroups were present, forming a large overlap zone. Main conclusions The mtDNA results allow us to infer post‐glacial colonization of large parts of Europe from a late glacial/early Holocene source population in the central or south‐central Balkans. The presence of Anatolian/Middle Eastern haplotypes in the large overlap zone in Bulgaria and north‐eastern Greece reveals gene flow from Anatolia to Europe across the late Pleistocene Bosporus land‐bridge. Although various restocking operations could be partly responsible for the presence of unexpected haplotypes in certain areas, we nevertheless trace a strong phylogeographic signal throughout all regions under study. Throughout Europe, mtDNA results indicate that brown hares are not separated into discernable phyletic groups.     

Beyond the Mediterranean peninsulas: evidence of central European glacial refugia for a temperate forest mammal species, the bank vole (Clethrionomys glareolus)

This study details the phylogeographic pattern of the bank vole, Clethrionomys glareolus, a European rodent species strongly associated with forest habitat. We used sequences of 1011 base pairs of the mitochondrial DNA cytochrome b gene from 207 bank voles collected in 62 localities spread throughout its distribution area. Our results reveal the presence of three Mediterranean (Spanish, Italian and Balkan) and three continental (western, eastern and 'Ural') phylogroups. The endemic Mediterranean phylogroups did not contribute to the post-glacial recolonization of much of the Palaearctic range of species. Instead, the major part of this region was apparently recolonized by bank voles that survived in glacial refugia in central Europe. Moreover, our phylogeographic analyses also reveal differentiated populations of bank voles in the Ural mountains and elsewhere, which carry the mitochondrial DNA of another related vole species, the ruddy vole (Clethrionomys rutilus). In conclusion, this study demonstrates a complex phylogeographic history for a forest species in Europe which is sufficiently adaptable that, facing climate change, survives in relict southern and northern habitats. The high level of genetic diversity characterizing vole populations from parts of central Europe also highlights the importance of such regions as a source of intraspecific genetic biodiversity.     

The genetic population structure of Thunnus thynnus (Linnaeus, 1758) in the Mediterranean Sea,a controversial issue

The Atlantic bluefin tuna (ABFT), Thunnus thynnus (Linnaeus, 1758), is an important commercial species managed as two different stocks, western and eastern Atlantic, with their spawning grounds in the Gulf of Mexico and in the Mediterranean Sea, respectively. The eastern Atlantic stock has been overexploited in the last decades, leading to the application of specific management measures introduced by the International Commission for the Atlantic Tuna (ICCAT). A clear understanding of the genetic structure of ABFT Mediterranean population should be pursued in order to support management decisions. To date the genetic studies on the Mediterranean ABFT, carried out with different molecular markers and sampling procedures, have produced unclear results. Here, we analysed ABFT samples from central and western Mediterranean Sea with mitochondrial sequences and 11 microsatellite loci to investigate, among the others, the area of the Strait of Messina, where environmental conditions seem to support a resident population of ABFT. Furthermore, genetic analyses of mitochondrial sequences were carried out including nucleotide sequences of Adriatic ABFT wild larvae retrieved from GenBank. Among the investigated areas a genetic differentiation was detected between the Strait of Messina and the Tyrrhenian Sea with microsatellite loci according to the exact G test, but not to the Bayesian analyses carried out with STRUCTURE. The analyses with mitochondrial sequences do not reveal any differentiation among sampled areas, however, a highly significant genetic divergence was observed between the Adriatic mitochondrial sequences retrieved from GenBank and the central‐western Mediterranean sequences obtained in the present work. Our results provide some evidence of population structure of Mediterranean ABFT adding pieces to a still unclear picture.     

Strong population structure despite evidence of recent migration in a selfing hermaphroditic vertebrate, the mangrove killifish (Kryptolebias marmoratus)

We employ a battery of 33 polymorphic microsatellite loci to describe geographical population structure of the mangrove killifish (Kryptolebias marmoratus), the only vertebrate species known to have a mixed-mating system of selfing and outcrossing. Significant population genetic structure was detected at spatial scales ranging from tens to hundreds of kilometres in Florida, Belize, and the Bahamas. The wealth of genotypic information, coupled with the highly inbred nature of most killifish lineages due to predominant selfing, also permitted treatments of individual fish as units of analysis. Genetic clustering algorithms, neighbour-joining trees, factorial correspondence, and related methods all earmarked particular killifish specimens as products of recent outcross events that could often be provisionally linked to specific migration events. Although mutation is the ultimate source of genetic diversity in K. marmoratus, our data indicate that interlocality dispersal and outcross-mediated genetic recombination (and probably genetic drift also) play key proximate roles in the local 'clonal' dynamics of this species.     

Genetic and morphological variation in a Mediterranean glacial refugium: evidence from Italian pygmy shrews,Sorex minutus (Mammalia: Soricomorpha)

At the Last Glacial Maximum (LGM), the southern European peninsulas were important refugia for temperate species. Current genetic subdivision of species within these peninsulas may reflect past population subdivision at the LGM, as in ‘refugia within refugia’, and/or at other time periods. In the present study, we assess whether pygmy shrew populations from different regions within Italy are genetically and morphologically distinct. One maternally and two paternally inherited molecular markers (cytochrome b and Y‐chromosome introns, respectively) were analysed using several phylogenetic methods. A geometric morphometric analysis was performed on mandibles to evaluate size and shape variability between populations. Mandible shape was also explored with a functional approach that considered the mandible as a first‐order lever affecting bite force. We found genetically and morphologically distinct European, Italian, and southern Italian groups. Mandible size increased with decreasing latitude and southern Italian pygmy shrews exhibited mandibles with the strongest bite force. It is not clear whether or not the southern Italian and Italian groups of pygmy shrews occupied different refugia within the Italian peninsula at the LGM. It is likely, however, that geographic isolation earlier than the LGM on islands at the site of present‐day Calabria was important in generating the distinctive southern Italian group of pygmy shrews, and also the genetic groups in other small vertebrates that we review here. Calabria is an important hotspot for genetic diversity, and is worthy of conservation attention. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 774–787.     

The declining Spadefoot toad, Pelobates fuscus (Pelobatidae): paleo and recent environmental changes as a major influence on current population structure and status

Populations of the European Spadefoot toad (Pelobates fuscus) have experienced recent declines all over Europe, but these appear to be more intense in north and western Europe. Due to the toad’s fossorial nature and specific habitat requirements, environmental conditions have played a major role in structuring current populations. We examined the phylogeographic structure in P. fuscus from 16 localities throughout Europe using mitochondrial cytochrome b gene sequence analysis. Sequence divergence among haplotypes was low (0.54±0.15%). Three very closely related haplotypes occupy northern and western parts of Europe whereas 12 others were observed among samples from south-eastern Europe, including the Balkans. Our results suggest that toads only recently colonized the northern and western parts of Europe following glacial retreat. This expansion probably took place in steppic-like areas during the younger Dryas cold interval, about 12,900–11,500 years ago. Restricted gene flow with an isolation-by-distance population structure characterises a major part of its distribution range. Based on our results we suggest that the northern and western lineages should be considered as distinct conservation units, while the south-eastern populations from the refugial areas, where nearly all genetic polymorphism occurs and populations appear less vulnerable, should receive special attention.     

Allotetraploid origin and divergence in Eleusine (Chloridoideae, Poaceae): evidence from low-copy nuclear gene phylogenies and a plastid gene chronogram

Background and Aims Eleusine

(Poaceae) is a small genus of the subfamily Chloridoideae exhibiting considerable morphological and ecological diversity in East Africa and the Americas. The interspecific phylogenetic relationships of Eleusine are investigated in order to identify its allotetraploid origin, and a chronogram is estimated to infer temporal relationships between palaeoenvironment changes and divergence of Eleusine in East Africa.

Methods

Two low-copy nuclear (LCN) markers, Pepc4 and EF-1α, were analysed using parsimony, likelihood and Bayesian approaches. A chronogram of Eleusine was inferred from a combined data set of six plastid DNA markers (ndhA intron, ndhF, rps16-trnK, rps16 intron, rps3, and rpl32-trnL) using the Bayesian dating method.

Key Results

The monophyly of Eleusine is strongly supported by sequence data from two LCN markers. In the cpDNA phylogeny, three tetraploid species (E. africana, E. coracana and E. kigeziensis) share a common ancestor with the E. indica–E. tristachya clade, which is considered a source of maternal parents for allotetraploids. Two homoeologous loci are isolated from three tetraploid species in the Pepc4 phylogeny, and the maternal parents receive further support. The A-type EF-1α sequences possess three characters, i.e. a large number of variations of intron 2; clade E-A distantly diverged from clade E-B and other diploid species; and seven deletions in intron 2, implying a possible derivation through a gene duplication event. The crown age of Eleusine and the allotetraploid lineage are 3·89 million years ago (mya) and 1·40 mya, respectively.

Conclusions

The molecular data support independent allotetraploid origins for E. kigeziensis and the E. africana–E. coracana clade. Both events may have involved diploids E. indica and E. tristachya as the maternal parents, but the paternal parents remain unidentified. The habitat-specific hypothesis is proposed to explain the divergence of Eleusine and its allotetraploid lineage.     

High population differentiation and unusual haplotype structure in a shade-intolerant pioneer tree species, Zanthoxylum ailanthoides (Rutaceae) revealed by analysis of DNA polymorphism at four nuclear loci

Differences in demographic history, life-history traits, and breeding systems affect nucleotide variation patterns. It is expected that shade-intolerant pioneer tree species have different patterns of genetic polymorphism and population structure than climax species. We studied patterns of nucleotide polymorphism at four putative starch pathway loci (agpSA, agpSB, agpL, and GBSSI) in Zanthoxylum ailanthoides , a shade-intolerant pioneer tree species that occupies forest gaps in warm-temperate forests of East Asia. Genetic diversity was lower within each population than among populations, and differentiation among populations was high across the loci ( F _ST = 0.32–0.64), as expected from the insect-pollinated breeding system and the metapopulation structure of this pioneer species. Numbers of haplotypes were smaller than those expected from the observed numbers of segregating sites. Single haplotypes accounted for more than 47% of all the sampled genes at the respective loci. These variation patterns were incompatible with neutral predictions for populations of a finite island model. Complex population dynamics, such as bottleneck and/or admixture, in the history of this pioneer tree species might have resulted in the observed patterns of genetic variation and population structure, which are different from those of climax wind-pollinated tree species, such as conifers. In contrast to the other loci investigated in this study, agpL showed nearly no variation in Z. ailanthoides (one singleton only), but there was some extent of variation in a closely related species, Zanthoxylum schinifolium . This suggests possibly a recent selective sweep at or near the locus in Z. ailanthoides .