Divergent mechanisms of reduced growth performance in Betula ermanii saplings from high-altitude and low-latitude range edges.

The reduced growth performance of individuals from range edges is a common phenomenon in various taxa, and considered to be an evolutionary factor that limits the species' range. However, most studies did not distinguish between two mechanisms that can lead to this reduction: genetic load and adaptive selection to harsh conditions. To address this lack of understanding, we investigated the climatic and genetic factors underlying the growth performance of Betula ermanii saplings transplanted from 11 populations including high-altitude edge and low-latitude edge population. We estimated the climatic position of the populations within the overall B. ermanii's distribution, and the genetic composition and diversity using restriction-site associated DNA sequencing, and measured survival, growth rates and individual size of the saplings. The high-altitude edge population (APW) was located below the 95% significance interval for the mean annual temperature range, but did not show any distinctive genetic characteristics. In contrast, the low-latitude edge population (SHK) exhibited a high level of linkage disequilibrium, low genetic diversity, a distinct genetic composition from the other populations, and a high relatedness coefficient. Both APW and SHK saplings displayed lower survival rates, heights and diameters, while SHK saplings also exhibited lower growth rates than the other populations' saplings. The low heights and diameters of APW saplings was likely the result of adaptive selection to harsh conditions, while the low survival and growth rates of SHK saplings was likely the result of genetic load. Our findings shed light on the mechanisms underlying the reduced growth performance of range-edge populations.

The evolutionary history of rice azaleas (Rhododendron tschonoskii alliance) involved niche evolution to a montane environment.

PREMISE: The formation of isolated montane geography on islands promotes evolution, speciation, and then radiation if there are ecological changes. Thus, investigating evolutionary histories of montane species and associated ecological changes may help efforts to understand how endemism formed in islands' montane floras. To explore this process, we investigated the evolutionary history of the Rhododendron tschonoskii alliance, which grows in montane environments of the Japanese archipelago and the Korean Peninsula. METHODS: We studied the five species in the R. tschonoskii alliance and 30 outgroup species, using genome-wide single-nucleotide polymorphisms and cpDNA sequences, in association with environmental analyses. RESULTS: The monophyletic R. tschonoskii alliance diverged since the late Miocene. Species in the alliance currently inhabit a cold climatic niche that is largely different from that of the outgroup species. We observed clear genetic and niche differentiations between the taxa of the alliance. CONCLUSIONS: The association of the alliance's evolution with the formation of cooler climates on mountains indicates that it was driven by global cooling since the mid-Miocene and by rapid uplift of mountains since the Pliocene. The combination of geographic and climatic isolation promoted high genetic differentiation between taxa, which has been maintained by climatic oscillations since the Quaternary.

5S-IGS rDNA in wind-pollinated trees (Fagus L.) encapsulates 55 million years of reticulate evolution and hybrid origins of modern species.

Standard models of plant speciation assume strictly dichotomous genealogies in which a species, the ancestor, is replaced by two offspring species. The reality in wind-pollinated trees with long evolutionary histories is more complex: species evolve from other species through isolation when genetic drift exceeds gene flow; lineage mixing can give rise to new species (hybrid taxa such as nothospecies and allopolyploids). The multi-copy, potentially multi-locus 5S rDNA is one of few gene regions conserving signal from dichotomous and reticulate evolutionary processes down to the level of intra-genomic recombination. Therefore, it can provide unique insights into the dynamic speciation processes of lineages that diversified tens of millions of years ago. Here, we provide the first high-throughput sequencing (HTS) of the 5S intergenic spacers (5S-IGS) for a lineage of wind-pollinated subtropical to temperate trees, the Fagus crenata - F. sylvatica s.l. lineage, and its distant relative F. japonica. The observed 4963 unique 5S-IGS variants reflect a complex history of hybrid origins, lineage sorting, mixing via secondary gene flow, and intra-genomic competition between two or more paralogous-homoeologous 5S rDNA lineages. We show that modern species are genetic mosaics and represent a striking case of ongoing reticulate evolution during the past 55 million years.

Assessments of fine-scale spatial patterns of SNPs in an old-growth beech forest.

The spatial patterns of non-neutral genetic variations at fine spatial scales and their possible associations with microenvironments have not been well-documented for tree populations. Based on 25-32 SNP markers, we examine whether non-neutral SNPs and their associations with microenvironments can be detected in FcMYB1603, a gene homologous to that encoding a protein induced by drought stress in Arabidopsis thaliana for the 166 adult trees in a 1-ha plot in a mature population of Fagus crenata. In the 83 individuals of a younger cohort of below canopy trees, the nonsynonymous SNP at locus FcMYB1603_684 exhibited a spatial signature representing a departure from the expected spatial patterns of neutral genetic variation. Evaluations of non-neutrality for this locus were robust against the potential risks of false positives due to the low number of SNP loci, a low criterion set for minor allele frequency, and any edge effect on the trees' spatial structure. An older cohort exhibited no signal of the existence of non-neutral genetic variation, suggesting that temporal fluctuation in the microenvironmental conditions on the forest floor may have exposed different cohorts to different magnitudes of selection pressure. Although genotypes of the locus showed a spatial association with a microenvironmental variable potentially related to soil moisture, the present study was subject to a limitation due to the generally low polymorphism of nonsynonymous loci within the single plot, which suggests that it will be important to replicate the study design in order to carry out research on fine-scale non-neutral genetic variations.

Correction to: Population genetic structure and demography of Magnolia kobus: variety borealis is not supported genetically.

The original article was updated.

Spatial variation in bird pollination and its mitigating effects on the genetic diversity of pollen pools accepted by Camellia japonica trees within a population at a landscape level.

Bird pollination can vary spatially in response to spatial fluctuations in flowering even within plant populations. In this study, we examined the hypothesis that the spatial variation in bird pollination may induce mitigating effects, which maintains or increases genetic diversity of pollen pools at local sites with low flowering densities. To test this hypothesis, we analyzed the landscape-level genetic effects within a population of Camellia japonica on the pollen pools accepted by individuals in two reproductive years by using genotypes at eight microsatellite loci of 1323 seeds from 19 seed parents. Regression analyses using the quadratic models of correlated paternity between pollen pools against spatial distances between the seed-parent pairs revealed not only local pollination but also some amount of long-distance pollen dispersal. The genetic diversity of pollen pools accepted by seed parents tended to be negatively related to the densities of flowering individuals near the seed parents during winter (when the effective pollination of C. japonica is mediated mostly by Zosterops japonica). We show that the low density of flowering individuals may induce the expansion of the foraging areas of Z. japonica and consequently increase the genetic diversity of pollen pools. This spatial variation in bird pollination may induce the mitigating effects on the C. japonica population. The comparisons between the two study years indicate that the overall pattern of bird pollination and the genetic effects described here, including the mitigating effects, may be stable over time.

Population genetic structure and demography of Magnolia kobus: variety borealis is not supported genetically.

Species delimitations by morphological and by genetic markers are not always congruent. Magnolia kobus consists of two morphologically different varieties, kobus and borealis. The latter variety is characterized by larger leaves than the former. For the conservation of M. kobus genetic resources in natural forests, the relationships between morphological and genetic variation should be clarified. We investigated variations in nuclear microsatellites, chloroplast DNA (cpDNA) sequences and leaf morphological traits in 23 populations of M. kobus over the range of species. Two genetically divergent lineages, northern and southern were detected and their geographical boundary was estimated to be at 39°N. The northern lineage consisted of two genetic clusters and a single cpDNA haplotype, while the southern one had multiple genetic clusters and cpDNA haplotypes. The northern lineage showed significantly lower genetic diversity than the southern. Approximate Bayesian computation indicated that the northern and southern lineages had experienced, respectively, population expansion and long-term stable population size. The divergence time between the two lineages was estimated to be 565,000 years ago and no signature of migration between the two lineages after divergence was detected. Ecological niche modeling showed that the potential distribution area in northern Japan at the last glacial maximum was very small. It is thus considered that the two lineages have experienced different population histories over several glacial-inter-glacial cycles. Individuals of populations in the central to northern part of Honshu on the Sea of Japan side and in Hokkaido had large leaf width and area. These leaf characteristics corresponded with those of variety borealis. However, the delimitation of the northern and southern lineages detected by genetic markers (39°N) was not congruent with that detected by leaf morphologies (36°N). It is therefore suggested that variety borealis is not supported genetically and the northern and southern lineages should be considered separately when identifying conservation units based not on morphology but on genetic markers.

The complete chloroplast genome of Fagus crenata (subgenus Fagus) and comparison with F. engleriana (subgenus Engleriana).

This study reports the whole chloroplast genome of Fagus crenata (subgenus Fagus), a foundation tree species of Japanese temperate forests. The genome has a total of 158,227 bp containing 111 genes, including 76 protein-coding genes, 31 tRNA genes and 4 ribosomal RNA genes. Comparison with the only other published Fagus chloroplast genome, F. engeleriana (subgenus Engleriana) shows that the genomes are relatively conserved with no inversions or rearrangements observed while the proportion of nucleotide sites differing between the two species was equal to 0.0018. The six most variable regions were, in increasing order of variability, psbK-psbI, trnG-psbfM, rpl32, trnV, ndhI-ndh and ndhD-psaC. These highly variable chloroplast regions in addition to 160 chloroplast microsatellites identified (of which 46 were variable between the two species) will provide useful genetic resources for studies of the inter- and intra-specific genetic structure and diversity of this important northern hemisphere tree genus.

Sex change in the subdioecious shrub Eurya japonica (Pentaphylacaceae).

Sex change affects the sex ratios of plant populations and may play an essential role in the evolutionary shift of sexual systems. Sex change can be a strategy for increasing fitness over the lifetime of a plant, and plant size, environmental factors, and growth rate may affect sex change. We described frequent, repeated sex changes following various patterns in a subdioecious Eurya japonica population over five successive years. Of the individuals, 27.5% changed their sex at least once, and these changes were unidirectional or bidirectional. The sex ratio (females/males/all hermaphrodite types) did not fluctuate over the 5 years. In our study plots, although the current sex ratio among the sexes appears to be stable, the change in sex ratio may be slowly progressing toward increasing females and decreasing males. Sex was more likely to change with higher growth rates and more exposure to light throughout the year. Among individuals that changed sex, those that were less exposed to light in the leafy season and had less diameter growth tended to shift from hermaphrodite to a single sex. Therefore, sex change in E. japonica seemed to be explained by a response to the internal physiological condition of an individual mediated by intrinsic and abiotic environmental factors.

Reduced incompatibility in the production of second generation hybrids between two Magnolia species revealed by Bayesian gene dispersal modeling.

PREMISE OF THE STUDY: Hybrid zones are areas where gene flow between related species is currently occurring, so information on the compatibility between related species and their hybrids is essential for predicting the dynamics of such zones generated by introgressive hybridization. In this study, we quantified the compatibility among Magnolia stellata, M. salicifolia, and their hybrids in a hybrid zone using gene dispersal modeling. METHODS: After determining the genealogical classes of adult trees in the hybrid zone, the paternity of 574 open-pollinated seeds from 37 known maternal trees was analyzed with microsatellite markers. A neighborhood-based Bayesian gene dispersal model developed by us for estimating compatibility was then applied to the paternity data. KEY RESULTS: When M. stellata or M. salicifolia were mothers, interspecific mating to produce F1 hybrids yielded significant incompatibility, but backcrossing with F1 hybrids did not. Furthermore, when F1 hybrids became mothers, no significant incompatibility resulted from backcrossing to parental species or intra-F1 mating to produce F2 hybrids. The estimated proportion of F1 hybrids in the outcrossed seeds (1.7%) in the hybrid zone was much lower than that in the adult trees (14.0%). CONCLUSIONS: While it is difficult to obtain F1 hybrids, their low incompatibility makes it easy to produce advanced generation hybrids, once they have been successfully obtained. Although the production of F1 seeds is rare, heterosis and/or weak selection pressure in an empty niche between the parental species' niches may have contributed to the increased proportion of adult F1 hybrids in the hybrid zone.

Population demographic history of a temperate shrub, Rhododendron weyrichii (Ericaceae), on continental islands of Japan and South Korea.

Continental islands provide opportunities for testing the effects of isolation and migration on genetic variation in plant populations. In characteristic of continental islands is that the geographic connections between these islands, which are currently distinguished by seaways, have experienced fluctuations caused by sea-level changes due to climate oscillations during the Quaternary. Plant populations on the islands have migrated between these islands via the exposed seafloors or been isolated. Here, we examined the demographic history of a temperate shrub, Rhododendron weyrichii, which is distributed in the southwestern parts of the Japanese archipelago and on an island of South Korea, using statistical phylogeographic approaches based on the DNA sequences of two chloroplast and eight nuclear loci in samples analyzed from 18 populations on eight continental islands, and palaeodistribution modeling. Time estimates for four island populations indicate that the durations of vicariance history are different between these populations, and these events have continued since the last glacial or may have predated the last glacial. The constancy or expansion of population sizes on the Japanese islands, and in contrast a bottleneck in population size on the Korean island Jeju, suggests that these islands may have provided different conditions for sustaining populations. The result of palaeodistribution modeling indicates that the longitudinal range of the species as a whole has not changed greatly since the last glacial maximum. These results indicate that exposed seafloors during the glacial period formed both effective and ineffective migration corridors. These findings may shed light on the effects of seafloor exposure on the migration of plants distributed across continental islands.

High male fertility in males of a subdioecious shrub in hand-pollinated crosses.

Female reproductive success in females versus hermaphrodites has been well documented. However, documenting a potential advantage in male fertility of male versus hermaphrodite individuals in subdioecious species is also essential for understanding the evolutionary pathway toward dioecy from hermaphroditism via gynodioecy. Siring success in terms of fruit set, fruit mass, number of seeds and mean seed mass was compared by hand-pollinated crosses in the subdioecious shrub Eurya japonica The pollen was from male and hermaphrodite individuals, and the pollen recipients were females and hermaphrodites. Seed quality was also evaluated in terms of seed germination rate, seed germination day and seedling survival. Overall, pollen from males sired more fruits of larger size and more seeds than did pollen from hermaphrodites. The male advantage was observed when pollen recipients were females, whereas no effect was found in hermaphrodite recipients. Pollen from males also produced better quality seeds with higher germination rate and sooner germination day. Although these results could also be explained by a higher pollen load for crosses with male pollen donors, we took care to saturate the stigma regardless of the pollen donor. Therefore, these results suggest that male individuals of E. japonica have advantages in male fertility in terms of both quantity and quality. Our previous studies indicated that females exhibit higher female reproductive success compared with hermaphrodites. Thus, both the female and male functions of hermaphrodites are outperformed by females and males, respectively, raising the possibility that the subdioecious E. japonica at this study site is entering the transitional phase to dioecy along the gynodioecy-dioecy pathway.

Pollen flow and effects of population structure on selfing rates and female and male reproductive success in fragmented Magnolia stellata populations.

BACKGROUND: Fragmentation of plant populations may affect mating patterns and female and male reproductive success. To improve understanding of fragmentation effects on plant reproduction, we investigated the pollen flow patterns in six adjacent local populations of Magnolia stellata, an insect-pollinated, threatened tree species in Japan, and assessed effects of maternal plant (genet) size, local genet density, population size and neighboring population size on female reproductive success (seed production rates), and effects of mating distance, paternal genet size, population size and separation of populations on male reproductive success. RESULTS: The seed production rate, i.e. the proportion of ovules that successfully turned into seeds, varied between 1.0 and 6.5%, and increased with increasing population size and neighboring population size, and with decreasing maternal genet size and local genet density. The selfing rate varied between 3.6 and 28.9%, and increased with increasing maternal genet size and with declining local genet density. Male reproductive success increased with increasing paternal genet size, and decreased with increasing mating distance and separation of population. Pollen flow between the populations was low (6.1%) and highly leptocurtic. CONCLUSIONS: Our results indicate that habitat fragmentation, separation and reduced size of populations, affected mating patterns and reproductive success of M. stellata. Local competition for pollinators and plant display size were likely to alter the reproductive success.

Patterns of pollen flow in a dense population of the insect-pollinated canopy tree species Castanopsis sieboldii.

Insect pollinations of tree species with high-density populations have rarely been studied. Since the density of adults can affect effective pollen dispersal, short-distance pollination, even by insects, may frequently occur in high-density populations. To test this prediction, we investigated pollination patterns in a high-density population of the insect-pollinated canopy tree species Castanopsis sieboldii by paternity analysis using genotypes at 8 microsatellite loci of 145 adult trees and 439 seeds from 11 seed parents in a 4-ha plot. We then explored their genetic effects on the population by calculating other population genetics parameters. Although C. sieboldii has high potential for long-distance dispersal of pollen (as indicated by a fat-tailed dispersal kernel), the cumulative pollination at the local scale was spatially limited and strongly dependent on the distance between parents due to the high density of adults. Genetic diversity estimates for pollen pools accepted by each seed parent converged on a maximum as the effective number of pollen parents increased. The genetic diversity of pollen pool bulked over all the seed parents from inside the plot did not differ from that of the total pollen pools. Therefore, although pollen flow from distant pollen parents may help to maintain the genetic diversity of offspring, pollen parents neighboring seed parents may be the main contributors to the genetic diversity of the offspring at the seed stage.

Genetic structure of an endemic Japanese conifer, Sciadopitys verticillata (Sciadopityaceae), by using microsatellite markers.

We investigated the genetic structure of Sciadopitys verticillata, an endemic conifer in Japan, using 11 microsatellite markers. Average expected heterozygosity varied from 0.282 to 0.450, with between 2 and 17 alleles per locus, and allelic richness varied from 2.14 to 2.94 within a population. Wilcoxon signed-rank tests suggested that there was a recent bottleneck in the southwestern populations in the Kinki district and in isolated northeastern populations in the Tohoku district. Bayesian clustering analysis showed that the populations of the Chubu, Kinki, and Shikoku districts were admixed from several clusters. Archaeological data suggest past anthropogenic exploitation for building and coffin material.

Interpopulation variation in mating system and late-stage inbreeding depression in Magnolia stellata.

Inbreeding has the potential to cause evolutionary changes in populations, although these changes are likely to drive populations to extinction through inbreeding depression and reductions in genetic diversity. We investigated the mating system and late-stage inbreeding depression (delta) in 10 populations of Magnolia stellata using nine microsatellite markers and evaluated the effects of population size and the degree of population isolation through inbreeding and inbreeding depression on the persistence of populations. The outcrossing rates were very similar (approximately 0.7) among populations, but the correlations of paternity, fractions of biparental inbreeding and inbreeding coefficients at the seed stage (F(S)) varied among populations, suggesting that the level of outcrossing was similar among populations, while the quality of it was not. A significant negative correlation was detected between F(S) and population size. The average value of delta was 0.709, and the values in six of the 10 populations were significant. The values of delta differed among populations, although clear relationships with population size and the degree of population isolation were not detected. However, in one population, which was very small and located in the edge of the species' range, we obtained a very low value of delta (-0.096), which may be indicative of purging or the fixation of deleterious alleles. Existing M. stellata populations that are small (and thus might be expected to have higher frequencies of inbreeding) and have large values of delta may be in danger of declining, even if the populations are located within the central region of the species' range.

Genetic divergence in nuclear genomes between populations of Fagus crenata along the Japan Sea and Pacific sides of Japan.

Genetic diversity and structure in Fagus crenata were studied by analyzing 14 nuclear microsatellite loci in 23 populations distributed throughout the species' range. Although population differentiation was very low (F (ST) = 0.027; R (ST) = 0.041), both neighbor-joining tree and Bayesian clustering analyses provided clear evidence of genetic divergence between populations along the Japan Sea (Japan Sea lineage) and Pacific (Pacific lineage) sides of Japan, indicating that physical barriers to migration and gene flow, notably the mountain ranges separating the populations along the Japan Sea and Pacific sides, have promoted genetic divergence between these populations. The two lineages of the nuclear genome are generally consistent with those of the chloroplast genome detected in a previous study, with several discrepancies between the two genomes. Within-population genetic diversity was generally very high (average H (E) = 0.839), but decreased in a clinal fashion from southwest to northeast, largely among populations of the Japan Sea lineage. This geographical gradient may have resulted from the late-glacial and postglacial recolonization to the northeast, which led to a loss of within-population genetic diversity due to cumulative founder effects.

Eighteen microsatellite loci in Salix arbutifolia (Salicaceae) and cross-species amplification in Salix and Populus species.

Salix arbutifolia is a riparian dioecious tree species that is of conservation concern in Japan because of its highly restricted distribution. Eighteen polymorphic loci of dinucleotide microsatellites were isolated and characterized. Among these, estimates of the expected heterozygosity ranged from 0.350 to 0.879. Cross-species amplification was successful at 9-13 loci among six Salix species and at three loci in one Populus species.

Population differentiation and gene flow within a metapopulation of a threatened tree, Magnolia stellata (Magnoliaceae).

We examined genetic differentiation among eight local populations of a metapopulation of Magnolia stellata using 10 nuclear and three chloroplast microsatellite (nSSR and cpSSR) markers and evaluated the influence of historical gene flow on population differentiation. The coefficient of genetic differentiation among populations for nSSR (F(ST) = 0.053) was less than half that for cpSSR (0.137). An isolation-by-distance pattern was detected for nSSRs, but not cpSSRs. These results suggest that pollen flow, as well as seed dispersal, has significantly reduced genetic differentiation among populations. We also examined patterns of contemporary pollen flow by paternity analysis of seeds from nine seed parents in one of the populations using the nSSR markers and found it to be greatly restricted by the distance between parents. Although most pollen flow occurred within the population, pollen flow from outside the population accounted for 2.5% of the total. When historical and contemporary pollen flows among populations were compared, the levels of pollen flow seem to have declined recently. We conclude that to conserve M. stellata, it is important to preserve the whole population by maintaining its metapopulation structure and the gene flow among its populations.

Effects of pollen shortage and self-pollination on seed production of an endangered tree, Magnolia stellata.

BACKGROUND AND AIMS: Pollen limitation is a significant determinant of seed production, and can result from both insufficient pollen quantity (pollen shortage) and quality (mainly relating to self-pollination). For animal-pollinated tree species with large floral displays, pollen limitation may be determined by a balance between increased pollen quantity due to increased attractiveness for pollinators, countered by increased self-pollination due to increased geitonogamy. The contributions of pollen shortage and self-pollination on seed production were quantitatively examined in the natural pollination of an insect-pollinated, dichogamous, endangered tree, Magnolia stellata, which has a large, showy floral display. METHODS: Manual self- and cross-pollinations were conducted to determine the effects of selfing on seed production. The outcrossing rate was measured using microsatellite analyses of open-pollinated seeds, and the embryo mortality rate caused by self-pollination was indirectly estimated. The frequency of ovule mortality due to pollen shortage was also inferred using the embryo mortality and ovule survival rates from natural pollination. KEY RESULTS: The average fruit set, seed set per fruit, and ovule survival rate per tree from hand cross-pollination were 1.37, 3.15, and 3.34 times higher than those from hand self-pollination, respectively, indicating that self-pollination causes inbreeding depression for fruit and seed set. The multilocus-outcrossing rate (t(m)) was intermediate, 0.632, and the primary selfing rate was 0.657. This indicates that frequent geitonogamous selfing occurs. The ovule mortality rate due to pollen shortage and the embryo mortality rate due to self-pollination were estimated to be 80.8 % and 45.9 %, respectively. CONCLUSIONS: It is concluded that seed production of M. stellata is strongly limited by both pollen shortage and self-pollination. Inefficient beetle-pollination and the automimicry system via asynchronous flowering might be responsible for the high level of pollen shortage and frequent geitonogamy. This is despite a large, showy floral display and the dichogamous system of the species.