Evidence for cryptic northern refugia in the last glacial period in Cryptomeria japonica

Genetic consequences of Last Glacial-Holocene changes in snowfall regime in Arnica mallotopus populations: A plant confined to heavy-snow areas of Japan

PREMISE

Snow is an important environmental factor affecting plant distribution. Past changes in snowfall regimes may have controlled the demographies of snow-dependent plants. However, our knowledge of changes in the distribution and demographies of such plants is limited because of the lack of fossil records.

METHODS

Population genetic and landscape genetic analyses were used to investigate the response of population dynamics of Arnica mallotopus (Asteraceae)-a plant confined to heavy-snow areas of Japan-to changes in snowfall regimes from the Last Glacial Period to the Holocene.

RESULTS

The population genetic analysis suggested that the four geographic lineages diverged during the Last Glacial Period. The interaction between reduced snowfall and lower temperatures during this period likely triggered population isolation in separate refugia. Subpopulation differentiation in the northern group was lower than in the southern group. Our ecological niche model predicted that the current distribution was patchy in the southern region; that is, the populations were isolated by topologically flat and climatically unsuitable lowlands. The landscape genetic analysis suggested that areas with little snowfall acted as barriers to the Holocene expansion of species distribution and continued limiting gene flow between local populations.

CONCLUSIONS

These findings indicate that postglacial population responses vary among regions and are controlled by environmental and geographic factors. Thus, changes in snowfall regime played a major role in shaping the distribution and genetic structure of the snow-dependent plant.

Approximate Bayesian computation and ecological niche models elucidate the demographic history and current fragmented population distribution of a Korean endemic shrub

Climatic fluctuations and geological events since the LGM are believed to have significantly impacted the population size, distribution, and mobility of many species that we observe today. In this paper, we determined the processes driving the phylogeographic structure of the Korean endemic white forsythia by combining the use of genome-wide SNPs and predicting paleoclimatic habitats during the LGM (21 kya), Early Holocene (10 kya), Mid-Holocene (6 kya), and Late Holocene (3 kya). Using a maximum of 1897 SNPs retrieved from 124 samples across nine wild populations, five environmental predictors, and the species' natural occurrence records, we aimed to infer the species' demographic history and reconstruct its possible paleodistributions with the use of approximate Bayesian computation and ecological niche models, respectively. Under this integrated framework, we found strong evidence for patterns of range shift and expansion, and population divergence events from the onset of the Holocene, resulting in the formation of its five distinct genetic units. The most highly supported model inferred that after the split of an ancestral population into the southern group and a larger central metapopulation lineage, the latter gave rise to the eastern and northern clusters, before finally dividing into two sub-central groups. While the use of molecular data allowed us to identify and refine the (phylo)genetic relationships of the species' lineages and populations, the use of ecological data helped us infer a past LGM refugium and the directions of post-glacial range dynamics. The time frames of these demographic events were shown to be congruent with climatic and geological events that affected the central Korean Peninsula during these periods. These findings gave us a better understanding of the consequences of past spatiotemporal factors that may have resulted in the current fragmented population distribution of this endangered plant.

Genetic differentiation and evolution of broad-leaved evergreen shrub and tree varieties of Daphniphyllum macropodum (Daphniphyllaceae)

Tree form evolution is an important ecological specialization for woody species, but its evolutionary process with adaptation is poorly understood, especially on the microevolutionary scale. Daphniphyllum macropodum comprises two varieties: a tree variety growing in a warm temperate climate with light snowfall and a shrub variety growing in a cool temperate climate with heavy snowfall in Japan. Chloroplast DNA variations and genome-wide single-nucleotide polymorphisms across D. macropodum populations and D. teijsmannii as an outgroup were used to reveal the evolutionary process of the shrub variety. Population genetic analysis indicated that the two varieties diverged but were weakly differentiated. Approximate Bayesian computation analysis supported a scenario that assumed migration between the tree variety and the southern populations of the shrub variety. We found migration between the two varieties where the distributions of the two varieties are in contact, and it is concordant with higher tree height in the southern populations of the shrub variety than the northern populations. The genetic divergence between the two varieties was associated with snowfall. The heavy snowfall climate is considered to have developed since the middle Quaternary in this region. The estimated divergence time between the two varieties suggests that the evolution of the two varieties may be concordant with such paleoclimatic change.

CLIMATIC PREREQUISITES FOR NATURALIZATION OF <i>MAGNOLIA SIEBOLDII </i>S.L. IN RUSSIA

The article presents the results of the study of the invasive potential of K. Koch on the territory of the Botanical Garden - Institute of the Far Eastern Branch of the Russian Academy of Sciences (BSI FEB RAS). Since 2015, we have observed the emergence of self-seedlings of M. sieboldii in natural ecosystems. This invasive population has a heterogeneous structure and morphometric parameters. Our results show the dependence of the mass appearance of seedlings on phenological shifts caused by modern climatic changes. We used a method of modeling potential habitats to describe changes in the boundaries of the natural range for the period from the last glacial maximum to the future in the short term (2070). A tendency for expansion of M. siebo ldii range in the northeastern direction, including potential to distribution of this species over the territory of the Russian Far East has been established.

Phylogenetic origin of two Japanese Torreya taxa found in two regions with strongly contrasting snow depth

The Japanese archipelago exhibits a notable difference in snow depth in winter, deep snow on the Sea of Japan side and low snow cover on the Pacific Ocean side. This contrasting pattern has shaped the distribution of infraspecific taxon pairs in a range of woody plants, with taxa found on the Sea of Japan side typically exhibiting a stunted shrub form with multiple decumbent stems. The phylogenetic origin of these taxon pairs is unknown, i.e., whether the two taxa diverged from the same species or if they have different origins. This study aimed to reveal the phylogenetic origin of two varieties of Torreya nucifera (Taxaceae); var. nucifera is a tree found on the Pacific Ocean side, whereas var. radicans is a shrub found on the Sea of Japan side. We examined the phylogenetic relationships of the two varieties and worldwide Torreya taxa using whole chloroplast genomes, chloroplast DNA fragments, and the nuclear ribosomal internal transcribed spacer (ITS). The whole chloroplast genome phylogeny indicated that T. nucifera var. radicans was a sister taxon to Chinese T. grandis, rather than to var. nucifera. In contrast, the nuclear ITS phylogeny indicated that while several haplotypes of T. nucifera var. radicans were closely related to T. grandis, most haplotypes of T. nucifera var. radicans formed a single clade with those of var. nucifera. This implies that the homogenization of the ITS has occurred between the two taxa, while taxon-specific chloroplast DNA haplotypes were retained. These discordant phylogenies suggested that the two taxa have different phylogenetic origins, but have an intricate evolutionary history, involving inter-taxa hybridization and gene flow, possibly when their distributions were confined to sympatric refugia. Given the genetic evidence and distinct difference in growth form, we propose that T. nucifera var. radicans should be taxonomically treated as a distinct species, T. fruticosa.

Genetic Distinctiveness but Low Diversity Characterizes Rear-Edge Thuja standishii (Gordon) Carr. (Cupressaceae) Populations in Southwest Japan

Rear-edge populations are of significant scientific interest because they can contain allelic variation not found in core-range populations. However, such populations can differ in their level of genetic diversity and divergence reflecting variation in life-history traits, demographic histories and human impacts. Using 13 EST-microsatellites, we investigated the genetic diversity and differentiation of rear-edge populations of the Japanese endemic conifer Thuja standishii (Gordon) Carr. in southwest Japan from the core-range in northeast Japan. Range-wide genetic differentiation was moderate (Fst = 0.087), with northeast populations weakly differentiated (Fst = 0.047), but harboring high genetic diversity (average population-level Ar = 4.76 and Ho = 0.59). In contrast, rear-edge populations were genetically diverged (Fst = 0.168), but contained few unique alleles with lower genetic diversity (Ar = 3.73, Ho = 0.49). The divergence between rear-edge populations exceeding levels observed in the core-range and results from ABC analysis and species distribution modelling suggest that these populations are most likely relicts of the Last Glacial Maximum. However, despite long term persistence, low effective population size, low migration between populations and genetic drift have worked to promote the genetic differentiation of southwest Japan populations of T. standishii without the accumulation of unique alleles.

Factors Affecting the Number of Pollen Grains per Male Strobilus in Japanese Cedar (Cryptomeria japonica)

Japanese cedar (Cryptomeria japonica) is the most important timber species in Japan; however, its pollen is the primary cause of pollinosis in Japan. The total number of pollen grains produced by a single tree is determined by the number of male strobili (male flowers) and the number of pollen grains per male strobilus. While the number of male strobili is a visible and well-investigated trait, little is known about the number of pollen grains per male strobilus. We hypothesized that genetic and environmental factors affect the pollen number per male strobilus and explored the factors that affect pollen production and genetic variation among clones. We counted pollen numbers of 523 male strobili from 26 clones using a cell counter method that we recently developed. Piecewise Structural Equation Modeling (pSEM) revealed that the pollen number is mostly affected by genetic variation, male strobilus weight, and pollen size. Although we collected samples from locations with different environmental conditions, statistical modeling succeeded in predicting pollen numbers for different clones sampled from branches facing different directions. Comparison of predicted pollen numbers revealed that they varied >3-fold among the 26 clones. The determination of the factors affecting pollen number and a precise evaluation of genetic variation will contribute to breeding strategies to counter pollinosis. Furthermore, the combination of our efficient counting method and statistical modeling will provide a powerful tool not only for Japanese cedar but also for other plant species.

Diversification of terpenoid emissions proposes a geographic structure based on climate and pathogen composition in Japanese cedar

Biogenic volatile organic compounds emitted from plants are important constituents of atmospheric chemistry and play a major role in the resistance of plants against various environmental stresses. However, little is known about how abiotic and biotic environments on a geographic scale relate to diversifications of the emission. Here, we present variations of terpenes stored in and emitted from leaves of a single species in a common garden, using genetically differentiated local populations of Japanese cedar, the most dominant and widely distributed tree species in Japan. Furthermore, we determined the composition of fungal communities in 50 locations, based on the presence or absence of 158 fungal species inhabiting the cedar. The results showed that terpenoids, especially those that are emitted, were highly diversified and geographically structured among the 12 populations. The total amount of stored terpenes was negatively affected by warm and less-snow climates. On the other hand, variations in some emitted terpenoid species among the populations were correlated to antagonistic fungal species inhabiting the Japanese cedar. We propose that the diversification of composition and amount of stored and emitted terpenoids in the tree species is not only structured by climate, but also antagonistic fungal communities through biological interactions.

Genetic Diversity and Structure of Japanese Endemic Genus Thujopsis (Cupressaceae) Using EST-SSR Markers

The genus Thujopsis (Cupressaceae) comprises monoecious coniferous trees endemic to Japan. This genus includes two varieties: Thujopsis dolabrata (L.f.) Siebold et Zucc. var. dolabrata (southern variety, Td) and Thujopsis dolabrata (L.f.) Siebold et Zucc. var. hondae Makino (northern variety, Th). The aim of this study is to understand the phylogeographic and genetic population relationships of the genus Thujopsis for the conservation of genetic resources and future breeding. A total of 609 trees from 22 populations were sampled, including six populations from the Td distribution range and 16 populations from the Th distribution range. The genotyping results for 19 expressed sequence tag (EST)-based simple sequence repeat (SSR) markers, followed by a structure analysis, neighbor-joining tree creation, an analysis of molecular variance (AMOVA), and hierarchical F statistics, supported the existence of two genetic clusters related to the distribution regions of the Td and Th varieties. The two variants, Td and Th, could be defined by their provenance, in spite of the ambiguous morphological differences between the varieties. The distribution ranges of both variants, which have been defined from their morphology, was confirmed by genetic analysis. The Th populations exhibited relatively uniform genetic diversity, most likely because Th refugia in the glacial period were scattered throughout their current distribution area. On the other hand, there was a tendency for Td’s genetic diversity to decrease from central to southern Honshu island. Notably, the structure analysis and neighbor-joining tree suggest the hybridization of the two varieties in the contact zone. More detailed studies of the genetic structure of Td are required in future analyses.

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.

Inferring the demographic history of Japanese cedar, Cryptomeria japonica, using amplicon sequencing

The evolution of a species depends on multiple forces, such as demography and natural selection. To understand the trajectory and driving forces of evolution of a target species, it is first necessary to uncover that species' population history, such as past and present population sizes, subdivision and gene flow, by using appropriate genetic markers. Cryptomeria japonica is a long-lived monoecious conifer species that is distributed in Japan. There are two main lines (omote-sugi and ura-sugi), which are distinguished by apparent differences in morphological traits that may have contributed to their local adaptation. The evolution of these morphological traits seems to be related to past climatic changes in East Asia, but no precise estimate is available for the divergence time of these two lines and the subsequent population dynamics in this species. Here, we analyzed the nucleotide variations at 120 nuclear genes in 94 individuals by using amplicon sequencing in combination with high-throughput sequencing technologies. Our analysis indicated that the population on Yakushima Island, the southern distribution limit of C. japonica in Japan, diverged from the other populations 0.85 million years ago (MYA). The divergence time of the other populations on mainland Japan was estimated to be 0.32 MYA suggesting that the divergence of omote-sugi and ura-sugi might have occurred before the last glacial maximum. Although we found modest levels of gene flow between the present populations, the long-term isolation and environmental heterogeneity caused by climatic changes might have contributed to the differentiation of the lines and their local adaptation.

Potential of Genome-Wide Studies in Unrelated Plus Trees of a Coniferous Species, Cryptomeria japonica (Japanese Cedar)

A genome-wide association study (GWAS) was conducted on more than 30,000 single nucleotide polymorphisms (SNPs) in unrelated first-generation plus tree genotypes from three populations of Japanese cedar Cryptomeria japonica D. Don with genomic prediction for traits of growth, wood properties and male fecundity. Among the assessed populations, genetic characteristics including the extent of linkage disequilibrium (LD) and genetic structure differed and these differences are considered to be due to differences in genetic background. Through population-independent GWAS, several significant SNPs found close to the regions associated with each of these traits and shared in common across the populations were identified. The accuracies of genomic predictions were dependent on the traits and populations and reflected the genetic architecture of traits and genetic characteristics. Prediction accuracies using SNPs selected based on GWAS results were similar to those using all SNPs for several combinations of traits and populations. We discussed the application of genome-wide studies for C. japonica improvement.

Genetic and chemical differentiation characterizes top-geoherb and non-top-geoherb areas in the TCM herb rhubarb

Medicinal herbs of high quality and with significant clinical effects have been designated as top-geoherbs in traditional Chinese medicine (TCM). However, the validity of this concept using genetic markers has not been widely tested. In this study, we investigated the genetic variation within the Rheum palmatum complex (rhubarb), an important herbal remedy in TCM, using a phylogeographic (six chloroplast DNA regions, five nuclear DNA regions, and 14 nuclear microsatellite loci) and a chemical approach (anthraquinone content). Genetic and chemical data identified two distinct groups in the 38 analysed populations from the R. palmatum complex which geographically coincide with the traditional top-geoherb and non-top-geoherb areas of rhubarb. Molecular dating suggests that the two groups diverged in the Quaternary c. 2.0 million years ago, a time of repeated climate changes and uplift of the Qinghai-Tibetan Plateau. Our results show that the ancient TCM concept of top-geoherb and non-top-geoherb areas corresponds to genetically and chemically differentiated groups in rhubarb.

Response of spatial vegetation distribution in China to climate changes since the Last Glacial Maximum (LGM)

Analyzing how climate change affects vegetation distribution is one of the central issues of global change ecology as this has important implications for the carbon budget of terrestrial vegetation. Mapping vegetation distribution under historical climate scenarios is essential for understanding the response of vegetation distribution to future climatic changes. The reconstructions of palaeovegetation based on pollen data provide a useful method to understand the relationship between climate and vegetation distribution. However, this method is limited in time and space. Here, using species distribution model (SDM) approaches, we explored the climatic determinants of contemporary vegetation distribution and reconstructed the distribution of Chinese vegetation during the Last Glacial Maximum (LGM, 18,000 14C yr BP) and Middle-Holocene (MH, 6000 14C yr BP). The dynamics of vegetation distribution since the LGM reconstructed by SDMs were largely consistent with those based on pollen data, suggesting that the SDM approach is a useful tool for studying historical vegetation dynamics and its response to climate change across time and space. Comparison between the modeled contemporary potential natural vegetation distribution and the observed contemporary distribution suggests that temperate deciduous forests, subtropical evergreen broadleaf forests, temperate deciduous shrublands and temperate steppe have low range fillings and are strongly influenced by human activities. In general, the Tibetan Plateau, North and Northeast China, and the areas near the 30°N in Central and Southeast China appeared to have experienced the highest turnover in vegetation due to climate change from the LGM to the present.

Phylogeography of postglacial range expansion in Juglans mandshurica (Juglandaceae) reveals no evidence of bottleneck, loss of genetic diversity, or isolation by distance in the leading-edge populations

The past studies of postglacial recolonization patterns in high latitude regions have revealed a significant role of dispersal capacity in shaping the genetic diversity and population structure of temperate trees. However, most of these studies have focused on species with long-distance dispersal followed by exponential population growth and were therefore unable to reveal the patterns in the case of a gradual expansion. Here we studied the impacts of postglacial range expansions on the distribution of genetic diversity in the Manchurian walnut (Juglans mandshurica), a common tree of East Asian cool-temperate deciduous forests that apparently lacks long-distance seed dispersal ability. The genetic diversity and structure of 19 natural walnut populations in Northeast China and the Korean Peninsula were examined using 17 nuclear simple sequence repeat (SSR) loci. Potential habitats under current and past climatic conditions were predicted using the ecological niche modelling (ENM) method. Bayesian clustering analysis revealed three groups, which were inferred to have diverged through multiple glacial-interglacial cycles in multiple refugia during the Quaternary Period. ENM estimated a southward range shift at the LGM, but high suitability scores still occurred in the western parts of the Changbai Mountains (Northeast China), the Korean peninsula and the exposed seafloor of the Yellow Sea. In contrast to most other cool-temperate trees co-occurring in the same region, the Manchurian walnut did not show any evidence of a population bottleneck, loss of genetic diversity or isolation by distance during the postglacial expansion. Our study clearly indicates that current northern populations originated from one glacial lineage and recolonization via a gradually advancing front due to the lack of a long-distance seed dispersal mechanism led to no latitudinal decrease in genetic diversity.

Genetic differentiation in spite of high gene flow in the dominant rainforest tree of southeastern Australia, Nothofagus cunninghamii

Nothofagus cunninghamii is a long-lived, wind-pollinated tree species that dominates the cool temperate rainforests of southeastern Australia. The species' distribution is more or less continuous in western Tasmania but is fragmented elsewhere. However, it is unknown whether this fragmentation has affected the species' genetic architecture. Thus, we examined N. cunninghamii using 12 nuclear microsatellites and 633 individuals from 18 populations spanning the species' natural range. Typical of wind-pollinated trees, there was low range-wide genetic structure (FST=0.04) consistent with significant gene flow across most of the species' range. However, gene flow was not high enough to overcome the effects of drift across some disjunctions. Victorian populations (separated from Tasmania by the 240 km wide Bass Strait) formed a genetic group distinct from Tasmanian populations, had lower diversity (mean allelic richness (Ar)=5.4 in Victoria versus 6.9 in Tasmania) and were significantly more differentiated from one another than those in Tasmania (FST=0.045 in Victoria versus 0.012 in Tasmania). Evidence for bottlenecking was found in small populations that were at least 20 km from other populations. Interestingly, we found little divergence in microsatellite markers between the extremes of genetically based morphological and physiological altitudinal clines suggesting adaptive differentiation is strongly driven by selection because it is likely to be occurring in the presence of gene flow. Even though the cool temperate rainforests of Australia are highly relictual, the species is relatively robust to population fragmentation due to high levels of genetic diversity and gene flow, especially in Tasmania.

Genetic differentiation and evolutionary adaptation in Cryptomeria japonica

Local adaptation of plant species is a central issue for survival during global climate change, especially for long-lived forest trees, with their lengthy regeneration time and spatially limited gene flow. Identification of loci and/or genomic regions associated with local adaptation is necessary for knowledge of both evolution and molecular breeding for climate change. Cryptomeria japonica is an important species for forestry in Japan; it has a broad natural distribution and can survive in a range of different environments. The genetic structure of 14 natural populations of this species was investigated using 3930 SNP markers. Populations on the Pacific Ocean side of Japan are clearly different from those on the Japan Sea side, as discussed in previous studies. Structure analysis and population network trees show that peripheral populations, including the most northerly and southerly ones, have unique features. We found that the genetic differentiation coefficient is low, F_ST = 0.05, although it must account for the presence of important genes associated with adaptation to specific environments. In total, 208 outlier loci were detected, of which 43 were associated with environmental variables. Four clumped regions of outlier loci were detected in the genome by linkage analysis. Linkage disequilibrium (LD) was quite high in these clumps of outlier loci, which were found in linkage groups (LGs) 2, 7, 10, and 11, especially between populations of two varieties, and when interchromosomal LD was also detected. The LG7 region is characteristic of the Yakushima population, which is a large, isolated, peripheral population occupying a specific environment resulting from isolation combined with volcanic activity in the region. The detected LD may provide strong evidence for selection between varieties.