Response of stomatal density and size in Betula ermanii to contrasting climate conditions: The contributions of genetic and environmental factors.

As plant distribution and performance are determined by both environmental and genetic factors, clarifying the contribution of these two factors is a key for understanding plant adaptation and predicting their distribution under ongoing global warming. Betula ermanii is an ideal species for such research because of its wide distribution across diverse environments. Stomatal density and size are crucial traits that plants undergo changes in to adapt to different environments as these traits directly influence plant photosynthesis and transpiration. In this study, we conducted a multi-location common garden experiment using B. ermanii to (1) clarify the contribution of both environmental and genetic factors to the variation in stomatal density and size of B. ermanii, (2) demonstrate the differences in the plasticity of stomatal density and size among B. ermanii populations, and (3) understand how stomatal density and size of B. ermanii would respond to increased temperature and changing precipitation patterns. Genetic factors played a more significant role in stomatal size than environmental factors, suggesting that B. ermanii struggles to adjust its stomatal size in response to a changing environment. Our results also revealed a positive correlation between stomatal size plasticity and original habitat suitability, indicating that in B. ermanii populations in harsh environments exhibit lower adaptability to environmental shifts. Although stomatal density and size of B. ermanii showed the significant responses to increased temperature and shifting precipitation patterns, the response ranges of stomatal density and size to the environmental factors varied among populations. Our findings highlighted the interplay between genetic and environmental factors in determining the intraspecific variation in stomatal density and size in B. ermanii. This indicated that certain populations of B. ermanii exhibit limited stomatal plasticity and adaptability, which could directly affect photosynthesis and transpiration, suggesting potential population-specific fitness implications for B. ermanii under future climate change.

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.

Comparing modeling methods of genomic prediction for growth traits of a tropical timber species, Shorea macrophylla.

Introduction: Shorea macrophylla is a commercially important tropical tree species grown for timber and oil. It is amenable to plantation forestry due to its fast initial growth. Genomic selection (GS) has been used in tree breeding studies to shorten long breeding cycles but has not previously been applied to S. macrophylla. Methods: To build genomic prediction models for GS, leaves and growth trait data were collected from a half-sib progeny population of S. macrophylla in Sari Bumi Kusuma forest concession, central Kalimantan, Indonesia. 18037 SNP markers were identified in two ddRAD-seq libraries. Genomic prediction models based on these SNPs were then generated for diameter at breast height and total height in the 7th year from planting (D7 and H7). Results and discussion: These traits were chosen because of their relatively high narrow-sense genomic heritability and because seven years was considered long enough to assess initial growth. Genomic prediction models were built using 6 methods and their derivatives with the full set of identified SNPs and subsets of 48, 96, and 192 SNPs selected based on the results of a genome-wide association study (GWAS). The GBLUP and RKHS methods gave the highest predictive ability for D7 and H7 with the sets of selected SNPs and showed that D7 has an additive genetic architecture while H7 has an epistatic genetic architecture. LightGBM and CNN1D also achieved high predictive abilities for D7 with 48 and 96 selected SNPs, and for H7 with 96 and 192 selected SNPs, showing that gradient boosting decision trees and deep learning can be useful in genomic prediction. Predictive abilities were higher in H7 when smaller number of SNP subsets selected by GWAS p-value was used, However, D7 showed the contrary tendency, which might have originated from the difference in genetic architecture between primary and secondary growth of the species. This study suggests that GS with GWAS-based SNP selection can be used in breeding for non-cultivated tree species to improve initial growth and reduce genotyping costs for next-generation seedlings.

Genetic Structure of the Liriope muscari Polyploid Complex and the Possibility of Its Genetic Disturbance in Japan.

Anthropogenic activities, such as the movement of plants through greening, can result in genetic disturbance that can interfere with local adaptation in wild populations. Although research is underway to prevent genetic disturbance associated with greening, genetic disturbance of intraspecific polyploidy, which is estimated to be present in 24% of vascular plants, has not been well studied. Liriope muscari is a polyploid complex with known diploid (2n = 36), tetraploid (2n = 72), and hexaploid (2n = 108) forms. The plants of this species tolerate dry and hot conditions and are therefore frequently used for greening and gardening. However, the distribution of this polyploid in Japan, its genetic structure, and genetic disturbance are not known. In this study, we investigated the polyploidy distribution and genetic structure in naturally distributed L. muscari in Japan using chloroplast DNA (cpDNA) haplotypes and nuclear DNA (nDNA). Commercially produced individuals were also studied and compared with natural populations to assess any genetic disturbance of the ploidy complex in this species. Chromosome counts, cpDNA, and nDNA results showed three genetically and cytologically distinct groups in Japan: first, a tetraploid group in mainland Japan; second, a hexaploid group in the Ryukyu Islands; and third, a diploid and tetraploid group in the Ryukyu Islands. Significant isolation by distance was also detected within the three groups (p = 0.001). Genetic disturbance due to greening and gardening should be avoided among the three groups. Genetic disturbance can be reduced by using individuals derived from natural populations that are close to the sites used for greening and gardening. For commercially produced individuals, genetic disturbance is unlikely in the Kanto region, an area of high usage, while genetic disturbance is thought possible in the Ryukyu Islands.

Genotype-by-environment interaction and genetic dissection of heartwood color in Cryptomeria japonica based on multiple common gardens and quantitative trait loci mapping.

The heartwood color of a major plantation tree Cryptomeria japonica shows high variability among clones and cultivars, and brighter heartwood has higher value in the usage of non-laminated wood such as in traditional construction, which makes heartwood color an important trait in breeding of this species. However, the genetic basis of the interactions between genetics and the environment on heartwood color has been understudied while these are necessary for effective breeding programs in multiple environmental condition. The objectives of the present study were to evaluate the effects of genetics and environments on heartwood color and how they interact in contrasting environments, and to identify genomic regions controlling heartwood color in C. japonica across multiple environments. Heartwood color in terms of L*a*b* color space and spectral reflectance was measured in common gardens established in three contrasting sites. Quantitative trait loci (QTL) that affect heartwood color were identified using previously constructed highly saturated linkage maps. Results found that heartwood color was largely genetically controlled, and genotype-by-environment interaction explained one-third of the total genetic variance of heartwood color. The effect of the environment was small compared to the effect of genetics, whereas environmental effects largely varied among heartwood color traits. QTL analysis identified a large number of QTLs with small to moderate effects (phenotypic variation explained of 6.6% on average). Some of these QTLs were stably expressed in multiple environments or had pleiotropic effects on heartwood color and moisture content. These results indicated that genetic variation in phenotypic plasticity plays an important role in regulating heartwood color and that the identified QTLs would maximize the breeding efficiency of heartwood color in C. japonica in heterogeneous environments.

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.

Construction of a reference transcriptome for the analysis of male sterility in sugi (Cryptomeria japonica D. Don) focusing on MALE STERILITY 1 (MS1).

Sugi (Cryptomeria japonica D. Don) is an important conifer used for afforestation in Japan. As the genome of this species is 11 Gbps, it is too large to assemble within a short timeframe. Transcriptomics is one approach that can address this deficiency. Here we designed a workflow consisting of three stages to de novo assemble transcriptome using Oases and Trinity. The three transcriptomic stage used were independent assembly, automatic and semi-manual integration, and refinement by filtering out potential contamination. We identified a set of 49,795 cDNA and an equal number of translated proteins. According to the benchmark set by BUSCO, 87.01% of cDNAs identified were complete genes, and 78.47% were complete and single-copy genes. Compared to other full-length cDNA resources collected by Sanger and PacBio sequencers, the extent of the coverage in our dataset was the highest, indicating that these data can be safely used for further studies. When two tissue-specific libraries were compared, there were significant expression differences between male strobili and leaf and bark sets. Moreover, subtle expression difference between male-fertile and sterile libraries were detected. Orthologous genes from other model plants and conifer species were identified. We demonstrated that our transcriptome assembly output (CJ3006NRE) can serve as a reference transcriptome for future functional genomics and evolutionary biology studies.

Identification and genetic diversity analysis of a male-sterile gene (MS1) in Japanese cedar (Cryptomeria japonica D. Don).

Identifying causative genes for a target trait in conifer reproduction is challenging for species lacking whole-genome sequences. In this study, we searched for the male-sterility gene (MS1) in Cryptomeria japonica, aiming to promote marker-assisted selection (MAS) of male-sterile C. japonica to reduce the pollinosis caused by pollen dispersal from artificial C. japonica forests in Japan. We searched for mRNA sequences expressed in male strobili and found the gene CJt020762, coding for a lipid transfer protein containing a 4-bp deletion specific to male-sterile individuals. We also found a 30-bp deletion by sequencing the entire gene of another individual with the ms1. All nine breeding materials with the allele ms1 had either a 4-bp or 30-bp deletion in gene CJt020762, both of which are expected to result in faulty gene transcription and function. Furthermore, the 30-bp deletion was detected from three of five individuals in the Ishinomaki natural forest. From our findings, CJt020762 was considered to be the causative gene of MS1. Thus, by performing MAS using two deletion mutations as a DNA marker, it will be possible to find novel breeding materials of C. japonica with the allele ms1 adapted to the unique environment of each region of the Japanese archipelago.

Scanning RNA-Seq and RAD-Seq approach to develop SNP markers closely linked to MALE STERILITY 1 (MS1) in Cryptomeria japonica D. Don.

Cryptomeria japonica is a major forestry tree species in Japan. Male sterility of the species is caused by a recessive gene, which shows dysfunction of pollen development and results in no dispersed pollen. Because the pollen of C. japonica induces pollinosis, breeding of pollen-free C. japonica is desired. In this study, single nucleotide polymorphism (SNP) markers located at 1.78 and 0.58 cM to a male sterility locus (MS1) were identified from an analysis of RNA-Seq and RAD-Seq, respectively. SNPs closely linked to MS1 were first scanned by a method similar to MutMap, where a type of index was calculated to measure the strength of the linkage between a marker sequence and MS1. Linkage analysis of selected SNP markers confirmed a higher efficiency of the current method to construct a partial map around MS1. Allele-specific PCR primer pair for the most closely linked SNP with MS1 was developed as a codominant marker, and visualization of the PCR products on an agarose gel enabled rapid screening of male sterile C. japonica. The allele-specific primers developed in this study would be useful for establishing the selection of male sterile C. japonica.

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.

TodoFirGene: Developing Transcriptome Resources for Genetic Analysis of Abies sachalinensis.

Todo-matsu (Abies sachalinensis) is one of the most important forestry species in Hokkaido, Japan and is distributed from near sea level to the alpine zone. Due to its wide spatial distribution, the species adapts to its environment, displaying phenotypes of ecological relevance. In order to identify candidate genes under natural selection, we collected the transcriptome from the female and male flower, leaf and inner bark. De novo assembly with 34.7 Gb of sequencing reads produced 158,542 transcripts from 69,618 loci, whose estimated coverage reached 95.6% of conserved eukaryotic genes. Homology searches against publicly available databases identified 134,190 (84.6%) transcripts with at least one hit. In total, 28,944 simple sequence repeats (SSRs) and 80,758 single nucleotide variants (SNVs) were detected from 23,570 (14.9%) and 25,366 (16.0%) transcripts, which were valuable for use in genetic analysis of the species. All the annotations were included in a relational database, TodoFirGene, which provides an interface for various queries and homology search, and can be accessed at http://plantomics.mind.meiji.ac.jp/todomatsu/. This database hosts not only the A. sachalinensis transcriptome but also links to the proteomes of 13 other species, allowing a comparative genomic study of plant species.

Large contribution of clonal reproduction to the distribution of deciduous liana species (Wisteria floribunda) in an old-growth cool temperate forest: evidence from genetic analysis.

Background and Aims: Extensive clonal (vegetative) reproduction in lianas is a common and important life history strategy for regeneration and colonization success. However, few studies have evaluated the contribution of clonal reproduction to stand-level distribution of lianas in their natural habitat using genetic tools. The objectives of the present study were to investigate (1) the contribution of clonal reproduction to the distribution of Wisteria floribunda, (2) the size of clonal patches and (3) how the distribution patterns of W. floribunda clones are affected by micro-topography. Methods: The contribution of clonal reproduction to the distribution of the deciduous liana species W. floribunda was evaluated using genetic analysis across a 6-ha plot of an old-growth temperate forest in Japan and preference in landform between clonal ramets and non-clonal ramets was assessed. Key Results: Of the 391 ramets sampled, clonal reproduction contributed to 71 and 62 % of the total abundance and basal area, respectively, or 57 and 31 % when the largest ramet within a genet was excluded. The large contribution of clonal reproduction to the density and basal area of W. floribunda was consistent with previous observational studies. The largest genet included a patch size of 0.47 ha and ranged over 180 m. Preferred landforms of clonal and non-clonal ramets were significantly different when evaluated by both abundance and basal area. Non-clonal ramets distributed more on lower part of the slope than other landforms in comparison with clonal ramets and trees, possibly reflecting the limitation of clonal growth by stolons. Conclusions: Using genetic analysis, the present study found evidence of a large contribution of clonal reproduction on the distribution of W. floribunda in its natural habitat. The results indicate that clonal reproduction plays an important role not only in the formation of populations but also in determining the distribution patterns of liana species.

Application of a simplified method of chloroplast enrichment to small amounts of tissue for chloroplast genome sequencing.

PREMISE OF THE STUDY: High-throughput sequencing of genomic DNA can recover complete chloroplast genome sequences, but the sequence data are usually dominated by sequences from nuclear/mitochondrial genomes. To overcome this deficiency, a simple enrichment method for chloroplast DNA from small amounts of plant tissue was tested for eight plant species including a gymnosperm and various angiosperms. METHODS: Chloroplasts were enriched using a high-salt isolation buffer without any step gradient procedures, and enriched chloroplast DNA was sequenced by multiplexed high-throughput sequencing. RESULTS: Using this simple method, significant enrichment of chloroplast DNA-derived reads was attained, allowing deep sequencing of chloroplast genomes. As an example, the chloroplast genome of the conifer Callitris sulcata was assembled, from which polymorphic microsatellite loci were isolated successfully. DISCUSSION: This chloroplast enrichment method from small amounts of plant tissue will be particularly useful for studies that use sequencers with relatively small throughput and that cannot use large amounts of tissue (e.g., for endangered species).

Inferences of population structure and demographic history for Taxodium distichum, a coniferous tree in North America, based on amplicon sequencing analysis.

PREMISE OF THE STUDY: Studies of natural genetic variation can elucidate the genetic basis of phenotypic variation and the past population structure of species. Our study species, Taxodium distichum, is a unique conifer that inhabits the flood plains and swamps of North America. Morphological and ecological differences in two varieties, T. distichum var. distichum (bald cypress) and T. distichum var. imbricarium (pond cypress), are well known, but little is known about the level of genetic differentiation between the varieties and the demographic history of local populations. METHODS: We analyzed nucleotide polymorphisms at 47 nuclear loci from 96 individuals collected from the Mississippi River Alluvial Valley (MRAV), and Gulf Coastal populations in Texas, Louisiana, and Florida using high-throughput DNA sequencing. Standard population genetic statistics were calculated, and demographic parameters were estimated using a composite-likelihood approach. KEY RESULTS: Taxodium distichum in North America can be divided into at least three genetic groups, bald cypress in the MRAV and Texas, bald cypress in Florida, and pond cypress in Florida. The levels of genetic differentiation among the groups were low but significant. Several loci showed the signatures of positive selection, which might be responsible for local adaptation or varietal differentiation. CONCLUSIONS: Bald cypress was genetically differentiated into two geographical groups, and the boundary was located between the MRAV and Florida. This differentiation could be explained by population expansion from east to west. Despite the overlap of the two varieties' ranges, they were genetically differentiated in Florida. The estimated demographic parameters suggested that pond cypress split from bald cypress during the late Miocene.

Complex pollination of a tropical Asian rainforest canopy tree by flower-feeding thrips and thrips-feeding predators.

PREMISE OF THE STUDY: In tropical rainforests of Southeast Asia, a highly fecund thrips (Thrips spp.) responds rapidly to the mass flowering at multiple-year intervals characteristic of certain species such as the canopy tree studied here, Shorea acuminata, by feeding on flower resources. However, past DNA analyses of pollen adherent to thrips bodies revealed that the thrips promoted a very high level of self-pollination. Here, we identified the pollinator that contributes to cross-pollination and discuss ways that the pollination system has adapted to mass flowering. METHODS: By comparing the patterns of floral visitation and levels of genetic diversity in adherent pollen loads among floral visitors, we evaluated the contribution of each flower visitor to pollination. KEY RESULTS: The big-eyed bug, Geocoris sp., a major thrips predator, was an inadvertent pollinator, and importantly contributed to cross-pollination. The total outcross pollen adhering to thrips was approximately 30% that on the big-eyed bugs. Similarly, 63% of alleles examined in S. acuminata seeds and seedlings occurred in pollen adhering to big-eyed bugs; about 30% was shared with pollen from thrips. CONCLUSIONS: During mass flowering, big-eyed bugs likely travel among flowering S. acuminata trees, attracted by the abundant thrips. Floral visitation patterns of big-eyed bugs vs. other insects suggest that these bugs can maintain their population size between flowering by preying upon another thrips (Haplothrips sp.) that inhabits stipules of S. acuminata throughout the year and quickly respond to mass flowering. Thus, thrips and big-eyed bugs are essential components in the pollination of S. acuminata.

Evolutionary rate variation in two conifer species, Taxodium distichum (L.) Rich. var. distichum (baldcypress) and Cryptomeria japonica (Thunb. ex L.f.) D. Don (Sugi, Japanese cedar).

With the advance of sequencing technologies, large-scale data of expressed sequence tags and full-length cDNA sequences have been reported for several conifer species. Comparative analyses of evolutionary rates among diverse taxa provide insights into taxon-specific molecular evolutionary features and into the origin of variation in evolutionary rates within genomes and between species. Here, we estimated evolutionary rates in two conifer species, Taxodium distichum and Cryptomeria japonica, to illuminate the molecular evolutionary features of these species, using hundreds of genes and employing Chamaecyparis obtusa as an outgroup. Our results show that the mutation rates based on synonymous substitution rates (dS) of T. distichum and C. japonica are approximately 0.67 × 10(-9) and 0.59 × 10(-9)/site/year, respectively, which are 15-25 times lower than those of annual angiosperms. We found a significant positive correlation between dS and GC3. This implies that a local mutation bias, such as context dependency of the mutation bias, exists within the genomes of T. distichum and C. japonica, and/or that selection acts on synonymous sites in these species. In addition, the means of the ratios of synonymous to nonsynonymous substitution rate in the two species are almost the same, suggesting that the average intensity of functional constraint is constant between the lineages. Finally, we tested the possibility of positive selection based on the site model, and detected one candidate gene for positive selection.

Evidence of intense chromosomal shuffling during conifer evolution.

While recent advances have been gained on genome evolution in angiosperm lineages, virtually nothing is known about karyotype evolution in the other group of seed plants, the gymnosperms. Here we used high density gene-based linkage mapping to compare the karyotype structure of two families of conifers (the most abundant group of gymnosperms) separated around 290 million years ago: Pinaceae and Cupressaceae. We propose for the first time a model based on the fusion of 20 ancestral chromosomal blocks that may have shaped the modern karyotpes of Pinaceae (with n=12) and Cupressaceae (with n=11). The considerable difference in modern genome organization between these two lineages contrasts strongly with the remarkable level of synteny already reported within the Pinaceae. It also suggests a convergent evolutionary mechanism of chromosomal block shuffling that has shaped the genomes of the spermatophytes.

Development of nuclear and chloroplast microsatellite markers for the endangered conifer Callitris sulcata (Cupressaceae).

PREMISE OF THE STUDY: Microsatellite markers were developed for Callitris sulcata (Cupressaceae), an endangered conifer species in New Caledonia. METHODS AND RESULTS: Using sequencing by synthesis (SBS) of an RNA-Seq library, 15 polymorphic nuclear and chloroplast microsatellite markers were developed. When evaluated with 48 individuals, these markers showed genetic variations ranging from two to 15 alleles and expected heterozygosity ranging from 0 to 0.881. CONCLUSIONS: These markers will be useful for examining the genetic diversity and structure of remaining wild populations and improving the genetic status of ex situ populations.

Mixed Mating System Are Regulated by Fecundity in Shorea curtisii (Dipterocarpaceae) as Revealed by Comparison under Different Pollen Limited Conditions.

The maintenance of mixed mating was studied in Shorea curtisii, a dominant and widely distributed dipterocarp species in Southeast Asia. Paternity and hierarchical Bayesian analyses were used to estimate the parameters of pollen dispersal kernel, male fecundity and self-pollen affinity. We hypothesized that partial self incompatibility and/or inbreeding depression reduce the number of selfed seeds if the mother trees receive sufficient pollen, whereas reproductive assurance increases the numbers of selfed seeds under low amounts of pollen. Comparison of estimated parameters of self-pollen affinity between high density undisturbed and low density selectively logged forests indicated that self-pollen was selectively excluded from mating in the former, probably due to partial self incompatibility or inbreeding depression until seed maturation. By estimating the self-pollen affinity of each mother tree in both forests, mother trees with higher amount of self-pollen indicated significance of self-pollen affinity with negative estimated value. The exclusion of self-fertilization and/or inbreeding depression during seed maturation occurred in the mother trees with large female fecundity, whereas reproductive assurance increased self-fertilization in the mother trees with lower female fecundity.

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

BACKGROUND AND AIMS: Distribution shifts and natural selection during past climatic changes are important factors in determining the genetic structure of forest species. In particular, climatic fluctuations during the Quaternary appear to have caused changes in the distribution ranges of plants, and thus strongly affected their genetic structure. This study was undertaken to identify the responses of the conifer Cryptomeria japonica, endemic to the Japanese Archipelago, to past climatic changes using a combination of phylogeography and species distribution modelling (SDM) methods. Specifically, this study focused on the locations of refugia during the last glacial maximum (LGM). METHODS: Genetic diversity and structure were examined using 20 microsatellite markers in 37 populations of C. japonica. The locations of glacial refugia were assessed using STRUCTURE analysis, and potential habitats under current and past climate conditions were predicted using SDM. The process of genetic divergence was also examined using the approximate Bayesian computation procedure (ABC) in DIY ABC to test the divergence time between the gene pools detected by the STRUCTURE analysis. KEY RESULTS: STRUCTURE analysis identified four gene pools: northern Tohoku district; from Chubu to Chugoku district; from Tohoku to Shikoku district on the Pacific Ocean side of the Archipelago; and Yakushima Island. DIY ABC analysis indicated that the four gene pools diverged at the same time before the LGM. SDM also indicated potential northern cryptic refugia. CONCLUSIONS: The combined evidence from microsatellites and SDM clearly indicates that climatic changes have shaped the genetic structure of C. japonica. The gene pool detected in northern Tohoku district is likely to have been established by cryptic northern refugia on the coast of the Japan Sea to the west of the Archipelago. The gene pool in Yakushima Island can probably be explained simply by long-term isolation from the other gene pools since the LGM. These results are supported by those of SDM and the predicted divergence time determined using ABC analysis.