Nested clade analyses of phylogeographic data: testing hypotheses about gene flow and population history

Since the 1920s, population geneticists have had measures that describe how genetic variation is distributed spatially within a species' geographical range. Modern genetic survey techniques frequently yield information on the evolutionary relationships among the alleles or haplotypes as well as information on allele frequencies and their spatial distributions. This evolutionary information is often expressed in the form of an estimated haplotype or allele tree. Traditional statistics of population structure, such as F statistics, do not make use of evolutionary genealogical information, so it is necessary to develop new statistical estimators and tests that explicitly incorporate information from the haplotype tree. One such technique is to use the haplotype tree to define a nested series of branches (clades), thereby allowing an evolutionary nested analysis of the spatial distribution of genetic variation. Such a nested analysis can be performed regarding the geographical sampling locations either as categorical or continuous variables (i.e. some measure of spatial distance). It is shown that such nested phylogeographical analyses have more power to detect geographical associations than traditional, nonhistorical analyses and, as a consequence, allow a broader range of gene-flow parameters to be estimated in a precise fashion. More importantly, such nested analyses can discriminate between phylogeographical associations due to recurrent but restricted gene flow vs. historical events operating at the population level (e.g. past fragmentation, colonization, or range expansion events). Restricted gene flow and historical events can be intertwined, and the cladistic analyses can reconstruct their temporal juxtapositions, thereby yielding great insight into both the evolutionary history and population structure of the species. Examples are given that illustrate these properties, concentrating on the detection of range expansion events.     

Comparative phylogeography of three Leptocarabus ground beetle species in South Korea, based on the mitochondrial COI and nuclear 28S rRNA genes

We analyzed the intraspecific gene genealogies of three Leptocarabus ground beetle species (L. seishinensis, L. semiopacus, L. koreanus) in South Korea using sequence data from the mitochondrial cytochrome oxidase subunit I (COI) and nuclear 28S rRNA (28S) genes, and compared phylogeographical patterns among the species. The COI data detected significant genetic differentiation among local populations of all three species, whereas the 28S data showed genetic differentiation only for L. seishinensis. The clearest differentiation of L. seishinensis among local populations was between the northern and southern regions in the COI clades, whereas the 28S clade, which likely indicates relatively ancient events, revealed a range expansion across the northern and southern regions. Leptocarabus semiopacus had the most shallow differentiation of the COI haplotypes, and some clades occurred across the northern and southern regions. In L. koreanus, four diverged COI clades occurred in different regions, with partial overlaps. We discuss the difference in phylogeographical patterns among these Leptocarabus species, as well as between these and other groups of carabid beetles in South Korea.     

Species status and phylogeography of two closely related Coptolabrus species (Coleoptera: Carabidae) in South Korea inferred from mitochondrial and nuclear gene sequences

We investigated the species status and intraspecific phylogeography in South Korea of two ground beetle species, Coptolabrus jankowskii and Coptolabrus smaragdinus (Coleoptera: Carabidae), using statistical parsimony networks and nested clade analyses based on sequences from the mitochondrial cytochrome oxidase subunit I (COI) and nuclear phosphoenolpyruvate carboxykinase (PepCK) and wingless (Wg) genes. Although traditional parsimony tree construction generally failed to resolve interspecific relationships and construct biologically meaningful genealogies, analysis using statistical parsimony networks yielded statistically significant inter- and intraspecific genealogical structures. We found that although these two species represent a notable case of trans-species polymorphisms in both mitochondrial and nuclear gene sequences, their status as separate species was evidenced by the nonrandom association between species and nested clades at various nesting levels. The exceptional occurrence of shared identical or very similar COI sequences was considered to be the result of introgressive hybridization. In addition, range expansion and fragmentation events across the Korean Peninsula and adjacent islands were inferred from nested clade phylogeographical analyses. The COI gene revealed the geographical divergence of major eastern and western clades and historical biogeographical events within each major clade, whereas the nuclear PepCK gene, which did not reveal corresponding east-west clades, indicated past fragmentation and range expansion across wide areas that may have been the result of older biogeographical events. Thus, phylogeographical inferences drawn from analyses of mitochondrial and nuclear genes can reveal different and potentially complementary information about phylogeographical processes.     

Genetic differentiation of relictual populations of Alsophila spinulosa in southern China inferred from cpDNA trnL-F noncoding sequences

The genetic differentiation and phylogeographical pattern of 11 relictual populations of Alsophila spinulosa distributed across Hainan, Guangdong, and Guangxi in southern China were inferred from sequence variations of trnL-F noncoding regions of chloroplast DNA (cpDNA). The length of trnL-F noncoding sequences varied from 863 to 940 bp. The A + T content was 62.23-63.36%. Sequences were neutral in terms of evolution (Tajima's criterion D=-0.62417, P>0.10 and Fu and Li's test D*=-1.45455, P>0.10; F*=-1.32798, P>0.10). Thirty-four haplotypes were identified based on nucleotide variation. Relatively high levels of haplotype diversity (h=0.929) and nucleotide diversity (Dij=0.022263) were detected in A. spinulosa, probably associated with its long evolutionary history which allowed the accumulation of genetic variation within lineages. Both the minimum spanning network and the strict consensus tree of the most parsimonious trees generated for haplotypes demonstrated that the investigated populations of A. spinulosa were subdivided into two geographical groups: Hainan and Guangdong-Guangxi. An analysis of molecular variance (AMOVA) indicated that most of the genetic variation (87.48%, P<0.001) was partitioned among regions. Spatial structure measurements revealed that population genetic structure was not related to geographical distance. This research suggests that blocked gene flow by Qiongzhou strait and an inbreeding system might result in the geographical subdivision between Hainan and Guangdong-Guangxi (F(ST)=0.92, Nm=0.09). Within each region, the "star like" pattern of phylogeography of haplotypes implied a population expansion process during evolutionary history. Gene genealogies together with coalescent theory were useful tools for uncovering the phylogeography of A. spinulosa.     

Comparative phylogeography of the two pink salmon broodlines: an analysis based on a mitochondrial DNA genealogy

Over most of their natural northern Pacific Ocean range, pink salmon (Oncorhynchus gorbuscha) spawn in a habitat that was repeatedly and profoundly affected by Pleistocene glacial advances. A strictly two-year life cycle of pink salmon has resulted in two reproductively isolated broodlines, which spawn in alternating years and evolved as temporal replicates of the same species. To study the influence of historical events on phylogeographical and population genetic structure of the two broodlines, we first reconstructed a fine-scale mtDNA haplotype genealogy from a sample of 80 individuals and then determined the geographical distribution of the major genealogical assemblages for 718 individuals sampled from nine Alaskan and eastern Asian even- and nine odd-year pink salmon populations. Analysis of restriction site states in seven polymerase chain reaction (PCR)-amplified mtDNA regions (comprising 97% of the mitochondrial genome) using 13 endonucleases resolved 38 haplotypes, which clustered into five genealogical lineages that differed from 0.065 to 0.225% in net sequence divergence. The lineage sorting between broodlines was incomplete, which suggests a recent common ancestry. Within each lineage, haplotypes exhibited star-like genealogies indicating recent population growth. The depth of the haplotype genealogy is shallow ( approximately 0.5% of nucleotide sequence divergence) and probably reflects repeated decreases in population size due to Pleistocene glacial advances. Nested clade analysis (NCA) of geographical distances showed that the geographical distribution observed for mitochondrial DNA (mtDNA) haplotypes resulted from alternating influences of historical range expansions and episodes of restricted dispersal. Analyses of molecular variance showed weak geographical structuring of mtDNA variation, except for the strong subdivision between Asian and Alaskan populations within the even-year broodline. The genetic similarities observed among and within geographical regions probably originated from postglacial recolonizations from common sources rather than extensive gene flow. The phylogeographical and population genetic structures differ substantally between broodlines. This can be explained by stochastic lineage sorting in glacial refugia and perhaps different recolonization routes in even- and odd-year broodlines.     

Haplotype diversity of mt<Emphasis Type="Italic">COI</Emphasis> in the fall webworm <Emphasis Type="Italic">Hyphantria cunea</Emphasis> (Lepidoptera: Arctiidae) in introduced regions in China,Iran, Japan,Korea, and its homeland,the United States

Hyphantria cunea (Drury) has colonized many countries outside its native range of North America and has become a model species for studies of the colonization and subsequent adaptation of agricultural pests. Molecular genetic analyses can clarify the origin and subsequent adaptations to non-native habitats. Using the mitochondrial COI gene, we examined the genetic relationships between invasive populations (China, Iran, Japan, and Korea) and native populations (i.e., the United States). The Jilin (China) and Guilan (Iran) populations showed nine previously unknown haplotypes that differed from those found in the south–central United States, suggesting multiple colonization events and different regions of invasion. A dominant mtDNA haplotype in populations in the United States was shared by all of the populations investigated, suggesting that H. cunea with that haplotype have successfully colonized China, Iran, Japan, and Korea.     

Contrasting patterns of historical colonization in white oaks (Quercus spp.) in California and Europe

Phylogeography allows the inference of evolutionary processes that have shaped the current distribution of genealogical lineages across a landscape. In this perspective, comparative phylogeographical analyses are useful in detecting common historical patterns by either comparing different species within the same area within a continent or by comparing similar species in different areas. Here, we analyse one taxon (the white oak, genus Quercus, subgenus Quercus, section Quercus) that is widespread worldwide, and we evaluate its phylogeographical pattern on two different continents: western North America and Western Europe. The goals of the present study are: (i) to compare the chloroplast genetic diversity found in one California oak species vs. that found in the extensively studied European oak species (in France and the Iberian Peninsula); (ii) to contrast the geographical structure of haplotypes between these two taxa and test for a phylogeographical structure for the California species. For this purpose, we used the same six maternally inherited chloroplast microsatellite markers and a similar sampling strategy. The haplotype diversity within site as well as the differentiation among sites was alike in both taxa, but the Californian species has higher allelic richness with a greater number of haplotypes (39 vs. 11 in the European white oak complex). Furthermore, in California these 39 haplotypes are distributed locally in patches while in the European oaks haplotypes are distributed into lineages partitioned longitudinally. These contrasted patterns could indicate that gene movement in California oak populations have been more stable in response to past climatic and geological events, in contrast to their European counterparts.     

Geographic genetic architecture of pocket gopher (Thomomys bottae) populations in Baja California, Mexico

Phylogenetic analyses of complete mitochondrial cytochrome b sequences support the monophyly of pocket gopher (Thomomys bottae) populations from the 1000 km length of the Baja California peninsula of Mexico, relative to other geographical segments of the species range in western North America. The Baja California peninsula is an area that encompasses considerable ecomorphological and infraspecific diversity within this pocket gopher species. However, detailed population analyses encompassing 35 localities distributed over the southern half of the peninsula reveal only trivial phylogeographical structure. Rather, most of the 72 unique 500-base pair haplotypes examined from 142 individuals is restricted to single populations, although a few haplotypes are shared broadly across geography. Individual populations are typically comprised of haplotype sets from different branches in a network of relationships. Analysis of molecular variance (amova) indicates that approximately half of the total pool of variation is contained among individuals within local populations, and that only about 25% can be explained by the regional subdivisions of current subspecies distributions or physiographic realms. A hypothesized historical vicariant event that has been causally linked to the phylogeographical structure of other, codistributed species has had little influence on these pocket gopher populations, explaining only 13% of the total variation. The temporal depth, estimated by coalescence parameters, of the haplotype lineage in Baja California is relatively recent, approximately 300,000 generations; both the mismatch distribution of pairwise comparisons and a significantly positive exponential growth estimate support a recent history of expanding populations; but current, or recent past, migration estimates have remained small, are largely unidirectional from north to south, and weak isolation by distance is present. All data suggest that pocket gophers have relatively recently invaded the southern half of peninsular Baja California, with the genetic signature of expansion still evident but with sufficient time having lapsed to result in a weak isolation by distance pattern. The geographical assemblage of sampled populations thus appears as a meta-population, with limited gene flow contrasting with random haplotype loss due to drift in small, localized populations.     

Biogeographical patterns of genetic differentiation in dung beetles of the genus Trypocopris (Coleoptera, Geotrupidae) inferred from mtDNA and AFLP analyses

Aim To examine the phylogeography and population structure of three dung beetle species of the genus Trypocopris (Coleoptera, Geotrupidae). We wanted to test whether genetic differences and genealogies among populations were in accordance with morphologically described subspecies and we aimed to establish times of divergence among subspecies to depict the appropriate temporal framework of their phylogeographical differentiation. We also wished to investigate the historical demographic events and the relative influences of gene flow and drift on the distribution of genetic variability of the different populations. Location Europe (mostly Italy). Methods We collected adult males from dung pats from 15 Italian localities over the period 2000–2002. For sequence analysis, some dried specimens from Albania, Croatia, Slovakia and Spain were also used. We applied cytochrome oxidase I mitochondrial DNA sequencing and the amplified fragment length polymorphism (AFLP) technique to determine whether phylogeographical patterns within the three species support the proposed hypotheses of subspecies designations, and to detect further structure among populations that might mediate diversification. Results and main conclusions The results show a high concordance between the distribution of mtDNA variation and the main morphological groups recognized as subspecies, which thus may represent independent evolutionary units. The degree of mitochondrial divergence suggests that speciation events occurred during the Pliocene, while diversification of the main subspecific lineages took place in the Pleistocene, from c. 0.3 to 1.5 Ma. Mitochondrial and nuclear data also reveal that there is phylogeographical structuring among populations within each of the main groups and that both contemporary and historical processes determined this pattern of genetic structure. Geographical populations form monophyletic clades in both phylogenetic and network reconstructions. Despite the high levels of intrapopulational diversity, F_ST values indicate moderate but significant genetic differentiation among populations, and a Bayesian clustering analysis of the AFLP data clearly separates the geographical populations. Nucleotide and gene diversity estimates reveal interspecific differences in the degree of diversification among populations that may be related to the different ecological requirements of the three species.     

Phylogeography and the geographic cline in the armament of a seed-predatory weevil: effects of historical events vs. natural selection from the host plant

Japanese camellia (Camellia japonica) and its seed predator, the camellia weevil (Curculio camelliae), provide a notable example of a geographic mosaic of coevolution. In the species interaction, the offensive trait of the weevil (rostrum length) and the defensive trait of the plant (pericarp thickness) are involved in a geographically-structured arms race, and these traits and selective pressures acting on the plant defence vary greatly across a geographical landscape. To further explore the geographical structure of this interspecific interaction, we tested whether the geographical variation in the weevil rostrum over an 800-km range along latitude is attributed to local natural selection or constrained by historical (phylogeographical) events of local populations. Phylogeographical analyses of the mitochondrial DNA sequences of the camellia weevil revealed that this species has experienced differentiation into two regions, with a population bottleneck and subsequent range and/or population expansion within each region. Although these phylogeographical factors have affected the variation in rostrum length, analyses of competing factors for the geographical variation revealed that this pattern is primarily determined by the defensive trait of the host plant rather than by the effects of historical events of populations and a climatic factor (annual mean temperature). Thus, our study suggests the overwhelming strength of coevolutionary selection against the effect of historical events, which may have limited local adaptation.     

Genetic variation and geographical structure of five mouse-sized opossums (Marsupialia, Didelphidae) throughout the Guiana Region

Aim and location To study the patterns of genetic variation and geographical structure of five forest‐dwelling didelphid species distributed throughout the Guiana Region (Marmosops parvidens, M. pinheiroi, Monodelphis brevicaudata, Marmosa murina, Micoureus demerarae) and place the results for wide‐ranging species (M. murina, M. demerarae) in the broader geographical context of the rest of the taxon ranges based on published data. Methods Variation in the mitochondrial cytochrome b gene were analysed in order to assess intra‐specific sequence divergence as well as haplotypic and nucleotide diversity among populations. Relationships between haplotypes were inferred by cladistic (maximum parsimony) and probabilistic (maximum likelihood) approaches, allowing comparisons between tree topologies and patterns of populations geographical distribution. Results Phylogenetic analysis of the Guiana Region populations suggest the lack of common patterns of geographical structure among mouse‐sized opossums except for both Marmosops species. Marmosops parvidens and M. pinheiroi showed individuals from Surinam strongly related with those from Guyana, and haplotypes from French Guiana having a basal position. Concerning Micoureus demerarae, haplotypes from Venezuela were positioned as ancestral, in contrast with Monodelphis brevicaudata, whose French Guiana haplotypes were clustered in a basal position. No evidence of geographical structure was observed for Marmosa murina. Genetic variation within Guiana Region populations fluctuated between 1% for M. murina and 7% in the case of M. brevicaudata. Regarding the widely distributed species, phylogeographical structure at the scale of the Amazon Basin suggests that Guiana Region populations are the sister group of south‐eastern Amazonian and Atlantic Forest populations, thus resulting in an eastern clade well separated from the western Amazonian regions. Main conclusions Monophyly of mouse‐sized opossums haplotypes over the Guiana Region confirms the presence of a well‐differentiated zoogeographical area relative to other Amazonia bioregions. Low genetic diversity among Guiana Region samples suggests a relatively recent origin of populations dating from later Miocene and Pleistocene periods. Moreover, nucleotide and haplotypic diversity values suggest an evolutionary scenario of rapid population growth and dispersal over the Guiana territory, from an ancestral population with small effective size. Dispersal events – rather than vicariance – seem to be responsible for the present genetic and phylogeographical patterns observed in the Guiana Region.     

Phylogeography of the endangered darkling beetle species of Pimelia endemic to Gran Canaria (Canary Islands)

Phylogenetic and geographical nested clade analysis (NCA) methods were applied to mitochondrial DNA sequences of Pimelia darkling beetles (Coleoptera, Tenebrionidae) endemic to Gran Canaria, an island in the Canary archipelago. The three species P. granulicollis, P. estevezi and P. sparsa occur on the island, the latter with three recognized subspecies. Another species, P. fernandezlopezi (endemic to the island of La Gomera) is a close relative of P. granulicollis based on partial Cytochrome Oxidase I mtDNA sequences obtained in a previous study. Some of these beetles are endangered, so phylogeographical structure within species and populations can help to define conservation priorities. A total of about 700 bp of Cytochrome Oxidase II were examined in 18 populations and up to 75 individuals excluding outgroups. Among them, 22 haplotypes were exclusive to P. granulicollis and P. estevezi and 31 were from P. sparsa. Phylogenetic analysis points to the paraphyly of Gran Canarian Pimelia, as the La Gomera P. fernandezlopezi haplotypes are included in them, and reciprocal monophyly of two species groups: one constituted by P. granulicollis, P. estevezi and P. fernandezlopezi (subgenus Aphanaspis), and the other by P. sparsa'sensu lato'. The two species groups show a remarkably high mtDNA divergence. Within P. sparsa, different analyses all reveal a common result, i.e. conflict between current subspecific taxonomic designations and evolutionary units, while P. estevezi and P. fernandezlopezi are very close to P. granulicollis measured at the mtDNA level. Geographical NCA identifies several cases of nonrandom associations between haplotypes and geography that may be caused by allopatric fragmentation of populations with some cases of restriction of gene flow or range expansion. Analyses of molecular variance and geographical NCA allow definition of evolutionary units for conservation purposes in both species-groups and suggest scenarios in which vicariance caused by geological history of the island may have shaped the pattern of the mitochondrial genetic diversity of these beetles.     

Diverse historical processes shape deep phylogeographical divergence in the pollinating seed parasite Greya politella

Understanding the historical framework in which species interactions have diversified across landscapes may help to partition the effects of vicariance and geographically variable selection in shaping the geographical mosaic of coevolving species. We used phylogeographical analyses of the pollinating seed parasite Greya politella (Lepidoptera: Prodoxidae) to define the historical processes that may have structured interactions of this species with its host plants across major biogeographical breaks in western North America. Using 648 bp of cytochrome oxidase I and amplified fragment length polymorphisims, we identified deep genetic breaks among some populations consistent with some definitions of cryptic species. A combination of phylogenetic and population genetic approaches indicates that different historical processes may have structured G. politella genetic diversity in four regions: northern Pacific Northwest, southern Oregon, southern Sierra Nevada, and the remainder of California. The northern Pacific Northwest had high genetic diversity likely due to glacial refugia and subsequent spatial expansion, concordant with some other taxa. Populations in southern Oregon possessed unique, closely related haplotypes with restricted gene flow, possibly indicating a long-standing set of populations in this endemic-rich region. Analyses of California populations showed evidence of restricted gene flow and spatial expansion with many closely related haplotypes that occupy a broad geographical range. Southern Sierra Nevada populations were genetically distinct and highly diverse, possibly due to a localized glacial refugium. Together, these results suggest that vicariance and population expansion, possibly in combination with geographically variable selection, have shaped the diversification of G. politella and its interactions with its host plants.     

Population Genetics of Aedes albopictus (Diptera: Culicidae) Invading Populations, Using Mitochondrial nicotinamide Adenine Dinucleotide Dehydrogenase Subunit 5 Sequences

Aedes albopictus (Skuse) (Diptera: Culicidae), the Asian tiger mosquito indigenous to Asia, now an invasive species worldwide, is an important vector for several arboviruses. Genetic analysis using the mitochondrial nicotinamide adenine dinucleotide dehydrogenase subunit 5 (ND5) gene was carried out in populations from Cameroon (n = 50), Hawaii (n = 38), Italy (n = 20), the continental United States, Brazil, and its native range. Data for Brazil, the continental United States, and the native range was obtained from Birungi and Munstermann (2002). Direct sequencing was used to identity unique haplotypes. The limited phylogeographic partitioning of haplotypes with low levels of sequence divergence in both Cameroon and Hawaii was consistent with the population structure of Ae. albopictus in the United States and Brazil. Four new haplotypes were identified from the samples from Cameroon and Hawaii, adding to previously described haplotypes. Hawaii shared a haplotype with Cameroon that was unique to these two regions. Hawaii also had higher overall haplotype diversity than seen in previous continental United States, Brazil, or native range populations. Hawaiian, Cameroon, and Italian populations did not share haplotypes with Brazil, which validates the earlier mitochondrial DNA studies indicating a separate introduction of this species into Brazil.     

Mapping the suturing of a continental biota

The present-day spatial distribution of interspecific contact zones and intraspecific phylogeographical breaks provides a window into the past ecological and evolutionary processes that underlie speciation and species ranges. The clustering of contact zones and/or phylogeographical breaks in space indicates the suturing of diverged biotas. The presence of such suture zones indicates that similar ecological and historical factors have influenced the past and present distributions of populations and their divergence. Thus, suture zones are ideal natural laboratories for studying divergence, secondary contact and speciation across many different taxa. The concept of suture zones was formalized decades ago by Remington (1968), but only a few detailed and quantitative investigations of suture zones exist (Swenson & Howard 2004, 2005; Whinnett et al. 2005; Moritz et al. 2009). This limited number of investigations is largely because of a lack of detailed geographical data and sophisticated analytical tools. In this issue of Molecular Ecology, Rissler & Smith (2010) have accomplished a detailed investigation into the suturing of amphibian lineages in the United States which uses both detailed geographical data and sophisticated analytical methods. The work greatly enhances our knowledge of suture zones by extending previous work that has focused less on amphibians and by explicitly considering the relationship between species richness and suture zones.     

Variation in mitochondrial COII gene sequences among two species of Japanese knotweed-boring moths, Ostrinia latipennis and O. ovalipennis (Lepidoptera: Crambidae)

The Ostrinia latipennis group contains two species, O. latipennis (Warren) and O. ovalipennis Ohno. These two species commonly utilize perennial knotweeds (Fallopia spp.) as their host plants, which are serious invasive weeds in Europe and North America. Ostrinia latipennis is widely distributed across north-east Asia including Japan whereas O. ovalipennis is restricted to north Japan (Hokkaido Is.) and highland areas of central Japan (Nagano Prefecture in Honshu Is.). To estimate the phylogenetic relatedness and geographical differentiation of the two species, mitochondrial COII gene sequences were determined for specimens covering their distribution ranges in Japan. The uncorrected sequence divergence between O. latipennis and O. ovalipennis was 0.6-0.7%, supporting a close relationship. According to the standard molecular clock proposed for arthropod mtDNA, the two species are speculated to have diverged about 0.3 Myr ago. A single COII gene haplotype was found in O. latipennis irrespective of collection locality. In contrast, two haplotypes were found in O. ovalipennis, and their frequencies were significantly different between the Hokkaido and Honshu populations. The patterns of geographical variation in the COII gene within the two species were in agreement with previously reported patterns of geographical differentiation in morphology of the two species in Japan. The present results support the hypothesis that gene flow among local populations of O. ovalipennis has been limited by geographical isolation.     

High genetic divergence of hormogastrid earthworms (Annelida, Oligochaeta) in the central Iberian Peninsula: evolutionary and demographic implications

Hormogastridae earthworms are highly important for the functioning of the Mediterranean soil system. However, little is known about the species distribution and genetic diversity of these soil invertebrates. In the present study, the genetic differentiation and gene flow were studied among populations of hormogastrids from the central Iberian Peninsula. A 648-bp portion of the mitochondrial cytochrome c oxidase I gene was sequenced for 82 individuals from 7 localities, resulting in the identification of 38 haplotypes exclusive to localities. All of the individuals were morphologically identified as Hormogaster elisae , but the high genetic divergence found among populations (up to 20.20%) suggests the occurrence of more than one cryptic species within this region. Further analysis of the phylogenetic relationships revealed six different evolutionary lineages coincident with geographical location, including the two nearest populations Molar and Redueña as one evolutionary unit. From these results, at least three new species could be inferred, in addition to the morphospecies H. elisae s.s . Partitioning of genetic variance among populations indicated that isolation by distance was the primary agent for differentiation of the investigated hormogastrid populations. Our data suggest that the evolutionary lineages for H. elisae s.l. originated between the late Miocene and the early Pleistocene, but that mtDNA genealogies coalesce on a more recent scale of a few thousand years.     

Geological barriers and restricted gene flow in the holarctic skipper Hesperia comma (Hesperiidae)

Patterns of genetic variation within a species may be a consequence of historical factors, such as past fragmentation, as well as current barriers to gene flow. Using sequence data from the mitochondrial cytochrome oxidase subunit II region (COII) and the nuclear gene wingless, we conducted a phylogeographical study of the holarctic skipper Hesperia comma to elucidate patterns of genetic diversity and to infer historical and contemporary processes maintaining genetic variation. One hundred and fifty-one individuals were sampled from throughout North America and Eurasia, focusing on California and adjacent regions in the western United States where morphological diversity is highest compared to the rest of the range. Analyses of sequence data obtained from both genes revealed a well-supported division between the Old and New World. Within western North America, wingless shows little geographical structure, while a hierarchical analysis of genetic diversity of COII sequences indicates three major clades: a western clade in Oregon and Northern California, an eastern clade including the Great Basin, Rocky Mountains and British Columbia, and a third clade in southern California. The Sierra Nevada and the Transverse Ranges appear to be the major barriers to gene flow for H. comma in the western United States. Relatively reduced haplotype diversity in Eurasia compared to North America suggests that populations on the two continents have been affected by different historical processes.     

Chloroplast phylogeography of Dipentodon (Dipentodontaceae) in southwest China and northern Vietnam

The evolutionary history of plants in the southeast Tibetan Plateau might be the most complicated around the world because of the area's extremely complex topography and climate induced by strong tectonic activity in recent history. In this research, we implemented a phylogeographical study using chloroplast sequences (psbA-trnH and trnQ-rps16 intergenic spacer) on Dipentodon, a monotypic or ditypic genus (D. sinicus and D. longipedicellatus) distributed in southwest China and adjacent areas including Myanmar (Burma), northeast India and northern Vietnam. A total of 257 samples from 16 populations from the southeast Tibetan Plateau (D. longipedicellatus) and the Yungui Plateau (D. sinicus) were collected. The results revealed that Dipentodon had 11 haplotypes for the two intergenic spacers, high genetic diversity (h(T) = 0.902) and high genetic differentiation (N(ST) = 0.987 and G(ST) = 0.948). amova analysis showed that the component of among-population within region/species (55.25%) was unexpectedly larger than the among-species/region component (43.69%), which indicates that there is no justification for recognizing two species in Dipentodon. Correlation of pairwise genetic and geographical distances showed that Dipentodon populations in the southeast Tibetan Plateau may have suffered more habitat fragmentation than populations in the Yungui Plateau because of the uplift of the Tibetan Plateau than populations in the Yungui Plateau have. Nested clade analysis showed that 11 haplotypes formed two 3-level, three 2-level and seven 1-level clades, with eight clades showing significant geographical association. However, clade 2-1 and 2-2 did not cluster together, although they are distributed in the same region (Yungui Plateau) and belong to the same species (D. sinicus). This led not only to incongruence between haplotype network and geographical distribution of 2-level clades, but also to paraphyly of D. sinicus to D. longipedicellatus. We concluded that the incongruence and paraphyly may result from incomplete lineage sorting during the rapid and extreme tectonic events of the Tibetan Plateau. The results reported here will no doubt provide new insights into the evolution of biodiversity on the Tibetan Plateau and adjacent areas, and a historical framework for the conservation of biodiversity in this area, including Dipentodon.