Genetic relationship of the melon fly,Bactrocera cucurbitae (Diptera: Tephritidae) inferred from mitochondrial DNA

Abstract The melon fruit fly, Bactrocera cucurbitae (Coquillett) (Diptera: Tephritidae), has been the subject of worldwide quarantine and management efforts due to its widespread agricultural impact and potential for rapid range expansion. From its presumed native distribution in India, this species has spread throughout the hot‐humid regions of the world. We provide information that reveals population structure, invasion history and population connectivity from 23 locations covering nine countries based on DNA sequences of the mitochondrial cytochrome oxidase I (COI) gene. Forty‐two polymorphic sites were described among 38 haplotypes. The most common haplotype, H1, was observed in 73% of the samples distributed among all populations. Highest genetic diversity was seen within populations, and no isolation‐by‐distance was detected. The western regions (Nepal, Bangladesh, Thailand, Burma and China‐west) showed higher haplotype diversity than eastern regions (China‐east). China‐Yunnan showed highest levels of genetic diversity in China. Haplotype diversity decreased with longitude from west to east. Together, these analyses suggest that B. cucurbitae has expanded from west to east within a limited geographic scale and recently invaded China through Yunnan Province.     

Population genetic structure and phylogeographical pattern of rice grasshopper, <Emphasis Type="Italic">Oxya hyla intricata</Emphasis>, across Southeast Asia

The rice grasshopper, Oxya hyla intricata, is a rice pest in Southeast Asia. In this study, population genetic diversity and structure of this Oxya species was examined using both DNA sequences and AFLP technology. The samples of 12 populations were collected from four Southeast Asian countries, among which 175 individuals were analysed using mitochondrial DNA cytochrome c oxidase subunit I (COI) sequences, and 232 individuals were examined using amplified fragment length polymorphisms (AFLP) to test whether the phylogeographical pattern and population genetics of this species are related to past geological events and/or climatic oscillations. No obvious trend of genetic diversity was found along a latitude/longitude gradient among different geographical groups. Phylogenetic analysis indicated three deep monophyletic clades that approximately correspond to three geographical regions separated by high mountains and a deep strait, and TCS analysis also revealed three disconnected networks, suggesting that spatial and temporal separations by vicariance, which were also supported by AMOVA as a source of the molecular variance presented among groups. Gene flow analysis showed that there had been frequent historical gene flow among local populations in different regions, but the networks exhibited no shared haplotype among populations. In conclusion, the past geological events and climatic fluctuations are the most important factor on the phylogeographical structure and genetic patterns of O. hyla intricata in Southeast Asia. Habitat, vegetation, and anthropogenic effect may also contribute to gene flow and introgression of this species. Moreover, temperature, abundant rainfall and a diversity of graminaceous species are beneficial for the migration of O. hyla intricata. High haplotype diversity, deep phylogenetic division, negative Fu’s F _s values and unimodal and multimodal distribution shapes all suggest a complicated demographic expansion pattern of these O. hyla intricata populations, which might have been caused by climatic oscillations during glacial periods in the Quaternary.     

Mitochondrial DNA variation in the caramote prawn <Emphasis Type="Italic">Penaeus (Melicertus) kerathurus</Emphasis> across a transition zone in the Mediterranean Sea

In this study we analysed mitochondrial DNA variation in Penaeus kerathurus prawns collected from seven locations along a transect across the Siculo–Tunisian region in order to verify if any population structuring exists over a limited geographical scale and to delineate the putative transition zone with sufficient accuracy. Partial DNA sequences of COI and 16S genes were analysed. In contrast to the highly conservative 16S gene, the COI sequences exhibited sufficient diversity for population analysis. The COI gene revealed low levels of haplotype and nucleotide diversities. The size of the annual landings of this commercial species suggests large population sizes. Hence, the low genetic diversity detected in this study could indicate a possible reduction in effective population sizes in the past. We detected significant genetic differentiation between eastern and western populations likely due to restricted gene flow across the Siculo–Tunisian boundary. We discuss the different evolutionary forces that may have shaped the genetic variation and suggest that the genetic divide is probably maintained by present-day dispersal limitation. R. Zitari-Chatti and N. Chatti are contributed equally to the work.     

Genetic structure of the star sea squirt, Botryllus schlosseri, introduced in southern European harbours

The introduction of new genetic variants or species is often caused by maritime transport between harbours. Botryllus schlosseri is a cosmopolitan ascidian species that is found in both harbours and open shore habitats. In order to determine the influence of ship traffic on the genetic structure and phylogeography of B. schlosseri in southern Europe, we analyzed the variability of a fragment of the mitochondrial gene cytochrome c oxidase subunit I (COI). We sampled seven Atlanto-Mediterranean harbour populations and three open-shore populations. In addition, we sequenced some colonies from the US-Atlantic coast and from other Mediterranean localities to perform phylogenetic analyses. Although the number of polymorphic sites recorded (25.8%) was within the range observed in other population studies based on ascidian COI sequences, the haplotypic diversity (16 haplotypes out of 181 sequences) was much lower. Moreover, a lack of intermediate haplotypes was observed. This pattern of high nucleotide diversity and low haplotype diversity was consistent with introduction events of a few divergent haplotypes. We found a strong genetic structure in the study populations. Gene flow was only appreciable between some harbour populations. Harbour- and open-shore populations were well differentiated, although there was no evidence for isolation by distance. A nested clade analysis pointed to long-distance colonization, possibly coupled with subsequent fragmentation, as the underlying process. Our results suggest that B. schlosseri entered the study area via harbour-hopping, possibly through recurrent introduction events. The haplotypes from North America and most of the European ones were grouped in the same phylogenetic clade. This suggests occasional gene flow between both continents, probably through ship transport.     

Genetic Diversity of the root‐knot nematode Meloidogyne enterolobii in Mulberry Based on the Mitochondrial COI Gene

This study explores the genetic diversity and structure of Meloidogyne enterolobii in mulberry in China. The COI mitochondrial gene (mtCOI) in M.enterolobii populations in Guangdong, Guangxi, and Hunan Provinces was PCR‐amplified, sequenced, and analyzed for genetic diversity. The total number of variations, haplotypes (Hap), the average number of nucleotide differences (k), haplotype diversity (H), and nucleotide diversity (π) of mtCOI were 25, 11, 4.248, 0.900, and 0.00596, respectively. Insignificant differences in Fst value (0.0169) and a high level of gene flow (7.02) were detected among the 19‐mulberry root‐knot nematode populations, and high genetic variation within each population and a small genetic distance among populations were observed. Both phylogenetic analyses and network mapping of the 11 haplotypes revealed a dispersed distribution pattern of 19 mulberry root‐knot nematode populations and an absence of branches strictly corresponding to the 19 range sampling sites. The neutrality test and mismatch analysis indicated that mulberry root‐knot nematode populations experienced a population expansion in the past. The analysis of molecular variance (AMOVA) revealed that the genetic differentiation of M. enterolobii was mainly contributed by the variation within each group. No significant correlation was found between the genetic distance and geographical distance of M. enterolobii populations. The findings of this study provide a profound understanding of the M. enterolobii population and will inform the development of strategies to combat and manage root‐knot nematodes in mulberry.     

Phylogeography of a mite,Halozetes fulvus,reflects the landscape history of a young volcanic island in the sub‐Antarctic

Previous studies of the microarthropods of Marion Island, Southern Ocean, documented high mitochondrial COI (cytochrome c oxidase subunit I) haplotype diversity and significant genetic structure, which were ascribed to landscape subdivision. In this paper we revisit these ideas in light of new geomorphological evidence indicating a major lineament orientated along N26.5°E. Using the microarthropod Halozetes fulvus, we test the hypothesis that the eastern and western sides of the island show different population genetic patterns, corresponding to the previously unrecognized geological separation of these regions, and perhaps also with differences in climates across the island and further landscape complexity. Mitochondrial COI data were collected for 291 H. fulvus individuals from 30 localities across the island. Notwithstanding our sampling effort, haplotype diversity was under‐sampled as indicated by rarefaction analyses. Overall, significant genetic structure was found across the island as indicated by Φ_ST analyses. Nested clade phylogeographical analyses suggested that restricted gene flow (with isolation‐by‐distance) played a role in shaping current genetic patterns, as confirmed by Mantel tests. At the local scale, coalescent modelling revealed two different genetic patterns. The first, characterizing populations on the south‐western corner of the island, was that of low effective population size and high gene flow. The converse was found on the eastern side of Marion Island. Taken together, substantial differences in spatial genetic structure characterize H. fulvus populations across Marion Island, in keeping with the hypothesis that the complex history of the island, including the N26.5°E geological lineament, has influenced population genetic structure. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 105 , 131–145.     

A phylogeographical analysis across three biogeographical provinces of the south-eastern Pacific: the case of the marine gastropod Concholepas concholepas

Aim The present‐day population structure of a species reflects the influence of population history as well as contemporary processes. Little is known about the mechanisms that have shaped the geographical distribution of genetic diversity in marine species present on the south‐eastern Pacific (SEP) coast. Here we provide the first comprehensive phylogeographical study of a species distributed along the SEP coast by analysing the endemic and emblematic muricid gastropod Concholepas concholepas. Location The study localities were distributed along the SEP coast ranging from Matarani (11° S) to Puerto Eden (49° S), crossing three major biogeographical provinces: Peruvian Province, Intermediate Area and Magellanic Province. Methods A total of 337 specimens of C. concholepas were collected from 14 localities in the three biogeographical provinces/areas. Mitochondrial cytochrome oxidase I (COI) gene partial sequences (658 bp) were obtained and analysed using coalescence‐based methods to infer molecular diversity and phylogeographical patterns. Results Across the 337 individuals, we found a large diversity, with a total of 179 haplotypes at the COI gene fragment. Although a slight decrease in gene diversity was observed from north to south, an analysis of molecular variance did not reveal any significant spatial population differentiation from Peru to the tip of Chile, not even across the recognized biogeographical boundaries at 30° S and 42° S. In addition, a star‐like haplotype network suggested the past occurrence of a rapid demographic and geographical expansion over the total range examined. Calculations of the onset of this expansion suggest that it might be due to climatic conditions during the period of the marine isotope stage 11 (MIS 11, 400,000 years ago), the longer and warmer interglacial episode during the Pleistocene epoch. Main conclusions Our phylogeographical analyses indicate that in the recent past C. concholepas mitochondrial DNA lineages underwent a sudden population expansion event. In addition, our data do not support the hypothesis of concordance between biogeographical barriers and phylogeographical breaks along the SEP coast. These two results are in accordance with the paradigm of high larval dispersal ability in marine species with an extended pelagic larval phase.     

Genetic diversity and population connectivity of the Asian green mussel Perna viridis in South China Sea,inferred from mitochondria DNA markers

The Asian green mussel Perna viridis is ecologically and economically important in the coastal regions of China. In order to characterize the genetic diversity and population connectivity of P. viridis in South China Sea, a 664 bp region of mitochondrial COI gene and a 293 bp region of 16S rRNA gene were sequenced and analyzed for 78 and 92 individuals from four populations in South China Sea, respectively. A total of 15 haplotypes were defined by 14 variable nucleotide sites in COI gene, and 7 haplotypes by 6 variable nucleotide sites in 16S rRNA gene. High haplotype diversity and low nucleotide diversity were observed in COI gene, while moderate haplotype diversity and low nucleotide diversity were observed in 16S rRNA gene. Pairwise F_ST values of COI gene were all negative and those of 16S rRNA gene ranged from −0.01409 to 0.10289. The results showed that no significant genetic divergence (or shallow genetic structure) and high levels of population connectivity among the four populations of P. viridis in South China Sea.     

Phylogeography of Hyalomma (Euhyalomma) lusitanicum (Acarina,Parasitiformes, Ixodidae) in Andalusia based on mitochondrial cytochrome oxidase I gene

The genetic population structure relationships of Hyalomma (Euhyalomma) lusitanicum in Andalusia (the south of the Iberian Peninsula) were examined using mtDNA sequence data from 887 bp of cytochrome oxidase subunit I (COI) gene. The sequence for the COI region was determined for 84 individuals collected in several localities of Andalusia, and 10 for other localities (i.e., five from Toledo, central Iberian Peninsula, four from Sicily (Italy) and one from Canary Island). Seventeen haplotypes were detected, including 27 polymorphic sites. The number of amino acid substitutions per site from mean diversity calculations for the entire population was 0.017. AMOVA analysis revealed a low gene flow that characterises the genetic population structure of this species in South Iberian Peninsula, with a haplotype diversity (h) value of 0.815. No geographically induced differentiation was observed, and separate evolutionary units were not detected. Our results indicate low genetic diversity across the geographical range of H. lusitanicum tick in Andalusia. Our data do not show any genetic discontinuity between the tick populations studied, including specimens from Canary Island and Sicily (Italy).

     

Genetic variation and population structure of the carpet shell clam <Emphasis Type="Italic">Ruditapes decussatus</Emphasis> along the Tunisian coast inferred from mtDNA and ITS1 sequence analysis

Surveys of allozyme polymorphisms in the carpet shell clam Ruditapes decussatus have revealed sharp genetic differentiation of populations. Analysis of population structure in this species has now been extended to include nuclear and mitochondrial genes. A partial sequence of a mitochondrial COI gene and of the internal transcribed spacer region (ITS-1) were used to study haplotype distribution, the pattern of gene flow, and population genetic structure of R. decussatus. The samples were collected from twelve populations from the eastern and western Mediterranean coasts of Tunisia, one from Concarneau and one from Thau. A total of twenty and twenty-one haplotypes were detected in the examined COI and ITS1 regions respectively. The study revealed higher levels of genetic diversity for ITS1 compared to COI. The analysis of haplotype frequency distribution and molecular variation indicated that the majority of the genetic variation was distributed within populations (93% and 86% for COI and ITS1 respectively). No significant differentiation was found among eastern and western groups on either side of the Siculo-Tunisian strait. However, distinct and significant clinal changes in haplotypes frequencies between eastern and western samples were found at the most frequent COI haplotype and at three out of five major ITS1 haplotypes. These results suggest the relative importance of historical processes and contemporary hydrodynamic features on the observed patterns of genetic structure.     

BIODIVERSITY RESEARCH: Genetic diversity in two introduced biofouling amphipods (Ampithoe valida & Jassa marmorata) along the Pacific North American coast: investigation into molecular identification and cryptic diversity

Aim We investigated patterns of genetic diversity among invasive populations of Ampithoe valida and Jassa marmorata from the Pacific North American coast to assess the accuracy of morphological identification and determine whether or not cryptic diversity and multiple introductions contribute to the contemporary distribution of these species in the region. Location Native range: Atlantic North American coast; Invaded range: Pacific North American coast. Methods We assessed indices of genetic diversity based on DNA sequence data from the mitochondrial cytochrome c oxidase subunit I (COI) gene, determined the distribution of COI haplotypes among populations in both the invasive and putative native ranges of A. valida and J. marmorata and reconstructed phylogenetic relationships among COI haplotypes using both maximum parsimony and Bayesian approaches. Results Phylogenetic inference indicates that inaccurate species‐level identifications by morphological criteria are common among Jassa specimens. In addition, our data reveal the presence of three well supported but previously unrecognized clades of A. valida among specimens in the north‐eastern Pacific. Different species of Jassa and different genetic lineages of Ampithoe exhibit striking disparity in geographic distribution across the region as well as substantial differences in genetic diversity indices. Main conclusions Molecular genetic methods greatly improve the accuracy and resolution of identifications for invasive benthic marine amphipods at the species level and below. Our data suggest that multiple cryptic introductions of Ampithoe have occurred in the north‐eastern Pacific and highlight uncertainty regarding the origin and invasion histories of both Jassa and Ampithoe species. Additional morphological and genetic analyses are necessary to clarify the taxonomy and native biogeography of both amphipod genera.     

MtDNA COI‐COII marker and drone congregation area: An efficient method to establish and monitor honeybee (Apis mellifera L.) conservation centres

Honeybee subspecies have been affected by human activities in Europe over the past few decades. One such example is the importation of nonlocal subspecies of bees which has had an adverse impact on the geographical repartition and subsequently on the genetic diversity of the black honeybee Apis mellifera mellifera. To restore the original diversity of this local honeybee subspecies, different conservation centres were set up in Europe. In this study, we established a black honeybee conservation centre Conservatoire de l'Abeille Noire d'Ile de France (CANIF) in the region of Ile‐de‐France, France. CANIF's honeybee colonies were intensively studied over a 3‐year period. This study included a drone congregation area (DCA) located in the conservation centre. MtDNA COI‐COII marker was used to evaluate the genetic diversity of CANIF's honeybee populations and the drones found and collected from the DCA. The same marker (mtDNA) was used to estimate the interactions and the haplotype frequency between CANIF's honeybee populations and 10 surrounding honeybee apiaries located outside of the CANIF. Our results indicate that the colonies of the conservation centre and the drones of the DCA show similar stable profiles compared to the surrounding populations with lower level of introgression. The mtDNA marker used on both DCA and colonies of the conservation centre seems to be an efficient approach to monitor and maintain the genetic diversity of the protected honeybee populations.     

Evidence of complex phylogeographic structure for the threatened rodent <Emphasis Type="Italic">Leopoldamys neilli</Emphasis>, in Southeast Asia

Leopoldamys neilli is a threatened murine rodent species endemic to limestone karsts of Thailand. We have studied the phylogeography of L. neilli using two mitochondrial markers (cytb, COI) and one nuclear fragment (bfibr), in order to assess the influence of its endemicity to karst habitat. One hundred fifteen individuals of L. neilli were collected in 20 localities throughout the geographic range of this species in Thailand. Our study revealed strong geographic structure of the mtDNA genetic diversity: six highly differentiated, allopatric genetic lineages were observed in our dataset. They exhibit a very high degree of genetic divergence, low gene flow among lineages and low levels of haplotype and nucleotide diversities within lineages. Our results suggest that L. neilli’s populations are highly fragmented due to the scattered distribution of its karst habitat. The most divergent lineage includes the populations from western Thailand, which have been separated from the other genetic lineages since at least the Early Pleistocene. The other lineages are more closely related and have diverged since the Middle Pleistocene. This study revealed an unexpected high level of genetic differentiation within L. neilli and highlighted the high endemicity of this species to limestone karsts. Our results enhance the importance of protecting limestone habitats to preserve not only the species but also intraspecific diversity.     

Low genetic diversity and lack of genetic structure in the giant jellyfish Nemopilema nomurai in Chinese coastal waters

The giant jellyfish Nemopilema nomurai is a scyphozoan species well-known in East Asian Marginal Seas for its damage to fisheries. The genetic diversity and population structure of N. nomurai, collected from five geographic regions in Chinese coastal seas, were examined based on mitochondrial cytochrome c oxidase subunit I (COI) gene and nuclear internal transcribed spacer one (ITS1) sequences. A total of 26 and five unique haplotypes were recovered from the COI and ITS1 genes, respectively. The overall genetic diversity of N. nomurai calculated by the COI and ITS1 sequences was low (haplotype diversity 0.727% and 0.108%, nucleotide diversity 0.212% and 0.039%). The median-joining network analysis revealed a star-like haplotype network. The hierarchical Analysis of Molecular Variance (AMOVA) of COI haplotypes showed that N. nomurai populations form a single population, with a low F_ST (0.0149, p?=?0.1036). The dispersal ability, together with the biological characteristics, could be important factors for the lack of a geographically structured pattern in N. nomurai in Chinese coastal waters.     

Long-term changes in population genetic features of a rapidly expanding marine invader: implication for invasion success

Large blooms of Rhopilema nomadica, a highly venomous rhizostamatid scyphozoan species introduced to the Mediterranean through the Suez Canal, have become ubiquitous in the summer and winter months along the Israeli coasts since the mid-1980s. This species has since spread across the eastern Mediterranean and was sighted as far west as Tunisia and Sardinia. For the past 12 years, we have studied changes in the mitochondrial COI haplotypes diversity of R. nomadica to investigate small scale fluctuactions of genetic diversity and to reveal possible genetic structuring of the fast spreading invader in the Eastern Mediterranean. The 1091 COI sequences analysed, revealed a highly diverse population displaying 89 haplotypes, 46 of which appeared as singletons, low frequency haplotypes. All the specimens analysed throughout the period belong to a single unstructured population. Though lacking data from the source population in the Red Sea, the high within-population diversity and the high diversity of COI haplotypes support the hypothesis of multiple introductions events, or an open corridor with a continuous influx of propagules. Tajima’s D and Fu’s Fst negative values and the increased numbers of COI singletons from early to late sampling periods, have verified that the Israeli population is characterized by a rapid expanding population. Further research is needed for the evaluation of COI diversity and patterns in R. nomadica populations across the eastern Mediterranean Sea and Red Sea, as well as any correlation of the high variability between COI locus and phenotypic diversity.

     

Variation in mitochondrial cytochrome c oxidase subunit I gene in Nezara viridula (Hemiptera: Pentatomidae) from Argentina

Here, we examine the genetic diversity in the agricultural pest Nezara viridula (Linnaeus, 1758 ) from populations of Argentina using mitochondrial cytochrome c oxidase subunit I (COI) gene sequences. The DNA sequence comparisons of 718 base pairs of the COI gene revealed seven haplotypes. The observed total haplotype diversity (Hd) value was of 0.138, and the nucleotide diversity was of 0.00039 and 0.00135 according to π and θ_W, respectively. Eight out of the 10 populations analysed, mostly from soya bean crops, only presented the more frequent haplotype, while 2 haplotypes were found in a mixed culture and 6 haplotypes in a peanut culture. Factors such as differential insecticide applications, as well as the surrounding habitat, and the host plant preference could be related to the genetic diversity differences observed among samples of N. viridula. The analysis of genetic diversity in samples collected in crops treated and non‐treated with insecticides, as well as in samples collected from different seasons, could help to clarify the role of the factors that led to the pattern of genetic diversity detected in this study. The result of a comparative analysis of COI gene sequences among populations from South America, Africa, Asia and Europe was consistent with the hypothesis of an African origin of N. viridula. On the other hand, the haplotypes of Europe were clustered with haplotypes from South America. In addition, specimens from Madeira (west of Europe) shared ancestry with South America and Europe. It has been suggested that a probable route of colonization of America could have been from Western Europe towards the eastern coasts of South America.     

Phylogeography of the European spiny lobster (Palinurus elephas): Influence of current oceanographical features and historical processes

The European spiny lobster (Palinurus elephas) is a suitable model organism to study the effects of past history and current oceanographic processes on the genetic diversity and population structure of marine species with a long-lived larval phase. A portion of the COI gene was sequenced in 227 individuals from 11 localities, covering most of the present distribution of the species. Divergence was found between Atlantic and Mediterranean regions, which could be explained by restricted gene flow between populations. Moreover, a principal component analysis detected differences within basins. The existence of genetic differentiation between Brittany and Ireland–Scotland populations could be accounted for by the large effect of the Gulf Stream, while mesoscale processes suffered by the incoming Atlantic waters could be responsible of genetic differentiation within the Mediterranean. Furthermore, historical processes could be responsible for a reduction on the overall genetic variability of P. elephas. The haplotypic distribution found in P. elephas, with the presence of one abundant haplotype and a large number of closely related haplotypes, is typical of species experiencing reduction in variability and subsequent expansions. Climatic fluctuations related to glacial cycles could explain the present level of variability and nucleotide diversity found. Interestingly, these glacial events do not seem to have the same impact in other species of the same genus. Our results indicate that recent glacial events could have had a lower impact on Palinurus mauritanicus, a congeneric species that presents an overlapping distribution area but is found in cooler waters than P. elephas.     

Contrasting phylogeographical patterns of two closely related species, Machilus thunbergii and Machilus kusanoi (Lauraceae), in Taiwan

Aim The purpose of this paper was to study the patterns of genetic variation, demographic history, haplotype relationships and potential location of diversity centres of two closely related species, Machilus thunbergii and Machilus kusanoi. Location The phylogeography of M. thunbergii and M. kusanoi was examined by sampling 110 and 106 individuals from 25 and 16 sampling sites, respectively, across their distributional range in Taiwan. Machilus thunbergii is distributed on the Asian mainland, South Korea, southern Japan, the Ryukyus, Taiwan and the Philippines, whereas M. kusanoi is endemic to Taiwan. These two species are closely related, and both are widely distributed in Taiwan but occupy different altitudinal zones and habitats. Methods The range‐wide variation of M. thunbergii and M. kusanoi in Taiwan was studied using chloroplast DNA (cpDNA) variations. A haplotype network was constructed with the computer program tcs . Nested clade analysis was conducted with the computer program ceodis , and various parameters of genetic diversity were calculated and neutrality tested by the computer program Dna SP. Population differentiation was estimated using the programs arlequin and hapstep . The contribution of the populations to gene diversity and to allelic richness was calculated using the software contrib . The level of divergence for each population from the remaining populations was calculated as the mean value of pairwise F_ST for each population against the rest of the populations. Results Extremely low levels of genetic differentiation were found for both species. This result suggested that these two species probably survived in multiple relict refugia with different population sizes throughout the island during low‐temperature periods of the Pleistocene. In addition, nested clade analysis (NCA) of cpDNA haplotypes indicated that restricted gene flow with isolation‐by‐distance characterized the recolonization after the Pleistocene by Tashueshan and Shiouhluan populations of M. thunbergii in the north‐central area west of the Central Mountain Range (CMR). In contrast, NCA analysis indicated that a major diversity centre on the southern tip of the island (Kending population) and contiguous range expansion characterized the recolonization by M. kusanoi of northern areas along the east side of the CMR. The major diversity centres found for the two species examined were further supported by the results of the mean F_ST for individual populations in comparison with other populations, and of the contribution of the divergence component to the total diversity. Main conclusions This research supports the multiple relict refugia hypothesis for both species investigated. Populations of M. thunbergii at Shiouhluan and Tashueshan in the north‐central area west of the CMR represent a diversity centre currently expanding its size. A diversity centre at the southern‐edge population of M. kusanoi, and a contiguous range expansion from Kending, were found. These results indicate that the M. thunbergii populations at Tashueshan and Shiouhluan and the M. kusanoi population at Kending, and even Soukar, are evolutionarily significant units for conservation programmes.     

The legacy of ice ages in mountain species: post‐glacial colonization of mountain tops rather than current range fragmentation determines mitochondrial genetic diversity in an endemic Pyrenean rock lizard

Aim The genetic impact of Quaternary climatic fluctuations on mountain endemic species has rarely been investigated. The Pyrenean rock lizard (Iberolacerta bonnali) is restricted to alpine habitats in the Pyrenees where it exhibits a highly fragmented distribution between massifs and between habitats within massifs. Using mitochondrial DNA markers, we set out: (1) to test whether several evolutionary units exist within the species; (2) to understand how the species persisted through the Last Glacial Maximum and whether the current range fragmentation originates from distribution shifts after the Last Glacial Maximum or from more ancient events; and (3) to investigate whether current mitochondrial diversity reflects past population history or current habitat fragmentation. Location The Pyrenees in south‐western France and northern Spain. Methods We used variation in the hypervariable left domain of the mitochondrial control region of 146 lizards collected in 15 localities, supplemented by cytochrome b sequences downloaded from GenBank to cover most of the species’ distribution range. Measures of population genetic diversity were contrasted with population isolation inferred from topography. Classical (F‐statistics) and coalescence‐based methods were used to assess the level of gene flow and estimate divergence time between populations. We used coalescence‐based simulations to test the congruence of our genetic data with a scenario of simultaneous divergence of current populations. Results Coalescence‐based analyses suggested that these peripheral populations diverged simultaneously at the end of the last glacial episode when their habitats became isolated on mountain summits. High mitochondrial diversity was found in peripheral, isolated populations, while the populations from the core of the species’ range were genetically impoverished. Where mitochondrial diversity has been retained, populations within the same massif exhibited high levels of genetic differentiation. Main conclusions As suggested for many other mountain species, the Pyrenean rock lizard survived glacial maxima through short‐distance range shifts instead of migration or contraction in distant southern refugia. Most of the main Pyrenean range has apparently been re‐colonized from a single or a few source populations, resulting in a loss of genetic diversity in re‐colonized areas. As a result, current levels of intra‐population mitochondrial diversity are better explained by post‐glacial population history than by current habitat fragmentation. Genetic population differentiation within massifs implies severe reduction in female‐mediated gene flow between patches of habitats.     

Phylogeography of the world's tallest angiosperm, Eucalyptus regnans: evidence for multiple isolated Quaternary refugia

Aim There is a need for more Southern Hemisphere phylogeography studies, particularly in Australia, where, unlike much of Europe and North America, ice sheet cover was not extensive during the Last Glacial Maximum (LGM). This study examines the phylogeography of the south‐east Australian montane tree species Eucalyptus regnans. The work aimed to identify any major evolutionary divergences or disjunctions across the species’ range and to examine genetic signatures of past range contraction and expansion events. Location South‐eastern mainland Australia and the large island of Tasmania. Methods We determined the chloroplast DNA haplotypes of 410 E. regnans individuals (41 locations) based on five chloroplast microsatellites. Genetic structure was examined using analysis of molecular variance (AMOVA), and a statistical parsimony tree was constructed showing the number of nucleotide differences between haplotypes. Geographic structure in population genetic diversity was examined with the calculation of diversity parameters for the mainland and Tasmania, and for 10 regions. Regional analysis was conducted to test hypotheses that some areas within the species’ current distribution were refugia during the LGM and that other areas have been recolonized by E. regnans since the LGM. Results Among the 410 E. regnans individuals analysed, 31 haplotypes were identified. The statistical parsimony tree shows that haplotypes divided into two distinct groups corresponding to mainland Australia and Tasmania. The distribution of haplotypes across the range of E. regnans shows strong geographic patterns, with many populations and even certain regions in which a particular haplotype is fixed. Many locations had unique haplotypes, particularly those in East Gippsland in south‐eastern mainland Australia, north‐eastern Tasmania and south‐eastern Tasmania. Higher haplotype diversity was found in putative refugia, and lower haplotype diversity in areas likely to have been recolonized since the LGM. Main conclusions The data are consistent with the long‐term persistence of E. regnans in many regions and the recent recolonization of other regions, such as the Central Highlands of south‐eastern mainland Australia. This suggests that, in spite of the narrow ecological tolerances of the species and the harsh environmental conditions during the LGM, E. regnans was able to persist locally or contracted to many near‐coastal refugia, maintaining a diverse genetic structure.