Dawn and dusk "chorus" in visually communicating Jamaican anole lizards

A key hypothesis explaining the existence of dawn and dusk choruses in acoustically communicating animals centers on the need to advertise continued territorial occupancy after and before a period of nocturnal inactivity. If this hypothesis is correct, it follows that similar dawn and dusk choruses should occur in territorial animals exploiting other signal modalities. Adult male Anolis lizards defend territories by using elaborate head-bobbing displays and extensions of a throat fan or dewlap. Males are inactive at night and return to their territories at dawn, remaining there until dusk. I quantified the production of visual displays as a function of time of day for four species on the island of Jamaica: Anolis lineatopus, Anolis sagrei, Anolis grahami, and Anolis opalinus. All exhibited dawn and/or dusk peaks in display behavior. These patterns have remarkable parallels with the dawn and dusk choruses reported for many acoustically communicating animals.     

Habitat light, colour variation, and ultraviolet reflectance in the Grand Cayman anole, Anolis conspersus

Data from a diversity of sources are consistent with the hypothesis that the Grand Cayman anole, Anolis conspersus, is descended directly from Anolis grahami of Jamaica. Although the two species have remained morphologically similar, coloration in A. conspersus has changed considerably from that of its ancestor. The most dramatic difference is seen in dewlap colour, where A. conspersus has evolved a blue and highly UV-reflective dewlap from the ancestral orange-and-yellow colour state. In addition, variation in normal (non-metachrosis) dorsum coloration in A. grahami populations is limited to shades of green (olive, emerald, teal), whereas in A. conspersus dorsum coloration varies from green to blue and to brown. This increased colour variation occurs despite Grand Cayman being a small, relatively featureless island only 35 km in length. Results of this study suggest that ambient light differences associated with precipitation-related vegetation structure may have played an important role in the evolution of A. conspersus body colour variation. Evidence is presented to show how geological, ecological, and physiological factors could have interacted to select for a short wavelength-reflective dewlap from a long wavelength-reflective precursor following the colonization of Grand Cayman from Jamaica by A. grahami between 2 and 3 Mya.     

Escape tactics and effects of perch height and habituation on flight initiation distance in two Jamaican anoles (Squamata: polychrotidae)

Escape by Anolis lizards is influenced by microhabitats and fight initiation distance increases with predation risk. Differences in microhabitat use among ecomorphs affect escape behavior, but only two studies have reported ecomorphological differences in flight initiation distance among Greater Antillean species. I studied effects of predation risk and microhabitats on escape behavior by conducting field experiments using two species of anoles, Anolis lineatopus and A. grahami, on the campus of the University of the West Indies at Mona, Jamaica. Because ecomorphological variation of anoles has evolved independently within each island of the Greater Antilles, but relationships between ecomorphs and escape behaviors are poorly known, I characterized microhabitat use and escape tactics, and determined relationships between flight initiation distance and two risk factors, habituation to human presence and perch height, in Anolis lineatopus, a trunk-ground anole and A. grahami, a trunk-crown anole. Sample sizes for A. lineatopus and A. grahami were 214 and 93, for microhabitat use and escape destinations, 74 and 34 for human presence and 125 and 34 for perch height. The two species occurred in similar microhabitats and exhibited similar escape tactics, but exhibited key differences expected for their ecomorphs. Both species were sighted frequently on the ground and on trees, but A. lineatopus were more frequently on ground and were perched lower than A. grahami. Both species escaped from ground to trees and when on trees hid on far sides and escaped without changing climbing direction with equal frequency. The frequency of fleeing upward was greater for A. grahami than A. lineatopus. Both species exhibited habituation by having shorter flight initiation distances in areas with more frequent exposure to people. In both species flight initiation distance increased as perch height decreased because, lizards had to climb farther to be out of reach when perched lower. The relationship between flight initiation distance and perch height may apply to other anole ecomorphs that flee upward when low perched on trees.     

Phylogeography of Arcterica nana (Ericaceae) suggests another range expansion history of Japanese alpine plants

We conducted a phylogeographic study on the alpine plant Arcterica nana based on haplotypes of chloroplast DNA. Using a sequence of approximately 1,071 bp of intergenic spacers of chloroplast DNA (trnT-L, psbB–psbF), we detected 13 haplotypes among 193 individuals sampled from 22 populations. Two dominant haplotypes were distributed over the entire range of this species in Japan, and we found several local private haplotypes. An analysis of molecular variance (AMOVA) indicated no geographic structure within the haplotype distribution. In addition, the genetic distance was not related to its corresponding geographic distance (Mantel test: r=−0.049, P=0.66), indicating a homogeneous geographic structure throughout the entire distribution range in the Japanese archipelago. The most parsimonious explanation for this geographic structure is that A. nana spread across its extant distribution range in the Japanese archipelago through a recent range expansion event. However, this pattern is inconsistent with the previous phylogeography of Japanese alpine plants, which reveals that haplotypes in central Honshu are differentiated from those in more northern regions. Arcterica nana may have experienced a different history from other alpine plants, suggesting that the history of Japanese alpine flora may include at least two different geographic radiation patterns.     

Phylogenetic analysis of ecological and morphological diversification in Hispaniolan trunk-ground anoles (Anolis cybotes group)

Abstract Anolis lizards in the Greater Antilles partition the structural microhabitats available at a given site into four to six distinct categories. Most microhabitat specialists, or ecomorphs, have evolved only once on each island, yet closely related species of the same ecomorph occur in different geographic macrohabitats across the island. The extent to which closely related species of the same ecomorph have diverged to adapt to different geographic macro-habitats is largely undocumented. On the island of Hispaniola, members of the Anolis cybotes species group belong to the trunk-ground ecomorph category. Despite evolutionary stability of their trunk-ground microhabitat, populations of the A. cybotes group have undergone an evolutionary radiation associated with geographically distinct macrohabitats. A combined phylogeographic and morphometric study of this group reveals a strong association between macrohabitat type and morphology independent of phylogeny. This association results from long-term morphological evolutionary stasis in populations associated with mesic-forest environments ( A. c. cybotes and A. marcanoi ) and predictable morphometric changes associated with entry into new macrohabitat types (i.e., xeric forests, high-altitude pine forest, rock outcrops). Phylogeographic analysis of 73 new mitochondrial DNA sequences (1921 aligned sites) sampled from 68 geographic populations representing 12 recognized species and subspecies diagnoses 16 allopatric or parapatric groupings of populations differing from each other by 5–18% sequence divergence. At least some of these groupings appear to have attained species-level divergence from others. Evolutionary specialization to different macrohabitat types may be a major factor in the evolutionary diversification of Greater Antillean anoles.     

Co-amplification of mitochondrial pseudogenes in Calomys musculinus (Rodentia, Cricetidae): a source of error in phylogeographic studies.

In a previous phylogeographic study of the rodent Calomys musculinus, 24 haplotypes of the mitochondrial DNA D-loop region were detected using the restriction fragment length polymorphism technique (PCR-RFLP). Seven percent of the individuals showed patterns in which the sum of the sizes of the restriction fragments exceeded the size of the original PCR product. In the present paper we analyze possible causes of these atypical haplotypes. PCR products were cloned, and two or three different clones from a single individual were detected by their RFLP patterns. Nine clones with different restriction patterns were selected for sequence analyses. A maximum parsimony phylogenetic analysis revealed two well-supported paraphyletic groups. One group comprised sequences showing low nucleotide divergence compared with the most common haplotypes detected in the phylogeographic study. The other group was basal to the three species of Calomys other than C. musculinus included in the study; the mutations in the short portion of the cytochrome b gene amplified corresponded to 12 amino acid substitutions. The results suggest that two independent insertions of mtDNA sequences into the nucleus occurred; these sequences would co-amplify in the PCR procedure. Identification of pseudogenes is crucial to obtain reliable reconstruction of the intraspecific genealogy in phylogeographic studies.     

Genetic subdivision and biogeography of the Danubian rheophilic barb Barbus petenyi inferred from phylogenetic analysis of mitochondrial DNA variation

The barb Barbus petenyi is a cyprinid widely distributed throughout the mountain regions in the Danube River basin. Phylogenetic analysis of the DNA sequence variation at the mitochondrial cytochrome b gene over much of this range yielded three deep-branching (5.9-9.4% average divergence), well-supported haplotype clades with mutually exclusive geographic distributions and divergence times estimated to be in the Tertiary. The clades did not form an altogether monophyletic group as the most divergent one coalesced more recently with haplotypes of phylogenetically close species than with the other B. petenyi haplotypes. This pattern was supported by bootstrap and log-likelihood Shimodaira-Hasegawa tests. The other two were sister clades, but their distinctiveness was supported by previous allozyme data. Hence, from a taxonomic point of view, the current recognition of B. petenyi is erroneous, as it does not represent a single evolutionary lineage, and we suggest that three species be recognized instead. Substantial phylogeographic differences were evident among the three putative species, the two more southerly ones displaying significant structure, which suggested that they each survived in several glacial refugia throughout the Pleistocene. The phylogeographic pattern of multiple populations of rheophilic barbs with a history of long-term persistence and separation within the Danube River basin is novel within fishes and provides a hypothesis against which phylogeographic patterns among other similarly distributed rheophilic species may be compared.     

Phylogeographic incongruence of codistributed amphibian species based on small differences in geographic distribution

Codistributed species may display either congruent phylogeographic patterns, indicating similar responses to a series of shared climatic and geologic events, or discordant patterns, indicating independent responses. This study compares the phylogeographic patterns of two similarly distributed salamander species within the Pacific Northwest of the United States: Cope's giant salamander (Dicamptodon copei) and Van Dyke's salamander (Plethodon vandykei). Previous studies of P. vandykei support two reciprocally monophyletic lineages corresponding to coastal populations, located from the Olympic Mountains to the mouth of the Columbia River, and inland populations within the Cascade Mountains. We hypothesized that D. copei would have a congruent phylogeographic pattern to P. vandykei due to similarity in distribution and dependence upon similar stream and stream-side habitats. We test this hypothesis by estimating the phylogeny of D. copei using approximately 1800bp of mitochondrial DNA and comparing it to that of P. vandykei. Sympatric populations of D. copei and of P. vandykei display an identical phylogeographic pattern, suggesting similar responses within their shared distribution. Populations of D. copei occurring outside the range of P. vandykei displayed high levels of genetic divergence from those sympatric to P. vandykei. Overall, phylogeographic patterns between the two species were ultimately incongruent due to the high divergence of these allopatric populations. These results provide an example of codistributed species displaying overall incongruent phylogeographic patterns while simultaneously displaying congruent patterns within portions of their shared geographic distribution. This pattern demonstrates that a simple dichotomy of congruent and incongruent phylogeographic patterns of codistributed species may be too simplistic and that more complex intermediate patterns can result even from minor differences in species' ranges.     

Deciphering the global phylogeography of a coastal shrub (Scaevola taccada) reveals the influence of multiple forces on contemporary population structure

The phylogeography of coastal plant species is heavily influenced by past sealevel fluctuations, dispersal barriers, and life-history traits, such as long-distance dispersal ability of the propagules. Unlike the widely studied mangroves, phylogeographic patterns have remained mostly obscure for other coastal plant species. In this study, we sampled 42 populations of Scaevola taccada (Gaertn.) Roxb., a coastal shrub of the family Goodeniaceae, from 17 countries across its distribution range. We used five chloroplast DNA (cpDNA) and 14 nuclear microsatellite (simple sequence repeat [SSR]) markers to assess the influence of abiotic factors and population genetic processes on the phylogeographic pattern of the species. Geographical distribution of cpDNA haplotypes suggests that the species originated in Australia, followed by historical dispersal and expansion of its geographic range. Multiple abiotic factors, including the sealevel changes during the Pleistocene, the presence of landmasses like the Malay Peninsula, and contemporary oceanic circulation patterns, restricted gene flow between geographically distinct populations, thereby creating low haplotype diversity and a strong population structure. Population genetic processes acted on these isolated populations, leading to high nuclear genetic diversity and population differentiation, as revealed from analyzing the polymorphic SSR loci. Although genetic divergence was mostly concordant between cpDNA and SSR data, asymmetrical gene flow and ancestral polymorphism could explain the discordance in the detailed genetic structure. Overall, our findings indicate that abiotic factors and population genetic processes interactively influenced the evolutionary history and current phylogeographic pattern of S. taccada across its distribution range.     

Genetic diversity in Australian ancient asexual <Emphasis Type="Italic">Vestalenula</Emphasis> (Ostracoda,Darwinulidae): little variability down under

Darwinulid ostracods are putative ancient asexuals, and are thus assumed to be unable to purge deleterious mutations from their genomes. Some darwinulids species can be found both above (epigeic) and below ground (hypogeic). We hypothesize that surface populations carry more mutations than their below-ground counterparts, which are buffered from mutagens such as UV-B. Given the age of the investigated area, the Pilbara in Western Australia, we also expect geographic patterning of observed haplotypes. We have used DNA sequence data from the nuclear ITS and the mitochondrial COI region to investigate a (limited) data set on two Australian species, the endemic Vestalenula matildae and V. marmonieri from the Pilbara region. We do not find differences in genetic variability between specimens from subterranean habitats as compared to those from habitats above ground. There was also no congruence between hydrological basins and distribution patterns of the haplotypes identified. Although our data indicate that the two species may have split from each other ca. 70 myr ago, this has not resulted in any clear phylogeographic patterns among the analysed specimens across the regions of the Pilbara.     

Origin of the disjunct tetraploid Cardamine amporitana (Brassicaceae) assessed with nuclear and chloroplast DNA sequence data

Seventy-four nucleotide sequences from the ITS regions of nuclear ribosomal DNA and 76 from the trnL-trnF spacer of chloroplast DNA were used to address the origin of tetraploid Cardamine amporitana, the conspecifity of central Italian and northeastern Spanish populations, and the possible cause for such geographic disjunction. Because of the complex lineage relationships in Cardamine, the sampling included 22 taxa. In the results, both data sets are highly congruent in supporting a close relationship of C. amporitana to the widespread Eurasian C. amara. Low genetic variability in northeastern Spanish populations of C. amporitana suggests long-distance dispersal from central Italy. The interior position of the single northeastern Spanish haplotype in a statistical parsimony network of trnL-trnF haplotypes however does not support this scenario and invokes other plausible phylogeographic explanations. The disappearance of geographically intermediate populations and genetic impoverishment by migration and isolation, both probably associated with Quaternary climatic oscillations, appears as an alternative hypothesis to explain the phylogeographic pattern. A recent hybridization event is reported between C. amporitana and a diploid from the C. pratensis group in central Italy on the basis of additive polymorphisms in ITS for all the 22 distinguishing nucleotides.     

Comparative phylogeography of red maple (Acer rubrum L.) and silver maple (Acer saccharinum L.): impacts of habitat specialization,hybridization and glacial history

Aim We analysed variation in chloroplast DNA (cpDNA) in red maple (Acer rubrum L.) and silver maple (Acer saccharinum L.) across a large part of their geographic ranges. Acer rubrum is one of the most common and morphologically variable deciduous trees of eastern North America, while its sister species A. saccharinum has a more restricted habitat distribution and displays markedly less morphological variation. Our objective was to infer the impact of biogeographic history on cpDNA diversity and phylogeographic structure in both species. Location Deciduous forests of eastern North America. Methods We sequenced 1289 to 1645 bp of non‐coding cpDNA from A. rubrum (n = 258) and A. saccharinum (n = 83). Maximum parsimony networks and spatial analysis of molecular variance (SAMOVA) were used to analyse phylogeographic structure. Rarefaction analyses were used to compare genetic diversity. Results A total of 40 cpDNA haplotypes were recovered from A. rubrum (38 haplotypes) and A. saccharinum (7 haplotypes). Five of the seven A. saccharinum haplotypes were shared with nearby samples of A. rubrum. SAMOVA recovered four phylogeographic groups for A. rubrum in: (1) south‐eastern USA, (2) the Gulf and south‐eastern Coastal Plain, (3) the lower Mississippi River Valley, and (4) the central and northern regions of eastern North America. Acer saccharinum had significantly lower haplotype diversity than A. rubrum, and novel haplotypes in post‐glaciated northern limits of its range were shared with A. rubrum. Main conclusions This is the first study of A. rubrum to report a distinct phylogeographic group centred on the lower Mississippi River, and the first to examine data comparatively with A. saccharinum. We hypothesized that A. rubrum would display stronger phylogeographic structure and greater haplotype diversity than A. saccharinum because of its greater geographic range, and ecological and morphological variation. This hypothesis was supported by the cpDNA analysis. The sharing of cpDNA and chloroplast simple sequence repeat (cpSSR) haplotypes in areas of geographic overlap provides evidence of introgression, which led to an increase in haplotype diversity in both species, and to novel phylogeographic structure in A. rubrum. We recommend that introgression be considered, along with other potential causes, as an explanation for the phylogeographic structure of cpDNA in plants.     

Geographic pattern of genetic variation in the fox tapeworm Echinococcus multilocularis

Intraspecific genetic variation of Echinococcus multilocularis, the etiologic agent of human alveolar echinococcosis, has been evaluated among 76 geographic isolates from Europe, Asia and North America by using sequence data of mitochondrial and nuclear DNA. Relatively low genetic variation was found only in the mitochondrial DNA sequence consisting of 3 protein-coding genes. Pairwise divergence among the resultant 18 haplotypes ranged from 0.03 to 1.91%. Phylogenetic trees and parsimony network of these haplotypes depicted a geographic division into European, Asian and North American clades, but 1 haplotype from Inner Mongolia was unrelated to other haplotypes. The coexistence of the Asian and North American haplotypes could be seen, particularly on the St. Lawrence Island in the Bering Sea. These data suggest an evolutionary scenario in which distinct parasite populations derived from glacial refugia have been maintained by indigenous host mammals. The nuclear DNA sequence for the immunodominant B cell epitope region of ezrin/radixin/moesin-like protein (elp) was extremely conservative, indicating that the elp antigen is available for immunodiagnosis in any endemic areas.     

Phylogeographic structure of the fire ant Solenopsis invicta in its native South American range: roles of natural barriers and habitat connectivity

We generated mitochondrial DNA (mtDNA) sequence data from 402 individuals of the fire ant Solenopsis invicta collected from 11 native populations and analyzed these data using a combination of demographic, phylogenetic, and phylogeographic methods to infer features of the evolutionary history of this species. Prior expectations regarding high levels of genetic structure and isolation by distance among populations were supported by the data, but we also discovered several unanticipated patterns. Our analyses revealed a major genetic break between S. invicta mtDNA haplotypes that coincides with the Mesopotamia wetlands region of South America, resulting in two higher level nested clade groupings. In addition, we identified contrasting patterns of genetic differentiation within these two major groups, which may reflect differences in connectivity of suitable habitat in different parts of the native range of S. invicta. Our study represents the first attempt to understand the phylogeographic history of S. invicta across its native range.     

Are subspecies of Anolis lizards that differ in dewlap color and pattern also genetically distinct? A mitochondrial analysis

Subspecies of Anolis lizards are often defined on the basis of geographic variation in the color and pattern of the dewlap, an extensible throat fan considered central to species recognition and sexual selection. Among the most impressive examples of this phenomenon are two species of trunk anoles found across Hispaniola and the Bahamas: Anolis distichus is divided into 16 subspecies with dewlap colors ranging from deep wine red to pale yellow while Anolis brevirostris is divided into three subspecies with dewlaps ranging from pale yellow to orange. Limited sampling of allozyme data indicates some genetic divergence among subspecies and suggests that they may deserve recognition at the species-level. Our goal here is to use more comprehensive geographic sampling of mtDNA haplotypes to test whether the five subspecies of A. distichus and three subspecies of A. brevirostris that occur in the Dominican Republic correspond with genetically distinct populations that may warrant recognition under the general lineage concept. We obtain an aligned dataset of 1462bp comprised of the genes encoding ND2 and adjacent tRNAs from 76 individuals of A. distichus from 28 localities and 12 individuals of A. brevirostris from five localities. We find that haplotypes sampled from each Dominican subspecies of A. distichus form well-supported and deeply divergent clades (>10% uncorrected sequence divergence). Strong concordance between mtDNA haplotype structure and previously diagnosed phenotypic variation in traits central to interspecific communication (i.e., the dewlap) leads us to hypothesize that each of the presently recognized Dominican subspecies of A. distichus and A. brevirostris deserves elevation to full species status under the general lineage concept.     

A molecular phylogenetic analysis of diversification in Amazonian Anolis lizards

We present a mitochondrial DNA (mtDNA) haplotype phylogeny for Amazonian Anolis lizards, including geographical sampling within four species distributed across the Amazon basin (A. fuscoauratus, A. nitens, A. ortonii and A. punctatus). Approximately 1500 bp of mtDNA encoding ND2, COI and four transfer RNAs (tRNAs) are reported for 39 specimens representing four to five populations of each widespread species, plus eight outgroups. These new sequences are aligned with eight previously published sequences, yielding 914 variable characters and 780 parsimony-informative characters. Phylogenetic analyses using maximum parsimony and maximum likelihood reject the hypothesis that Amazonian anoles form a monophyletic group excluding Central American and Caribbean anoles, and suggest multiple faunal exchanges among these regions. Haplotype divergence among geographical populations of A. nitens, whose variation was influential in formulating the Pleistocene refuge hypothesis of Amazonian speciation, is very large (13-22% sequence difference), suggesting that these populations separated well before the Pleistocene. Haplotype divergences among geographical populations of A. fuscoauratus (3-4%), A. punctatus (4-9%) and A. ortonii (6-8%) also indicate pre-Pleistocene differentiation within each species, but temporally incongruent patterns among species.     

Phylogeography and population structure of the Atlantic and Mediterranean green turtle Chelonia mydas: a mitochondrial DNA control region sequence assessment

Mitochondrial (mt) DNA sequences were analysed to resolve the phylogeography and population genetic structure of Atlantic and Mediterranean populations of green turtles ( Chelonia mydas ). Analysis of sequence variation over 487 base pairs of the control (D-loop) region identified 18 haplotypes among 147 individuals from nine nesting populations. Pairwise comparisons of haplotype frequencies distinguished most nesting colonies, indicating significant genetic differentiation among rookeries and a strong propensity for natal homing behaviour by nesting females. Comparison of control region sequence data to earlier restriction fragment length polymorphism (RFLP) data for the same individuals demonstrates approximately a sixfold higher substitution rate in the 5' end of the control region. The sequence data provide higher resolution both in terms of the number of mtDNA genotype variants and the phylogeographic relationships detected within the Atlantic region, and reveal a gene genealogy that distinguishes two groups of haplotypes corresponding to (i) the western Caribbean and Mediterranean, and (ii) eastern Caribbean, South Atlantic and West Africa. The data suggest that phylogeographic patterns in the Atlantic Ocean may be interpreted in terms of female nest site fidelity and episodic dispersal events. The distribution of mtDNA haplotypes within the region is thus explained by the geological and climatic alternations (glacial and interglacial) over the last million years.     

Phylogeography of a widespread lizard complex reflects patterns of both geographic and ecological isolation

A primary challenge for modern phylogeography is understanding how ecology and geography, both contemporary and historical, shape the spatial distribution and evolutionary histories of species. Phylogeographic patterns are the result of many factors, including geology, climate, habitat, colonization history and lineage‐specific constraints. Assessing the relative influences of these factors is difficult because few species, regions and environments are sampled in enough detail to compare competing hypotheses rigorously and because a particular phylogeographic pattern can potentially result from different evolutionary scenarios. The silky anoles (Anolis sericeus complex) of Central America and Mexico are abundant and found in all types of lowland terrestrial habitat, offering an excellent opportunity to test the relative influences of the factors affecting diversification. Here, we performed a range‐wide statistical phylogeographic analysis on restriction site‐associated DNA (RAD) markers from silky anoles and compared the phylogeographic patterns we recovered to historical and contemporary environmental and topographic data. We constructed niche models to compare niche overlap between sister lineages and conducted coalescent simulations to characterize how the major lineages of silky anoles have diverged. Our results revealed that the mode of divergence for major lineage diversification events was geographic isolation, resulting in ecological divergence between lineages, followed by secondary contact. Moreover, comparisons of parapatric sister lineages suggest that ecological niche divergence contributed to isolation by environment in this system, reflecting the natural history differences among populations in divergent environments.     

Phylogeography and genetic diversity of East Asian Neolitsea sericea (Lauraceae) based on variations in chloroplast DNA sequences

Neolitsea sericea is an evergreen broad leaved tree in the warm-temperate regions of East Asia. This area is a hotspot for plant species richness and endemism caused by dynamic changes in land configuration during the Quaternary. However, the historical migration of such evergreen tree species is still poorly understood. In an attempt to reconstruct the phylogeographic history of N. sericea during the Quaternary, we identified the chloroplast DNA haplotypes of 287 individuals from 33 populations covering almost all of its geographic range. Analyses were based on sequence data from the trnL–F, psbC–trnS, and rps16 regions. Nine haplotypes were identified. The majority included ancestral types in the southwestern part of the main islands of Japan, with other region-specific haplotypes being found in populations on the Korean Peninsula, Taiwan (Isl. Lanyu), and elsewhere in Japan. A statistical parsimony network revealed two lineages derived from Japanese main islands. One was represented on the Korean Peninsula, the other on Isl. Lanyu. The current distribution of N. sericea has been shaped by colonization via land bridges. During the glacial periods, two primary, but separate migration routes were followed—from the southwestern part of the Japanese main islands to either the Korean Peninsula or Taiwan. In addition, we believe the Zhoushan populations were shaped by post-glacial processes through an ECS land bridge (East China Sea basin) from northern refugia that existed during the late Pleistocene.     

Distinct Phylogeographic Structures of Wild Radish (Raphanus sativus L. var. raphanistroides Makino) in Japan

Coastal plants with simple linear distribution ranges along coastlines provide a suitable system for improving our understanding of patterns of intra-specific distributional history and genetic variation. Due to the combination of high seed longevity and high dispersibility of seeds via seawater, we hypothesized that wild radish would poorly represent phylogeographic structure at the local scale. On the other hand, we also hypothesized that wild radish populations might be geographically differentiated, as has been exhibited by their considerable phenotypic variations along the islands of Japan. We conducted nuclear DNA microsatellite loci and chloroplast DNA haplotype analyses for 486 samples and 144 samples, respectively, from 18 populations to investigate the phylogeographic structure of wild radish in Japan. Cluster analysis supported the existence of differential genetic structures between the Ryukyu Islands and mainland Japan populations. A significant strong pattern of isolation by distance and significant evidence of a recent bottleneck were detected. The chloroplast marker analysis resulted in the generation of eight haplotypes, of which two haplotypes (A and B) were broadly distributed in most wild radish populations. High levels of variation in microsatellite loci were identified, whereas cpDNA displayed low levels of genetic diversity within populations. Our results indicate that the Kuroshio Current would have contributed to the sculpting of the phylogeographic structure by shaping genetic gaps between isolated populations. In addition, the Tokara Strait would have created a geographic barrier between the Ryukyu Islands and mainland Japan. Finally, extant habitat disturbances (coastal erosion), migration patterns (linear expansion), and geographic characteristics (small islands and sea currents) have influenced the expansion and historical population dynamics of wild radish. Our study is the first to record the robust phylogeographic structure in wild radish between the Ryukyu Islands and mainland Japan, and might provide new insight into the genetic differentiation of coastal plants across islands.