Glacial bottleneck and postglacial recolonization of a seed parasitic weevil,Curculio hilgendorfi, inferred from mitochondrial DNA variation

Phylogeographical evidence for historical long-distance dispersal in the flightless stick insect Ramulus mikado

Exploring how organisms overcome geographical barriers to dispersal is a fundamental question in biology. Passive long-distance dispersal events, although infrequent and unpredictable, have a considerable impact on species range expansions. Despite limited active dispersal capabilities, many stick insect species have vast geographical ranges, indicating that passive long-distance dispersal is vital for their distribution. A potential mode of passive dispersal in stick insects is via the egg stage within avian digestive tracts, as suggested by experimental evidence. However, detecting such events under natural conditions is challenging due to their rarity. Therefore, to indirectly assess the potential of historical avian-mediated dispersal, we examined the population genetic structure of the flightless stick insect Ramulus mikado across Japan, based on a multifaceted molecular approach [cytochrome oxidase subunit I (COI) haplotypes, nuclear simple sequence repeat markers and genome-wide single nucleotide polymorphisms]. Subsequently, we identified unique phylogeographic patterns, including the discovery of identical COI genotypes spanning considerable distances, which substantiates the notion of passive long-distance genotypic dispersal. Overall, all the molecular data revealed the low and mostly non-significant genetic differentiation among populations, with identical or very similar genotypes across distant populations. We propose that long-distance dispersal facilitated by birds is the plausible explanation for the unique phylogeographic pattern observed in this flightless stick insect.

Evolutionary history and colonization patterns of the wing dimorphic grasshopper Dichroplus vittatus in two Argentinean biomes

Quaternary climate oscillations and modification of the environment by humans have played an important role in shaping species distribution and genetic structure of modern species. Here, population genetic parameters were inferred from the analysis of 168 individuals belonging to 11 populations of the South American grasshopper, Dichroplus vittatus, distributed in two Argentinean Biomes (Grassland and Savanna), by sequencing a 543 bp of the mitochondrial COI gene. Overall, we detected considerable haplotype diversity and low nucleotide diversity. AMOVA analyses showed a significant degree of differentiation among Biomes and between populations. Two major mitochondrial lineages can be distinguished. The haplogroup containing the most common haplotype split 17,000 years BP while the haplogroup including the second most common haplotype has a divergence date of about 11,700 years. Approximate Bayesian Computation (ABC) analyses showed that the palaeodemographic scenario that best fitted our data is consistent with a hypothesis of divergence from an ancestral population and subsequent admixture with Grassland-Savanna (South-North) direction. Our results suggest that populations located in both Biomes would derive from a single ancestral population that colonized the region after the Last Glacial Maximum and Grassland would have a more ancestral origin than Savanna. Further, our results emphasize the importance of human-mediated dispersal in the reconfiguration of genetic diversity of species with potential pest capacity.

Timing Is Everything. Temporal and Spatial Niche Segregation in Curculio spp. (Coleoptera: Curculionidae) Associated with Oak Trees

Oak seed predatory weevils occurring in Poland are prone to increased interspecific competition due to the limited number of Quercus species, compared to southern Europe, in which they can develop. Therefore, analyses on the preferences of three weevil species for acorn sizes chosen for reproduction, as well as on reproductive period duration, were performed. Cafeteria-type experiments were set for females of three species associated with one oak species. Females were allowed to choose and oviposit in acorns of different sizes and growth stages. Research revealed statistically significant differences between the masses of acorns chosen for oviposition by females of Curculio glandium (the biggest), C. pellitus (medium), and C. venosus (the smallest). Studied weevils also differed in terms of the beginning of the reproductive period, which corresponded with the increasing mass of growing acorns. Moreover, C. glandium was observed to be the only species to perform radial egg galleries and lay a considerably higher and varied number of eggs. The results support the hypothesis of a strategy aimed at reducing interspecific competition between Curculio spp. in terms of limited host plant species number.

Phylogeography of the stonefly Kamimuria tibialis: multiple glacial refugia and sympatric occurrence of different lineages in the southern islands of Japan

To elucidate the effect of Pleistocene climatic fluctuations on the historical distribution and geographical genetic structure of temperate Japanese species, we performed phylogeographical and demographic analyses using mitochondrial gene sequences obtained from the stonefly species Kamimuria tibialis, sampled from four main islands of the Japanese Archipelago (i.e. Hokkaido, Honshu, Shikoku and Kyushu) and Tsushima Island. We detected three main clades with distinct geographical distributions, including the Tsushima, Kyushu and Hokkaido–Honshu–Shikoku phylogroups. These groups were estimated to have diverged from one another 0.54–2.02 Mya, suggesting they have undergone several glacial cycles in different refugia. Our results showed that during the glacial epochs and with a fall in sea-level, gene flow was limited among Tsushima and Kyushu, and among Hokkaido and Honshu, probably because the straits between these islands are deep. The population in Kyushu and Shikoku, the southernmost islands, exhibited high genetic diversity, with two distinct haplotype lineages occurring sympatrically. These results suggest that the population division into multiple refugia and the existence of stable southern refugia have contributed to the high genetic diversity of the species in this region.

Inferring historical survivals of climate relicts: the effects of climate changes, geography, and population-specific factors on herbaceous hydrangeas

Climate relicts hold considerable importance because they have resulted from numerous historical changes. However, there are major interspecific variations among the ways by which they survived climate changes. Therefore, investigating the factors and timing that affected population demographics can expand our understanding of how climate relicts responded to historical environmental changes. Here, we examined herbaceous hydrangeas of genus Deinanthe in East Asia, which show limited distributions and a remarkable disjunction between Japan and central China. Chloroplast genome and restriction site-associated DNA sequencing revealed that speciation event occurred in the late Miocene (ca. 7-9 Mya) in response to global climate change. Two lineages apparently remained not branched until the middle Quaternary, and afterwards started to diverge to regional population groups. The narrow endemic species in central China showed lower genetic diversity (He = 0.082), as its population size rapidly decreased during the Holocene due to isolation in montane refugia. Insular populations in the three Japanese islands (He = 0.137-0.160) showed a genetic structure that was inconsistent with sea barriers, indicating that it was shaped in the glacial period when its range retreated to coastal refugia on the exposed sea floor. Demographic modelling by stairway-plot analysis reconstructed variable responses of Japanese populations: some experienced glacial bottlenecks in refugial isolation, while post-glacial range expansion seemingly exerted founder effects on other populations. Overall, this study demonstrated the involvement of not just one, but multiple factors, such as the interplay between climate changes, geography, and other population-specific factors, that determine the demographics of climate relicts.

Nuclear microsatellite and mitochondrial DNA analyses reveal the regional genetic structure and phylogeographical history of a sanguivorous land leech, Haemadipsa japonica, in Japan

Recent molecular studies have indicated that phylogeographical history of Japanese biota is likely shaped by geohistory along with biological events, such as distribution shifts, isolation, and divergence of populations. However, the genetic structure and phylogeographical history of terrestrial Annelida species, including leech species, are poorly understood. Therefore, we aimed to understand the genetic structure and phylogeographical history across the natural range of Haemadipsa japonica, a sanguivorous land leech species endemic to Japan, by using nine polymorphic nuclear microsatellites (nSSR) and cytochrome oxidase subunit one (COI) sequences of mitochondrial DNA (mtDNA). Analyses using nSSR revealed that H. japonica exhibited a stronger regional genetic differentiation among populations (G'ST = 0.77) than other animal species, probably because of the low mobility of land leech. Analyses using mtDNA indicated that H. japonica exhibited two distinct lineages (A and B), which were estimated to have diverged in the middle Pleistocene and probably because of range fragmentation resulting from climatic change and glacial and interglacial cycles. Lineage A was widely distributed across Japan, and lineage B was found in southwestern Japan. Analyses using nSSR revealed that lineage A was roughly divided into two population groups (i.e., northeastern and southwestern Japan); these analyses also revealed a gradual decrease in genetic diversity with increasing latitude in lineage A and a strong genetic drift in populations of northeastern Japan. Combined with the largely unresolved shallow polytomies from the mtDNA phylogeny, these results implied that lineage A may have undergone a rapid northward migration, probably during the Holocene. Then, the regional genetic structure with local unique gene pools may have been formed within each lineage because of the low mobility of this leech species.

Approximate Bayesian computation analysis of EST-associated microsatellites indicates that the broadleaved evergreen tree Castanopsis sieboldii survived the Last Glacial Maximum in multiple refugia in Japan

Climatic changes have played major roles in plants' evolutionary history. Glacial oscillations have been particularly important, but some of their effects on plants' populations are poorly understood, including the numbers and locations of refugia in Asian warm temperate zones. In the present study, we investigated the demographic history of the broadleaved evergreen tree species Castanopsis sieboldii (Fagaceae) during the last glacial period in Japan. We used approximate Bayesian computation (ABC) for model comparison and parameter estimation for the demographic modeling using 27 EST-associated microsatellites. We also performed the species distribution modeling (SDM). The results strongly support a demographic scenario that the Ryukyu Islands and the western parts in the main islands (Kyushu and western Shikoku) were derived from separate refugia and the eastern parts in the main islands and the Japan Sea groups were diverged from the western parts prior to the coldest stage of the Last Glacial Maximum (LGM). Our data indicate that multiple refugia survived at least one in the Ryukyu Islands, and the other three regions of the western and eastern parts and around the Japan Sea of the main islands of Japan during the LGM. The SDM analysis also suggests the potential habitats under LGM climate conditions were mainly located along the Pacific Ocean side of the coastal region. Our ABC-based study helps efforts resolve the demographic history of a dominant species in warm temperate broadleaved forests during and after the last glacial period, which provides a basic model for future phylogeographical studies using this approach.

Modeling the Multiple Facets of Speciation-with-Gene-Flow toward Inferring the Divergence History of Lake Whitefish Species Pairs (Coregonus clupeaformis)

Parallel divergence across replicated species pairs occurring in similar environmental contrasts may arise through distinct evolutionary scenarios. Deciphering whether such parallelism actually reflects repeated parallel divergence driven by divergent selection or a single divergence event with subsequent gene flow needs to be ascertained. Reconstructing historical gene flow is therefore of fundamental interest to understand how demography and selection jointly shaped genomic divergence during speciation. Here, we use an extended modeling framework to explore the multiple facets of speciation-with-gene-flow with demo-genetic divergence models that capture both temporal and genomic variation in effective population size and migration rate. We investigate the divergence history of replicate sympatric species pairs of Lake Whitefish (normal benthic and dwarf limnetic) characterized by variable degrees of ecological divergence and reproductive isolation. Genome-wide SNPs were used to document the extent of genetic differentiation in each species pair, and 26 divergence models were fitted and compared with the unfolded joint allele frequency spectrum of each pair. We found evidence that a recent (circa 3,000–4,000 generations) asymmetrical secondary contact between expanding postglacial populations has accompanied Whitefish diversification. Our results suggest that heterogeneous genomic differentiation has emerged through the combined effects of linked selection generating variable rates of lineage sorting across the genome during geographical isolation, and heterogeneous introgression eroding divergence at different rates across the genome upon secondary contact. This study thus provides a new retrospective insight into the historical demographic and selective processes that shaped a continuum of divergence associated with ecological speciation.

Phylogeography of Nasutitermes corniger (Isoptera: Termitidae) in the Neotropical Region

BACKGROUND

The Neotropical Region is known for its biodiversity and ranks third in number of known termite species. However, biogeographic and phylogeographic information of termites of this region is limited compared to other world geographic regions. Nasutitermes corniger is widely distributed in the region and is of considerable economic importance. The goal of this study was to describe the phylogeography of N. corniger in the Neotropical Region, to better understand its evolutionary processes.

RESULTS

The sampled populations of N. corniger showed high genetic variation. Results indicated strong geographic structure among N. corniger populations, with most haplotypes not broadly shared among separated locations. Phylogeographic analyses showed a dispersal route for N. corniger from Central America into South America via the Isthmus of Panama, with subsequent dispersal through the highlands east of the Andes and into eastern regions of the continent. The majority of haplotypes were limited in distribution to proximal regions, corresponding to particular biomes (Atlantic Forest, Amazonia, Chaco, Cerrado and Caatinga).

CONCLUSIONS

Nasutitermes corniger is suggested to be a good model for biogeographic and phylogeographic studies in the Neotropical Region. This study clarified the phylogeographic history of N. corniger and can contribute to the understanding of biogeographic dispersion processes in the Neotropical Region.

Environmental factors contribute to the formation and maintenance of the contact zone observed in deciduous broad-leaved tree species in Japan

Contact zones are defined as areas where populations from different refugia meet during a postglacial expansion and distinct DNA lineages are mixedly distributed. In Japan, contact zones of various plants and animals were reported from the Kinki-Chugoku region. These contact zones appear to be maintained without any drastic topographic barriers such as those observed in the Alps and Pyrenees Mountains. In this study, the mechanisms underlying the formation and/or maintenance of these contact zones were investigated using six deciduous broad-leaved tree species (Carpinus laxiflora, C. tschonoskii, C. japonica, Magnolia obovata, Padus grayana, and Euonymus oxyphyllus). First, the precise location of the contact zones was examined by intensive genetic analysis of the six species. Second, the relationships between the geographic location of the contact zone and various environmental factors, including climate and topography, were investigated by generalized additive models to reveal the mechanisms of the formation and maintenance of the contact zones. As a result, four of the six examined plant species clearly showed a geographically common contact zone in Hyogo Prefecture and its adjacent areas. The results of the generalized additive models indicate that the pattern of low habitat suitability estimated by ecological niche modeling was the most important factor for determining the location of the common contact zone. These results suggest that areas with low habitat suitability in Hyogo Prefecture restrict the migration and gene flow of the four species in this region, and thus, they maintain the pattern of the contact zones. This study suggests that there are major effects of habitat suitability on the formation and maintenance of the contact zones.

Little gene flow between domestic silkmoth Bombyx mori and its wild relative Bombyx mandarina in Japan, and possible artificial selection on the CAD gene of B. mori

We analyzed PCR-amplified carbamoyl-phosphate synthetase 2, aspartate transcarbamylase, and dihydroorotase (CAD) gene fragments from 146 Bombyx mori native strains and found extremely low levels of DNA polymorphism. Two haplotypes were identified, one of which was predominant. CAD haplotype analysis of 42 samples of Japanese B. mandarina revealed four haplotypes. No common haplotype was shared between the two species and at least five base substitutions were detected. This result was suggestive of low levels of gene flow between the two species. The nucleotide diversity (π) scores of the two samples differed markedly: lower π values were estimated for B. mori native strains than Japanese B. mandarina. We further analyzed 12 Chinese B. mandarina derived from seven areas of China, including Taiwan. The results clearly indicated that the π score was ~80-fold greater in Chinese B. mandarina than in B. mori. The extremely low level of DNA polymorphism in B. mori compared to its wild relatives suggested that the CAD gene itself or its tightly linked regions are possible targets for silkworm domestication.

Loss of genetic variability induced by Agroecosystems: Chrysoperla externa (Hagen) (Neuroptera: Chrysopidae) as a case study

Four species of green lacewings occur in Brazil, of which Chrysoperla externa (Hagen) (Neuroptera: Chrysopidae) exhibits the widest geographical distribution. Chrysoperla externa is a predatory insect that is potentially useful as a biological control agent of agricultural pests. Studies on the genetic diversity of lacewing populations are essential to reduce the environmental and economic harm that may be caused by organisms with a low ability to adapt to the adverse and/or different environmental conditions to which they are exposed. We used the cytochrome oxidase I mitochondrial gene as a molecular marker to investigate the genetic diversity of green lacewing species collected from native and agroecosystem environments. Populations derived from native areas showed higher rates of genetic variability compared to populations from agroecosystems. Demographic changes in the form of population expansion were observed in agroecosystems, whereas populations in the native environment appeared stable over time. A statistical analysis showed significant genetic structure between each of the sampled groups, combined with its complete absence within each group, corroborating each group's identity. We infer that the loss of variability exhibited by populations from the agroecosystems is the result of genetic drift by means of the founder effect, a similar effect that has been observed in other introduced populations. Agroecosystems might therefore function as exotic areas for green lacewings, even when these areas are within the normal range of the species.

Use of the Internal Transcribed Spacer (ITS) Regions to Examine Symbiont Divergence and as a Diagnostic Tool for Sodalis-Related Bacteria

Bacteria excel in most ecological niches, including insect symbioses. A cluster of bacterial symbionts, established within a broad range of insects, share high 16S rRNA similarities with the secondary symbiont of the tsetse fly (Diptera: Glossinidae), Sodalis glossinidius. Although 16S rRNA has proven informative towards characterization of this clade, the gene is insufficient for examining recent divergence due to selective constraints. Here, we assess the application of the internal transcribed spacer (ITS) regions, specifically the ITS(glu) and ITS(ala,ile), used in conjunction with 16S rRNA to enhance the phylogenetic resolution of Sodalis-allied bacteria. The 16S rRNA + ITS regions of Sodalis and allied bacteria demonstrated significant divergence and were robust towards phylogenetic resolution. A monophyletic clade of Sodalis isolates from tsetse species, distinct from other Enterobacteriaceae, was consistently observed suggesting diversification due to host adaptation. In contrast, the phylogenetic distribution of symbionts isolated from hippoboscid flies and various Hemiptera and Coleoptera were intertwined suggesting either horizontal transfer or a recent establishment from an environmental source. Lineage splitting of Sodalis-allied bacteria into symbiotic and free-living sister groups was also observed. Additionally, we propose an ITS region as a diagnostic marker for the identification of additional Sodalis-allied symbionts in the field. These results expand our knowledge of informative genome regions to assess genetic divergence since splitting from the last common ancestor, of this versatile insect symbiont clade that have become increasingly recognized as valuable towards our understanding of the evolution of symbiosis. These facultative and recently associated symbionts may provide a novel source of traits adaptable to the dynamic ecologies encountered by diverse host backgrounds.

Multiple origins of outbreak populations of a native insect pest in an agro-ecosystem

Native insects can become epidemic pests in agro-ecosystems. A population genetics approach was applied to analyze the emergence and spread of outbreak populations of native insect species. Outbreaks of the mirid bug, Stenotus rubrovittatus, have rapidly expanded over Japan within the last two decades. To characterize the outbreak dynamics of this species, the genetic structure of local populations was assessed using polymorphisms of the mtDNA COI gene and six microsatellite loci. Results of the population genetic analysis suggested that S. rubrovittatus populations throughout Japan were genetically isolated by geographic distance and separated into three genetic clusters occupying spatially segregated regions. Phylogeographic analysis indicated that the genetic structure of S. rubrovittatus reflected post-glacial colonization. Early outbreaks of S. rubrovittatus in the 1980s occurred independently of genetically isolated populations. The genetic structure of the populations did not fit the pattern of an outbreak expansion, and therefore the data did not support the hypothesis that extensive outbreaks were caused by the dispersal of specific pestiferous populations. Rather, the historical genetic structure prior to the outbreaks was maintained throughout the increase in abundance of the mirid bug. Our study indicated that changes in the agro-environment induced multiple outbreaks of native pest populations. This implies that, given suitable environmental conditions, local populations may have the potential to outbreak even without invasion of populations from other environmentally degraded areas.

Phylogeography of Phytophagous Weevils and Plant Species in Broadleaved Evergreen Forests: A Congruent Genetic Gap between Western and Eastern Parts of Japan

The Quaternary climate cycles played an important role in shaping the distribution of biodiversity among current populations, even in warm-temperate zones, where land was not covered by ice sheets. We focused on the Castanopsis-type broadleaved evergreen forest community in Japan, which characterizes the biodiversity and endemism of the warm-temperate zone. A comparison of the phylogeographic patterns of three types of phytophagous weevils associated with Castanopsis (a host-specific seed predator, a generalist seed predator, and a host-specific leaf miner) and several other plant species inhabiting the forests revealed largely congruent patterns of genetic differentiation between western and eastern parts of the main islands of Japan. A genetic gap was detected in the Kii Peninsula to Chugoku-Shikoku region, around the Seto Inland Sea. The patterns of western-eastern differentiation suggest past fragmentation of broadleaved evergreen forests into at least two separate refugia consisting of the southern parts of Kyushu to Shikoku and of Kii to Boso Peninsula. Moreover, the congruent phylogeographic patterns observed in Castanopsis and the phytophagous insect species imply that the plant-herbivore relationship has been largely maintained since the last glacial periods. These results reinforce the robustness of the deduced glacial and postglacial histories of Castanopsis-associated organisms.

Climate change and evolutionary adaptations at species' range margins

During recent climate warming, many insect species have shifted their ranges to higher latitudes and altitudes. These expansions mirror those that occurred after the Last Glacial Maximum when species expanded from their ice age refugia. Postglacial range expansions have resulted in clines in genetic diversity across present-day distributions, with a reduction in genetic diversity observed in a wide range of insect taxa as one moves from the historical distribution core to the current range margin. Evolutionary increases in dispersal at expanding range boundaries are commonly observed in virtually all insects that have been studied, suggesting a positive feedback between range expansion and the evolution of traits that accelerate range expansion. The ubiquity of this phenomenon suggests that it is likely to be an important determinant of range changes. A better understanding of the extent and speed of adaptation will be crucial to the responses of biodiversity and ecosystems to climate change.

Historical demography and phylogeography of a specialist bark beetle, Dendroctonus pseudotsugae Hopkins (Curculionidae: Scolytinae)

Contemporary distribution of North American species has been shaped by past glaciation events during the Quaternary period. However, their effects were not as severe in the southern Rocky Mountains and Northern Mexico as elsewhere in North America. In this context, we test hypotheses about the historical demography of Dendroctonus pseudotsugae, based on 136 haplotypes of mitochondrial cytochrome oxidase I. The phylogenetic analysis yielded four haplogroups corresponding to northwestern United States and southwestern Canada (NUS), southwestern United States (Arizona, SUS), northwestern Mexico (Sierra Madre Occidental, SMOC), and northeastern Mexico (Sierra Madre Oriental, SMOR). Predictions of demographic expansion were examined through neutrality tests against population growth and mismatch distribution. Results showed that the NUS and SMOC haplogroups have experienced demographic expansion events, whereas the SUS and SMOR haplogroups have not. Divergence times between pairs of haplogroups were estimated from early to middle Pleistocene. The longer divergence time between NUS and all other haplogroups could be the result of refugia within the Pacific Northwest and northern Rocky Mountains and long-term isolation from southernmost populations in Mexico. The results obtained in this study are in agreement with the evolutionary history of the host Douglas-fir, as the warmer climates of interglacial periods pushed conifers northward of Colorado, New Mexico, and Arizona, whereas environmental changes reduced the population size of Douglas-fir and forced fragmentation of distribution range southward into northern Mexico.

Revisiting the insect mitochondrial molecular clock: the mid-Aegean trench calibration

Phylogenetic trees in insects are frequently dated by applying a "standard" mitochondrial DNA (mtDNA) clock estimated at 2.3% My(-1), but despite its wide use reliable calibration points have been lacking. Here, we used a well-established biogeographic barrier, the mid-Aegean trench separating the western and eastern Aegean archipelago, to estimate substitution rates in tenebrionid beetles. Cytochrome oxidase I (cox1) for six codistributed genera across 28 islands (444 individuals) on both sides of the mid-Aegean trench revealed 60 independently coalescing entities delimited with a mixed Yule-coalescent model. One representative per entity was used for phylogenetic analysis of mitochondrial (cox1, 16S rRNA) and nuclear (Mp20, 28S rRNA) genes. Six nodes marked geographically congruent east-west splits whose separation was largely contemporaneous and likely to reflect the formation of the mid-Aegean trench at 9-12 Mya. Based on these "known" dates, a divergence rate of 3.54% My(-1) for the cox1 gene (2.69% when combined with the 16S rRNA gene) was obtained under the preferred partitioning scheme and substitution model selected using Bayes factors. An extensive survey suggests that discrepancies in mtDNA substitution rates in the entomological literature can be attributed to the use of different substitution models, the use of different mitochondrial gene regions, mixing of intraspecific with interspecific data, and not accounting for variance in coalescent times or postseparation gene flow. Different treatments of these factors in the literature confound estimates of mtDNA substitution rates in opposing directions and obscure lineage-specific differences in rates when comparing data from various sources.

Phylogeography of Nasonia vitripennis (Hymenoptera) indicates a mitochondrial-Wolbachia sweep in North America

Here we report evidence of a mitochondrial-Wolbachia sweep in North American populations of the parasitoid wasp Nasonia vitripennis, a cosmopolitan species and emerging model organism for evolutionary and genetic studies. Analysis of the genetic variation of 89 N. vitripennis specimens from Europe and North America was performed using four types of genetic markers: a portion of the mitochondrial cytochrome oxidase I gene, nine polymorphic nuclear microsatellites, sequences from 11 single-copy nuclear markers and six Wolbachia genes. The results show that the European populations have a sevenfold higher mitochondrial sequence variation than North American populations, but similar levels of microsatellite and nuclear gene sequence variation. Variation in the North American mitochondria is extremely low (pi=0.31%), despite a highly elevated mutation rate (approximately 35-40 times higher than the nuclear genes) in the mitochondria of Nasonia. The data are indicative of a mitochondrial sweep in the North American population, possibly due to Wolbachia infections that are maternally co-inherited with the mitochondria. Owing to similar levels of nuclear variation, the data could not resolve whether N. vitripennis originated in the New or the Old World.

Phylogeographical patterns of a generalist acorn weevil: insight into the biogeographical history of broadleaved deciduous and evergreen forests

Background

Climatic changes during glacial periods have had a major influence on the recent evolutionary history of living organisms, even in temperate forests on islands, where the land was not covered with ice sheets. We investigated the phylogeographical patterns of the weevil Curculio sikkimensis (Curculionidae), a generalist seed predator of Fagaceae plants living in both deciduous oak and evergreen forests of Japan. Its genetic structure was compared to that of another host-specific seed predator, C. hilgendorfi, inhabiting only evergreen forests.

Results

We examined 921 bp of mitochondrial DNA for 115 individuals collected from 33 populations of C. sikkimensis from 11 plant species of three genera, Quercus, Lithocarpus, and Castanopsis. An analysis of molecular variance revealed that a large proportion (almost 50%, P < 0.001) of the total genetic variance could be explained by differences between two geographical regions, the southwestern and northeastern parts of the main islands of Japan. In contrast, no significant genetic differentiation of the weevil was observed among vegetation types of their utilized host plant species. The phylogeographical patterns of the generalist and the host-specific seed predator exhibited a congruent genetic boundary in the Chugoku-Shikoku region.

Conclusion

Our results suggest that geology and historical environment have contributed to shaping the present genetic structure of C. sikkimensis. The geographical patterns of genetic differentiation in the Chugoku-Shikoku region observed in the two types of Fagaceae-associated Curculio in this study have also been observed in several plant species growing in warm and cool temperate zones of Japan. The occurrence of this common pattern suggests that deciduous oak and evergreen forests of Japan survived together, or adjacent to each other, in small refugia during glacial ages, in the southwestern and northeastern parts of the main islands, although these two types of forests are presently distributed in cool and warm temperate zones of Japan, respectively.