Yoshihiko Tsumura

Phylogenetic relationships of Eurasian pines (Pinus, Pinaceae) based on chloroplast rbcL, MATK, RPL20-RPS18 spacer, and TRNV intron sequences.

The sequence divergence of chloroplast rbcL, matK, trnV intron, and rpl20-rps18 spacer regions was analyzed among 32 Pinus species and representatives of six other genera in Pinaceae. The total aligned sequence length is 3570 bp. Of the four sequences examined, matK evolved much faster than rbcL in Pinus and in other Pinaceae genera. The two noncoding regions did not show more divergence than the two coding regions, especially within each Pinus subgenus. Phylogenetic analyses based on these four sequences gave consistent results and strongly supported the monophyly hypothesis for the genus Pinus and its two recognized subgenera. Pinus krempfii, the two-flat-needle pine endemic to Vietnam, was placed in subgen. Strobus and showed closer affinity to subsect. Gerardianae. The ancient character of sect. Parrya is further confirmed. However, monophyly of the sect. Parrya is not supported by our data. Among the Eurasian pines of subgen. Pinus, Mediterranean pines formed one clade and the Asian members of subsect. Sylvestres formed another. The Himalayan P. roxburghii showed considerable divergence from all the other hard pines from both regions. Pinus merkusii was distinctly separated from all the Asian members of subsect. Sylvestres. The implications of our results for Pinus classification are discussed.

Estimation of gene flow in the tropical-rainforest tree Neobalanocarpus heimii (Dipterocarpaceae), inferred from paternity analysis.

Pollen flow and population genetic structure among 30 potentially flowering individuals of Neobalanocarpus heimii, a tropical emergent tree, were investigated in a lowland tropical rainforest of Malaysia using microsatellite polymorphism. The 248 offspring in the vicinity of five reproductive trees of the 30 potentially flowering trees were used in paternity analysis for pollen‐flow study. Four primer pairs, developed in different species of dipterocarps, were adopted to detect microsatellite polymorphism. Based upon microsatellite polymorphism, pollen flow and seed migration were detected. Pollen‐flow events of more than 400 m were observed directly, based on paternity analysis in the study plot. The estimated average mating distance of the five reproductive trees was 524 m. This result suggests that reproduction of this species is mediated by a long‐distance pollinator. The haplotypes of some offspring were not compatible with the nearest reproductive tree. Thus, the results suggest that some seeds are dispersed by a seed dispersal vector. Investigation of genetic structure showed significant and negative correlation of genetic relatedness and spatial distances between the 30 potentially flowering trees, but this correlation was weak. We suggest that long‐distance gene flow and seed migration are responsible for the poorly developed genetic structure of this species.

Genetic Diversity and Outcrossing Rate between Undisturbed and Selectively Logged Forests of Shorea curtisii (Dipterocarpaceae) Using Microsatellite DNA Analysis

To evaluate the effect of flowering tree density on genetic variation in a dominant tree species, genetic diversities and outcrossing rates of Shorea curtisii (Dipterocarpaceae) were compared between undisturbed and selectively logged forests using microsatellite markers. The mean heterozygosity and the effective number of alleles per locus in the selectively logged plot were not significantly different from those of the undisturbed plot. However, the average number of alleles derived from the pollen donor was significantly lower in the selectively logged forests, indicating the possibility of a loss of genetic variation in future generations. The average outcrossing rate for 11 trees in an undisturbed plot was 96.3%, which was similar to other reports from dipterocarp species based on isozyme analysis. However, the outcrossing rate in the selectively logged plot was 52.2% for five trees, which was significantly lower than that in the undisturbed plot. This decline of outcrossing rate in the selectively logged plot may result from the much lower density of conspecific flowering trees, and the high rate of selfing may result in inbreeding depression and/or a decrease in genetic variation of this forest in future generations.

Expression analysis of ESTs derived from the inner bark of Cryptomeria japonica

To assist genetic research into Cryptomeria japonica, which is one of the most important forest tree species in Japan, expressed sequence tag (EST) analysis was carried out. The cDNA clones were isolated from a library derived from inner bark tissues. Partial sequences were obtained from 2231 clones, representing 1398 unique transcripts. Putative functions were assigned to 1583 clones, which represented 882 unique transcripts, by a Blast algorithm. Homology analysis suggested that ESTs related to cell wall formation represented about 3% of the clones. Transcripts of plant stress response genes were also abundant in the inner bark library, especially genes involved in wounding and drought responses. This indicates that the stress response systems of this tree species are similar to those of other plants, and that these systems are highly conserved among plant species. The remaining 648 clones, which represented 516 unique transcripts, did not show any significant homology to known sequences in the databases searched: these are expected to represent genes specific to Cryptomeria and, possibly, to related species.

Genome Scan to Detect Genetic Structure and Adaptive Genes of Natural Populations of Cryptomeria japonica

We investigated 29 natural populations of Cryptomeria japonica using 148 cleaved amplified polymorphic sequence markers to elucidate their genetic structure and identify candidate adaptive genes of this species. In accordance with the inferred evolutionary history of the species during and after the last glacial episode, the genetic diversity was higher in western populations than in northern populations. The results of phylogenetic and genetic structure analyses suggest that populations of the two main varieties of the species have clearly diverged from each other and that two of the examined loci are strongly associated with the differentiation between the two varieties. Using a coalescent simulation based on FST and He values, we detected five genes that had higher, and two that had lower, values than the respective 99% confidence intervals (C.I.s) that are theoretically expected intervals under a neutral infinite-island model. We also detected 13 outlier loci using a coalescent simulation based on the assumption that the 2 varieties originated from the splitting of an ancestral population. Four of these loci were detected by both methods, two of which were detected in a genetic structure analysis as loci associated with differentiation between the two varieties of the species, and are strong candidates for genes that have been subject to selection.

Selfing and inbreeding depression in seeds and seedlings of Neobalanocarpus heimii (Dipterocarpaceae).

We evaluated the degree of selfing and inbreeding depression at the seed and seedling stages of a threatened tropical canopy tree, Neobalanocarpus heimii, using microsatellite markers. Selection resulted in an overall decrease in the level of surviving selfed progeny from seeds to established seedlings, indicating inbreeding depression during seedling establishment. Mean seed mass of selfed progeny was lower than that of outcrossed progeny. Since the smaller seeds suffered a fitness disadvantage at germination in N. heimii, the reduced seed mass of selfed progeny would be one of the determinants of the observed inbreeding depression during seedling establishment. High selfing rates in some mother trees could be attributed to low local densities of reproductive individuals, thus maintenance of a sufficiently high density of mature N. heimii should facilitate regeneration and conservation of the species.

Development and polymorphism of simple sequence repeat DNA markers for Shorea curtisii and other Dipterocarpaceae species

Nine simple sequence repeat (SSR) markers were developed from Shorea curtisii using two different methods. One SSR locus was isolated by the commonly used method of screening by colony hybridization, and the other eight loci were isolated by a vectorette PCR method. Primer pairs were designed based on the sequences of all these SSR loci. Analysis of 40 individuals of S. curtisii from natural forest in Malaysia revealed that all SSR loci were polymorphic. Four SSR markers, Shc01, Shc04, Shc07 and Shc09, were highly polymorphic. We have also tested the applicability of these SSR printers to other species of Dipterocarpaceae using PCR amplification. Because the flanking region sequences of the S. curtisii SSRs were well conserved within this family, the SSR primers for S. curtisii can be applied to almost all species of Dipterocarpaceae.

Development and characterization of microsatellite markers for Cryptomeria japonica D.Don.

Abstract.Thirty four microsatellite markers for Cryptomeria japonica D. Don were developed by searching three types of library: a database of C. japonica cDNA sequences, a standard non-enriched genomic DNA library and a microsatellite-enriched library using magnetic particles. The enrichment of microsatellite sequences using magnetic particles is very efficient compared to the other two methods both in terms of the numbers of markers generated, and in the polymorphism they detect. The microsatellites developed from the genomic DNA library generally have longer repeat sequences and are more polymorphic than those from cDNA. All the developed microsatellite markers in this study showed polymorphism among 28 plus trees selected from locations scattered throughout Japan. The mean number of alleles per locus (MNA) detected in the 28 plus trees ranged from 2 to 21 with an average of 7.5. The Polymorphism Information Content (PIC) ranged from 0.160 to 0.936 with an average of 0.666. Co-dominant segregation of alleles in a three-generation pedigree of C. japonica was demonstrated at 34 microsatellite loci, and the segregation was not distorted from Mendelian expectation for all loci. In 12 out of 34 loci, a null allele was detected. Key relationships between polymorphic parameters, such as MNA and PIC, and the characteristics of microsatellite sequences, such as the longest repeat number, total repeat number and total number of nucleotides, were investigated using rank correlation coefficients, Kendalls τ. A positive correlation was found between repeat lengths and polymorphisms. The markers provide sufficient resolution for investigating gene flow within forests and seed orchards, and for genome mapping.

Development and polymorphism of microsatellite markers for Fagus crenata and the closely related species, F. japonica

Abstract We have developed microsatellite markers (SSRs) applicable to Fagus crenata using the RAHM method and investigated their polymorphisms. We also applied the SSRs in an analysis of a closely related species, F. japonica. Here we describe the isolation and characterization of nine polymorphic microsatellite markers, of which eight are applicable to both species. Among 30 individuals of each of F. crenata and F. japonica we detected a total of 79 and 77 alleles, respectively, with an average of 9.9 and 8.6 alleles per locus. The mean expected heterozygosity (He) was 0.615 (range: 0.216–0.925) in F. crenata and 0.660 in F. japonica (range: 0.259–0.827). The He values were considerably higher than those previously found for isozymes. Paternity exclusion probabilities for multiple loci, calculated over all loci, were extremely high (0.999 and 0.998 in F. crenata and F. japonica, respectively): sufficiently high to study pollen flow in both species.

Microsatellite markers reveal high allelic variation in natural populations of Cryptomeria japonica near refugial areas of the last glacial period

Using 11 microsatellite markers, we investigated the allelic variation and genetic structure of Cryptomeria japonica, across most of its natural distribution. The markers displayed high levels of polymorphism (average gene diversity=0.77, average number of alleles=24.0), in sharp contrast to the lower levels of polymorphism found in allozyme and cleaved amplified polymorphic sequence markers in previous studies. Little genetic differentiation was found among populations (FST=0.028, P<0.001), probably because the species is wind-pollinated and long-lived. No clear relationship between Nei’s genetic distances and geographical locations of the populations were found using the principal coordinate and unweighted pair-group method with arithmetic averaging analyses. The lack of such trends might be due partly to microsatellite homoplasy arising from mutation blurring the genealogical record. However, there was a trend towards high allelic diversity in five populations (Ashitaka, Ashiu, Oki-Island, Yakushima-Island-1 and -2), which are very close to, or in, refugial areas of the last glacial period as defined by Tsukada based on pollen analysis data and current climatic divisions. We postulate that these refugial populations might have been less affected by genetic drift than the other populations due to their relatively large size.