A DNA barcode reference library for the native woody seed plants of Japan.

DNA barcode databases are increasingly available for a range of organisms, facilitating the wide application of DNA barcode-based studies. Here we announce the development of a comprehensive DNA barcode reference library of Japanese native woody seed plants representing 43 orders, 99 families, 303 genera and 834 species, and comprising 77.3% of the genera and 72.2% of the species of native woody seed plants in Japan. A total of 6216 plant specimens were collected from 223 sites across the subtropical, temperate, boreal and alpine biomes in Japan with most species represented by multiple accessions. This reference library utilized three chloroplast DNA regions (rbcL, trnH-psbA and matK) and consists of 14,403 barcode sequences. Individual regions varied in their identification rates, with species-level and genus-level rates for rbcL, trnH-psbA and matK based on blast being 57.4%/96.2%, 78.5%/99.1% and 67.8%/98.1%, respectively. Identification rates were higher using region combinations, with total species-level rates for two region combinations (rbcL & trnH-psbA, rbcL & matK and trnH-psbA & matK) ranging between 90.6% and 95.8%, and for all three regions being equal to 98.6%. Genus-level identification rates were even higher, ranging between 99.7% and 100% for two region combinations and being 100% for the three regions. These results indicate that this DNA barcode reference library is an effective resource for investigations of native woody seed plants in Japan using DNA barcodes and provides a useful template for the development of libraries for other components of the Japanese flora.

Genetic diversity, structure, and demography of Pandanus boninensis (Pandanaceae) with sea drifted seeds, endemic to the Ogasawara Islands of Japan: Comparison between young and old islands.

Pandanus boninensis, endemic to the Ogasawara Islands, Japan, is distributed on both the older Bonin and younger Volcano Islands. In this study, we conducted population genetic analyses of P. boninensis on these islands to examine the population diversity and structure across old and young islands, to assess potential differences in population demography with island age, and to collect any evidence of migration between old and young islands. We found that the genetic diversity of expressed sequence tag (EST)-based microsatellite (SSR) markers, the nucleotide diversity of nuclear DNA sequences, and the haplotype diversity of chloroplast DNA on young islands were lower than those on old islands. Clustering analyses of EST-SSR indicated that populations on old islands were strongly diverged from those on young islands. Approximate Bayesian computation analysis of EST-SSR suggested that population expansion occurred on old islands while population reduction occurred on young islands. We also found evidence of migration among old islands (mostly from south to north), while it appears that there have been very few migration events between old and young islands. These differences could be due to the fact that young islands tend to be geographically isolated and support smaller populations that began a shorter time ago from limited founders. The P. boninensis populations on the Volcano Islands are interesting from an evolutionary perspective as they constitute a classic example of the early stages of progressive colonization on oceanic islands with small effective population sizes and low genetic diversity.

Genetic structure of Sakhalin spruce (Picea glehnii) in northern Japan and adjacent regions revealed by nuclear microsatellites and mitochondrial gene sequences.

The genetic structure of Sakhalin spruce (Picea glehnii) was studied across the natural range of the species, including two small isolated populations in south Sakhalin and Hayachine, by using six microsatellite loci and maternally inherited mitochondrial gene sequences. We also analyzed P. jezoensis, a sympatric spruce in the range. Genetic diversity of P. glehnii was higher in central Hokkaido and the lowest in the Hayachine. Bayesian clustering and principal coordinate analysis by using the microsatellites indicated that the Hayachine was clearly distinct from other populations, implying that it had undergone strong genetic drift since the last glacial period. P. glehnii harbored four mitochondrial haplotypes, two of which were shared with P. jezoensis. One of the two was observed without geographical concentration, suggesting its derivation from ancestral polymorphism. Another was observed in south Sakhalin and in P. jezoensis across Sakhalin. The Bayesian clustering--by using four microsatellite loci, including P. jezoensis populations--indicated unambiguous species delimitation, but with possible admixture of P. jezoensis genes into P. glehnii in south Sakhalin, where P. glehnii is abundantly overwhelmed by P. jezoensis; this might explain the occurrence of introgression of the haplotype of P. jezoensis into P. glehnii.

Genetic differentiation in Elaeocarpus photiniifolia (Elaeocarpaceae) associated with geographic distribution and habitat variation in the Bonin (Ogasawara) Islands.

Gene flow between populations in different environmental conditions can be limited due to divergent natural selection, thus promoting genetic differentiation. Elaeocarpus photiniifolia, an endemic tree species in the Bonin Islands, is distributed in two types of habitats, dry scrubs and mesic forests. We aim to elucidate the genetic differentiation in E. photiniifolia within and between islands and between the habitat types. We investigated genotypes of 639 individuals from 19 populations of E. photiniifolia and its closely-related E. sylvestris at 24 microsatellite loci derived from expressed sequence tags. The data revealed genetic differentiation (1) between E. photiniifolia and E. sylvestris (0.307 ≤ F ST ≤ 0.470), (2) between the E. photiniifolia populations of the Chichijima and Hahajima Island Groups in the Bonin Islands (0.033 ≤ F ST ≤ 0.121) and (3) between E. photiniifolia populations associated with dry scrubs and mesic forests in the Chichijima Island Group (0.005 ≤ F ST ≤ 0.071). Principal coordinate analysis and Bayesian clustering analysis also showed that genetically distinct groups were associated with the habitat types, and isolation by distance was not responsible for the genetic differentiation. These findings suggest that E. photiniifolia is divided into genetically differentiated groups associated with different environmental conditions in the Bonin Islands.

Clone identification in Japanese flowering cherry (Prunus subgenus Cerasus) cultivars using nuclear SSR markers.

Numerous cultivars of Japanese flowering cherry (Prunus subgenus Cerasus) are recognized, but in many cases they are difficult to distinguish morphologically. Therefore, we evaluated the clonal status of 215 designated cultivars using 17 SSR markers. More than half the cultivars were morphologically distinct and had unique genotypes. However, 22 cultivars were found to consist of multiple clones, which probably originate from the chance seedlings, suggesting that their unique characteristics have not been maintained through propagation by grafting alone. We also identified 23 groups consisting of two or more cultivars with identical genotypes. Most members of these groups were putatively synonymously related and morphologically identical. However, some of them were probably derived from bud sport mutants and had distinct morphologies. SSR marker analysis provided useful insights into the clonal status of the examined Japanese flowering cherry cultivars and proved to be a useful tool for cultivar characterization.

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.

Phylogeography of the Korean pine (Pinus koraiensis) in northeast Asia: inferences from organelle gene sequences.

Range-wide genetic variation of Korean pine (Pinus koraiensis) was assessed using maternally inherited mtDNA and paternally inherited cpDNA for 16 natural populations throughout northeast Asia in order to study its phylogeographical history during the Quaternary. The cpDNA variation indicated that there was no difference between populations on the Asian continent and those in the Japanese archipelago. In contrast, the mtDNA variation indicated that there was significant difference between the populations from the two regions, with each region having a different lineage. The continental populations exhibited no diversity in the mtDNA examined despite the species' current extensive range and large populations. Conversely, while the Korean pine is rare in Japan, the Japanese populations exhibited greater levels of mtDNA diversity (H (T) = 0.502). The higher mtDNA diversity and evidence from numerous Korean pine macrofossil remains dated to the Pleistocene and recovered various sites in Japan suggest that the Japanese archipelago once served as a refugium to a much larger Korean pine population with a more extensive range than is the case today. The presence of the single mtDNA haplotype across the Asian continent suggests that the present widespread populations could have expanded from a single refugium population after the last glacial periods.

Development of EST-SSR markers for Elaeocarpus photiniifolia (Elaeocarpaceae), an endemic taxon of the Bonin Islands.

PREMISE OF THE STUDY: Expressed sequence tag (EST)-derived microsatellite markers were developed for Elaeocarpus photiniifolia, an endemic taxon of the Bonin Islands. METHODS AND RESULTS: Initially, a complementary DNA (cDNA) library was constructed by de novo pyrosequencing of total RNA extracted from a seedling. A total of 267 primer pairs were designed from the library. Of the 48 tested loci, 25 loci were polymorphic among 41 individuals representing the entire geographical range of the species, with the number of alleles per locus and expected heterozygosity ranging from two to 14 and 0.09 to 0.86, respectively. Most loci were transferable to a related species, E. sylvestris. CONCLUSIONS: The developed markers will be useful for evaluating the genetic structure of E. photiniifolia.

Isolation and characterization of EST-SSR markers in Schima mertensiana (Theaceae) using pyrosequencing technology.

PREMISE OF THE STUDY: Microsatellite or simple sequence repeat (SSR) markers from expressed sequence tags (ESTs) in Schima mertensiana, an evergreen tree that is endemic to the Bonin Islands, were developed to facilitate evaluation of the species's population genetic structure. METHODS AND RESULTS: Using 149831 ESTs derived by pyrosequencing S. mertensiana complementary DNA (cDNA), 683 EST-SSRs were developed, 48 of which were tested for amplification. Thirty of the 48 showed clear amplification and detected polymorphism in sampled populations, with mean expected heterozygosity at the amplified loci ranging from 0.045 to 0.874 (average 0.545). CONCLUSIONS: The markers developed in this study can be used for future studies of the genetic structure of S. mertensiana and related species.

Rapid development of microsatellite markers for Pandanus boninensis (Pandanaceae) by pyrosequencing technology.

PREMISE OF THE STUDY: To facilitate rapid development of microsatellite or simple sequence repeat (SSR) markers, an expressed sequence tags (EST) database was constructed for Pandanus boninensis, an evergreen tree endemic to the Bonin Islands, using pyrosequencing technology. METHODS AND RESULTS: We designed primers for 340 EST-SSRs identified from 109620 pyrosequencing reads, 48 of which were tested for PCR amplification. Thirty-four primers provided clear amplification, and 26 of those 34 displayed clear polymorphic patterns in sampled populations, with mean expected heterozygosity at the amplified loci ranging from 0.022 to 0.742 (average 0.262). CONCLUSIONS: The developed markers are promising tools for future genetic studies of P. boninensis and related species.

Microsatellite markers derived from Calophyllum inophyllum (Clusiaceae) expressed sequence tags.

PREMISE OF THE STUDY: Robust markers are required (inter alia) for assessing origins of Calophyllum inophyllum populations on the Bonin Islands, Japan. Therefore, informative expressed sequence tag (EST)-based microsatellite or simple sequence repeat (SSRs) markers in the species were sought. METHODS AND RESULTS: Using 135378 ESTs derived from de novo pyrosequencing, primers for 475 EST-SSRs were developed, 48 of which were tested for PCR amplification. Thirty-six of the 48 primers showed clear amplification, with 23 displaying polymorphism in sampled populations. Expected heterozygosity in the samples from the Bonin Islands and Ryukyu Islands populations ranged from 0.041 to 0.697 and from 0.041 to 0.773, respectively. CONCLUSIONS: As EST-SSRs are potentially tightly linked with functional genes, and reportedly more transferable to related species than anonymous genomic SSRs, the developed primers have utility for future studies of the origins, genetic structure, and conservation of C. inophyllum and related species.

Isolation and characterization of microsatellite loci in an endangered tree species, Afzelia xylocarpa (Kurz) Craib (Caesalpinioideae).

Afzelia xylocarpa is listed as an endangered species on the International Union for the Conservation of Nature and Natural Resources World list of Threatened Trees, due to overexploitation for its valuable timber and habitat loss. We isolated eight polymorphic microsatellite loci from A. xylocarpa using dual suppression polymerase chain reaction technique. These loci provide microsatellite markers with high polymorphism, as the number of alleles ranged from two to seven, and the estimate of gene diversity was between 0.285 and 0.795. The markers are now available for more detailed investigation of population genetic structure and gene flow among A. xylocarpa populations.

Adaptive and slightly deleterious evolution in a conifer, Cryptomeria japonica.

In order to evaluate effects of the population structure and natural selection on organisms having long generation times, we surveyed DNA polymorphisms at five loci encoding 9-cis-epoxycarotenoid dioxygenase (NCED), ammonium transporter, calmodulin, aquaporin, and the second major allergen with polymethylgalacturonase enzyme activity in the pollen (Cryj2) in a conifer, Cryptomeria japonica. The average nucleotide diversity at silent sites across 12 loci including the previously analyzed seven loci was 0.0044. The population recombination rate (4Nr, where N and r are the effective population size and recombination rate per base per generation, respectively) was estimated as 0.00046 and a slow reduction in the population size was indicated, according to the maximum likelihood method implemented in LAMARC. At NCED, the McDonald-Kreitman (MK) test revealed an excess of replacement polymorphisms, suggesting contributions of slightly deleterious mutations. In contrast, the MK test revealed an excess of replacement divergence at Cryj2 and a maximum likelihood approach using the PAML package revealed that certain amino acid sites had a nonsynonymous/synonymous substitution rate ratio (omega) > 4.0, indicating adaptive evolution at this locus. The overall analysis of the 12 loci suggested that adaptive, neutral, and slightly deleterious evolution played important roles in the evolution of C. japonica.

Development and characterization of microsatellite markers for Hibiscus glaber Matsum. ex Nakai, an endemic tree species of the oceanic Bonin Islands, Japan.

Polymorphic microsatellite markers were developed for Hibiscus glaber, an endemic tree of the Bonin Islands. Eighty-seven of the 208 sequences from an enriched library were unique and containing microsatellites. Ten loci were proved to be highly polymorphic among 78 individuals from the Nishi-jima Island. Total exclusionary powers for the first and the second parents were 99.989% and 99.999%, respectively. Nine loci also amplified single fragment from genomic DNA of H. tiliaceus, a related and widespread congener. Our markers can be reliably used for the estimation of current gene flow within/among populations of the two woody Hibiscus species.

Determination of the genetic structure of remnant Morus boninensis Koidz. trees to establish a conservation program on the Bonin Islands, Japan.

BACKGROUND: Morus boninensis, is an endemic plant of the Bonin (Ogasawara) Islands of Japan and is categorized as "critically endangered" in the Japanese red data book. However, little information is available about its ecological, evolutionary and genetic status, despite the urgent need for guidelines for the conservation of the species. Therefore, we adopted Moritz's MU concept, based on the species' current genetic structure, to define management units and to select mother tree candidates for seed orchards. RESULTS: Nearly all individuals of the species were genotyped on the basis of seven microsatellite markers. Genetic diversity levels in putative natural populations were higher than in putative man-made populations with the exception of those on Otouto-jima Island. This is because a limited number of maternal trees are likely to have been used for seed collection to establish the man-made populations. A model-based clustering analysis clearly distinguished individuals into nine clusters, with a large difference in genetic composition between the population on Otouto-jima Island, the putative natural populations and the putative man-made populations. The Otouto-jima population appeared to be genetically differentiated from the others; a finding that was also supported by pairwise FST and RST analysis. Although multiple clusters were detected in the putative man-made populations, the pattern of genetic diversity was monotonous in comparison to the natural populations. CONCLUSION: The genotyping by microsatellite markers revealed strong genetic structures. Typically, artificial propagation of this species has ignored the genetic structure, relying only on seeds from Otouto-jima for replanting on other islands, because of a problem with inter-specific hybridization on Chichi-jima and Haha-jima Islands. However, this study demonstrates that we should be taking into consideration the genetic structure of the species when designing a propagation program for the conservation of this species.

Contrasting patterns of DNA variation in natural populations of two related conifers, Cryptomeria japonica and Taxodium distichum (Cupressaceae sensu lato).

We investigated DNA variation within and between two closely related conifers, Cryptomeria japonica and Taxodium distichum, at nuclear loci encoding ferredoxin, glutamyl-tRNA reductase, lycopene beta cyclase, and phosphoribosylanthranilate transferase. Average nucleotide diversity at silent sites was estimated to be 0.0035 (SE 0.0012) in C. japonica and 0.0058 (SE 0.0006) in T. distichum. One population in C. japonica was differentiated from the others but generally there was not much differentiation among populations or varieties within the two species. However, the two species seemed to differ in frequency spectra of DNA polymorphisms. Excesses of intermediate-frequency variants were found in C. japonica, whereas excesses of both rare and high-frequency variants were found in T. distichum, which suggested different histories of population structures in the two species. Deviations from the standard neutral expectations in DNA polymorphisms were found by applications of neutrality tests. The results show that actions of selection to respective loci seem to differ between the two species, indicating differences of interaction among evolutionary factors.

A consensus linkage map for sugi (Cryptomeria japonica) from two pedigrees, based on microsatellites and expressed sequence tags.

A consensus map for sugi (Cryptomeria japonica) was constructed by integrating linkage data from two unrelated third-generation pedigrees, one derived from a full-sib cross and the other by self-pollination of F1 individuals. The progeny segregation data of the first pedigree were derived from cleaved amplified polymorphic sequences, microsatellites, restriction fragment length polymorphisms, and single nucleotide polymorphisms. The data of the second pedigree were derived from cleaved amplified polymorphic sequences, isozyme markers, morphological traits, random amplified polymorphic DNA markers, and restriction fragment length polymorphisms. Linkage analyses were done for the first pedigree with JoinMap 3.0, using its parameter set for progeny derived by cross-pollination, and for the second pedigree with the parameter set for progeny derived from selfing of F1 individuals. The 11 chromosomes of C. japonica are represented in the consensus map. A total of 438 markers were assigned to 11 large linkage groups, 1 small linkage group, and 1 nonintegrated linkage group from the second pedigree; their total length was 1372.2 cM. On average, the consensus map showed 1 marker every 3.0 cM. PCR-based codominant DNA markers such as cleaved amplified polymorphic sequences and microsatellite markers were distributed in all linkage groups and occupied about half of mapped loci. These markers are very useful for integration of different linkage maps, QTL mapping, and comparative mapping for evolutional study, especially for species with a large genome size such as conifers.

DNA variation in a conifer, Cryptomeria japonica (Cupressaceae sensu lato).

We investigated the nucleotide variation of a conifer, Cryptomeria japonica, and the divergence between this species and its closest relative, Taxodium distichum, at seven nuclear loci (Acl5, Chi1, Ferr, GapC, HemA, Lcyb, and Pat). Samples of C. japonica were collected from three areas, Kantou-Toukai, Hokuriku, and Iwate. No apparent geographic differentiation was found among these samples. However, the frequency spectrum of the nucleotide polymorphism revealed excesses of intermediate-frequency variants, which suggests that the population was not panmictic and a constant size in the past. The average nucleotide diversity, pi, for silent sites was 0.00383. However, values of pi for silent sites vary among loci. Comparisons of polymorphism to divergence among loci (the HKA test) showed that the polymorphism at the Acl5 locus was significantly lower. We also observed a nearly significant excess of replacement polymorphisms at the Lcyb locus. These results suggested possibilities of natural selection acting at some of the loci. Intragenic recombination was detected only once at the Chi1 locus and was not detected at the other loci. The low level of population recombination rate, 4Nr, seemed to be due to both low level of recombination, r, and small population size, N.

Molecular evolution of nuclear genes in Cupressacea, a group of conifer trees.

We surveyed the molecular evolutionary characteristics of 11 nuclear genes from 10 conifer trees belonging to the Taxodioideae, the Cupressoideae, and the Sequoioideae. Comparisons of substitution rates among the lineages indicated that the synonymous substitution rates of the Cupressoideae lineage were higher than those of the Taxodioideae. This result parallels the pattern previously found in plastid genes. Likelihood-ratio tests showed that the nonsynonymous-synonymous rate ratio did not change significantly among lineages. In addition, after adjustments for lineage effects, the dispersion indices of synonymous and nonsynonymous substitutions were considerably reduced, and the latter was close to 1. These results indicated that the acceleration of evolutionary rates in the Cupressoideae lineage occurred in both the nuclear and plastid genomes, and that generally, this lineage effect affected synonymous and nonsynonymous substitutions similarly. We also investigated the relationship of synonymous substitution rates with the nonsynonymous substitution rate, base composition, and codon bias in each lineage. Synonymous substitution rates were positively correlated with nonsynonymous substitution rates and GC content at third codon positions, but synonymous substitution rates were not correlated with codon bias. Finally, we tested the possibility of positive selection at the protein level, using maximum likelihood models, assuming heterogeneous nonsynonymous-synonymous rate ratios among codon (amino acid) sites. Although we did not detect strong evidence of positively selected codon sites, the analysis suggested that significant variation in nonsynonymous-synonymous rate ratio exists among the sites. The most likely sites for action of positive selection were found in the ferredoxin gene, which is an important component of the apparatus for photosynthesis.