Phylogeography of the temperate grassland plant Tephroseris kirilowii (Asteraceae) inferred from multiplexed inter-simple sequence repeat genotyping by sequencing (MIG-seq) data.

A group of temperate grassland plant species termed the "Mansen elements" occurs in Japan and is widely distributed in the grasslands of continental East Asia. It has been hypothesized that these species are continental grassland relicts in Japan that stretch back to a colder age, but their migration history has not been elucidated. To assess the migration history of the Mansen elements, we performed phylogeographic analyses of Tephroseris kirilowii, a member of this group, using single-nucleotide polymorphisms (SNPs) obtained from multiplexed inter-simple sequence repeat genotyping by sequencing (MIG-seq). It was estimated that the Japanese populations of T. kirilowii were divided from those of continental East Asia at 25.2 thousand years ago (ka) with 95% highest probability density interval (HPD) of 15.3-40.0 ka and that Japanese clades first diverged at 20.2 ka with 95% HPD of 10.4-30.1 ka. As the climatically suitable range during the last glacial maximum (LGM) estimated using ecological niche modeling (ENM) was limited in Japan and there was a slight genetic differentiation among Japanese populations, a post-glacial expansion of T. kirilowii in the Japanese Archipelago was indicated.

Effects of fruit dimorphism on genetic structure and gene flow in the coastal shrub Scaevola taccada.

BACKGROUND AND AIMS: Plant propagules often possess specialized morphologies that facilitate dispersal across specific landscapes. In the fruit dimorphism of a coastal shrub, Scaevola taccada, individual plants produce either cork-morph or pulp-morph fruits. The former is buoyant and common on sandy beaches, whereas the latter does not float, is bird-dispersed, and is common on elevated sites such as slopes on sea cliffs and behind rocky shores. We hypothesized that beach populations bridge the heterogeneous landscapes by serving as a source of both fruit types, while dispersal is biased for the pulp morph on elevated sites within the islands and for the cork morph between beaches of different islands. Based on this hypothesis, we predicted that populations in elevated sites would diverge genetically over time due to isolation by distance, whereas beach populations would maintain high genetic similarity via current gene flow. METHODS: The genetic structure and gene flow in S. taccada were evaluated by investigating genome-wide single nucleotide polymorphisms in plants from 17 sampling sites on six islands (belonging to the Ryukyu, Daito and Ogasawara Islands) in Japan. KEY RESULTS: Geographical isolation was detected among the three distant island groups. Analyses within the Ryukyu Islands suggested that sandy beach populations were characterized by genetic admixture, whereas populations in elevated sites were relatively isolated between the islands. Pairwise FST values between islands were lowest between sandy beaches, intermediate between sandy beaches and elevated sites, and highest between elevated sites. CONCLUSIONS: Dispersal across the ocean by cork morphs is sufficiently frequent to prevent genetic divergence between beaches of different islands. Stronger genetic isolation of elevated sites between islands suggests that bird dispersal by pulp morphs is restricted mainly within islands. These contrasting patterns of gene flow realized by fruit dimorphism provide evidence that fruit characteristics can strongly mediate genetic structure.

Phylogeography of a canopy-forming kelp, Eisenia bicyclis (Laminariales, Phaeophyceae), based on a genome-wide sequencing analysis.

Analyses of phylogeographic patterns and genetic diversity provide fundamental information for the management and conservation of species. However, little is published about these patterns in Japanese kelp species. In this study, we conducted phylogeographic analyses of a canopy-forming kelp, Eisenia bicyclis, based on genome-wide SNPs identified by ddRAD-seq. We obtained 1,299 SNPs for 76 samples from nine localities across the distribution. STRUCTURE, NeighborNet, and discriminant analysis of principal components consistently showed high genetic differentiation among the Eastern Pacific, Central Pacific, and Sea of Japan coastal regions. Relatively strong gene flow was detected only within populations in the Eastern Pacific and in the Sea of Japan. Genetic diversity and genetic uniqueness were high in the Central Pacific and low in the Sea of Japan. These results suggest that there were at least three independent refugia corresponding to the three regions during the Last Glacial Maximum (LGM). Furthermore, relatively larger populations in the Central Pacific and smaller populations in the Sea of Japan have been maintained in the demographic history from before the LGM to the present. These phylogeographic histories were supported by an Approximate Bayesian Computation analysis. From a conservation genetics perspective, the loss of southern populations in the Central Pacific would greatly reduce the total genetic diversity of the species. Southern populations in the Sea of Japan, which have relatively low genetic diversity, may be highly vulnerable to environmental change, such as heat waves and increased feeding. Therefore, careful monitoring and conservation are needed in the two regions.

Phylogeography of the East Asian grassland plant, Viola orientalis (Violaceae), inferred from plastid and nuclear restriction site-associated DNA sequencing data.

To elucidate the origin and migration history of the "Mansen elements," a group of temperate grassland plants mainly distributed in northeastern Asia, phylogeographic analyses based on chloroplast DNA markers and double-digest restriction site-associated DNA sequencing (ddRAD-seq) data were performed on Viola orientalis, one of the representative species of the group. Phylogenetic analyses using ddRAD-seq data revealed that the populations of V. orientalis were clustered into five clades, among which the continental clades made of populations from Russia and Korea diverged more than 100,000 years earlier than the Japanese clades. The Japanese clade likely diverged during the last glacial period, followed by a further post-glacial divergence into the Kyushu and the Honshu subclades. Our study demonstrated that V. orientalis originated in the continental area of northeastern Asia and, during the last glacial period, has spread southward through the Korean Peninsula across the Japanese Islands. This finding supports the previously proposed evolutionary hypothesis regarding the origin and migration routes of the Mansen elements.

Geographical co-occurrence of butterfly species: the importance of niche filtering by host plant species.

The relevance of interspecific resource competition in the context of community assembly by herbivorous insects is a well-known topic in ecology. Most previous studies focused on local species assemblies that shared host plants. Few studies evaluated species pairs within a single taxon when investigating the effects of host plant sharing at the regional scale. Herein, we explore the effect of plant sharing on the geographical co-occurrence patterns of 232 butterflies distributed across the Japanese archipelago; we use two spatial scales (10 × 10 and 1 × 1 km grids) to this end. We considered that we might encounter one of two predictable patterns in terms of the relationship between co-occurrence and host sharing among butterflies. On the one hand, host sharing might promote distributional exclusivity attributable to interspecific resource competition. On the other hand, sharing of host plants may promote co-occurrence attributable to filtering by resource niche. At both grid scales, we found significant negative correlations between host use similarity and distributional exclusivity. Our results support the hypothesis that the butterfly co-occurrence pattern across the Japanese archipelago is better explained by filtering via resource niche rather than interspecific resource competition.

A Genome Scan for Genes Underlying Microgeographic-Scale Local Adaptation in a Wild Arabidopsis Species.

Adaptive divergence at the microgeographic scale has been generally disregarded because high gene flow is expected to disrupt local adaptation. Yet, growing number of studies reporting adaptive divergence at a small spatial scale highlight the importance of this process in evolutionary biology. To investigate the genetic basis of microgeographic local adaptation, we conducted a genome-wide scan among sets of continuously distributed populations of Arabidopsis halleri subsp. gemmifera that show altitudinal phenotypic divergence despite gene flow. Genomic comparisons were independently conducted in two distinct mountains where similar highland ecotypes are observed, presumably as a result of convergent evolution. Here, we established a de novo reference genome and employed an individual-based resequencing for a total of 56 individuals. Among 527,225 reliable SNP loci, we focused on those showing a unidirectional allele frequency shift across altitudes. Statistical tests on the screened genes showed that our microgeographic population genomic approach successfully retrieve genes with functional annotations that are in line with the known phenotypic and environmental differences between altitudes. Furthermore, comparison between the two distinct mountains enabled us to screen out those genes that are neutral or adaptive only in either mountain, and identify the genes involved in the convergent evolution. Our study demonstrates that the genomic comparison among a set of genetically connected populations, instead of the commonly-performed comparison between two isolated populations, can also offer an effective screening for the genetic basis of local adaptation.

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.

The role of oceanic currents in the dispersal and connectivity of the mangrove Rhizophora mangle on the Southwest Atlantic region.

Dispersal is a crucial mechanism to living beings, allowing them to reach new resources such that populations and species can occupy new environments. However, directly observing the dispersal mechanisms of widespread species can be costly or even impractical, which is the case for mangrove trees. The influence of ocean currents on mangrove dispersal is increasingly evident; however, few studies mechanistically relate the patterns of population distribution with the dispersal by oceanic currents under an integrated framework. Here, we evaluate the role of oceanic currents on connectivity of Rhizophora mangle along the Southwest Atlantic. We inferred population genetic structure and migration rates, simulated the displacement of propagules and tested our hypotheses with Mantel tests and redundancy analysis. We observed populations structured in two major groups, north and south, which is corroborated by other studies with Rhizophora and other coastal plants. Inferred recent migration rates do not indicate ongoing gene flow between sites. Conversely, long-term migration rates were low across groups and contrasting dispersal patterns within each one, which is consistent with long-distance dispersal events. Our hypothesis tests suggest that both isolation by distance and isolation by oceanography (derived from the oceanic currents) can explain the neutral genetic variation of R. mangle in the region. Our findings expand current knowledge of mangrove connectivity and highlight how the association of molecular methods with oceanographic simulations improve the interpretation of the dispersal process. This integrative approach is a cost- and time-efficient strategy to include dispersal and connectivity data into marine protected areas planning and management.

Comparative phylogeography of four component species of deciduous broad-leaved forests in Japan based on chloroplast DNA variation.

A phylogeographic study of four tree species (Padus grayana, Euonymus oxyphyllus, Magnolia hypoleuca, and Carpinus laxiflora) growing in Japanese deciduous broad-leaved forests was conducted based on chloroplast DNA (cpDNA) variations. Using nucleotide sequences of 702-1,059 bp of intergenic spacers of cpDNA, 20, 27, eight, and eight haplotypes were detected among 251, 251, 226, and 262 individuals sampled from 67, 79, 75, and 71 populations of the above species, respectively. The geographical pattern of the cpDNA variations was highly structured in each species, and the following three regional populations were genetically highly differentiated among all four species: (1) the Sea of Japan-side area, (2) the Kanto region, and (3) southwestern Japan. Based on some interspecific similarities among the phylogeographic patterns, the following migration scenario of Japanese deciduous broad-leaved forests was postulated. During the last glacial maximum (LGM), the forests were separately distributed in six regions. After LGM, as the climate warmed, the forests in eastern Japan separately expanded from each of the refugia along the Sea of Japan-side or along the Pacific Ocean-side. In contrast, those in southwestern Japan retreated and moved to high altitudes from each of the continuous forests.

Biogeographic Patterns of Ectomycorrhizal Fungal Communities Associated With Castanopsis sieboldii Across the Japanese Archipelago.

Biogeographic patterns in ectomycorrhizal (ECM) fungal communities and their drivers have been elucidated, including effects of host tree species and abiotic (climatic and edaphic) conditions. At these geographic scales, genotypic diversity and composition of single host tree species change with spatial and environmental gradients, reflecting their historical dispersal events. However, whether the host genotypes can be associated with the biogeographic patterns of ECM communities remains unclear. We investigated the biogeographic pattern of ECM fungal community associated with the single host species Castanopsis sieboldii (Fagaceae), whose genotypic diversity and composition across the Japanese archipelago has already been evaluated. ECM communities were investigated in 12 mature Castanopsis-dominated forests covering almost the entire distribution range of C. sieboldii, and we quantified the effect of host genotypes on the biogeographic pattern of ECM fungal communities. Richness and community composition of ECM fungi changed with latitude and longitude; these biogeographic changes of ECM community were significantly correlated with host genotypic variables. Quantitative analyses showed a higher relative explanatory power of climatic and spatial variables than that of host genotypic variables for the biogeographic patterns in the ECM community. Our results suggest historical events of host dispersal can affect the biogeographic patterns of the ECM fungal community, while their explanation power was lower than that for climatic filtering and/or fungal dispersal.

ddRAD-seq based phylogeographic study of Sargassum thunbergii (Phaeophyceae, Heterokonta) around Japanese coast.

Using genome-wide SNP data obtained from high-throughput techniques based on double digest restriction site-associated DNA sequencing (ddRAD-seq), we elucidated the migration history and genetic diversity of the Japanese population of the ecologically important brown seaweed Sargassum thunbergii (Mertens ex Roth) Kuntze. STRUCTURE and NeighborNet analyses showed a clear genetic differentiation among populations of four geographic regions: Kyushu (POP1); Sea of Japan (POP2); Hokkaido and Tohoku (POP3); and Pacific coast from Kyushu to Kanto (POP4). Approximate Bayesian Computation (ABC) analysis indicated that POP4 diverged first, followed by the separation between POP2 (the largest effective population size) and POP3; POP1 was the last to form, shaped by the mixture of POP2 (73%) and POP4 (27%). High genetic diversity was detected in POP1 and POP2, whereas low genetic diversity was detected in POP3 and POP4. These results indicated that S. thunbergii populations of Kyushu and the Sea of Japan might have been maintained as large and stable populations gathered different lineages from China, Korea and/or Japan.

Intraspecific sequence variation of chloroplast DNA among the component species of deciduous broad-leaved forests in Japan.

To select appropriate plant materials for a phylogeography of deciduous broad-leaved forests in Japan, we surveyed intraspecific chloroplast DNA variation in 34 species found in these forests. A relatively large number of intraspecific cpDNA variations were detected in ten species: Carpinus japonica (nucleotide diversity pi=0.00083), C. laxiflora (pi=0.00221), Magnolia obovata (pi=0.00134), Lindera triloba (pi=0.00255), L. obtusiloba (pi=0.00289), Pourthiaea villosa var. leavis (pi=0.00263), Acer japonicum (pi=0.00170), A. micranthum (pi=0.00237), Euonymus oxyphyllus (pi=0.00322) and Styrax obassia (pi=0.00100).

Earthquake source fault beneath Tokyo.

Devastating earthquakes occur on a megathrust fault that underlies the Tokyo metropolitan region. We identify this fault with use of deep seismic reflection profiling to be the upper surface of the Philippine Sea plate. The depth to the top of this plate, 4 to 26 kilometers, is much shallower than previous estimates based on the distribution of seismicity. This shallower plate geometry changes the location of maximum finite slip of the 1923 Kanto earthquake and will affect estimations of strong ground motion for seismic hazards analysis within the Tokyo region.

High pore fluid pressure may cause silent slip in the Nankai Trough.

Silent-slip events have been detected at several subduction zones, but the cause of these events is unknown. Using seismic imaging, we detected a cause of the Tokai silent slip, which occurred at a presumed fault zone of a great earthquake. The seismic image that we obtained shows a zone of high pore fluid pressure in the subducted oceanic crust located down-dip of a subducted ridge. We propose that these structures effectively extend a region of conditionally stable slips and consequently generate the silent slip.