Differences in female reproductive success between female and hermaphrodite individuals in the subdioecious shrub Eurya japonica (Theaceae).

Subdioecy is thought to occupy a transitional position in the gynodioecy-dioecy pathway, explaining one of the evolutionary routes from hermaphroditism to dioecy. Quantifying any female reproductive advantage of females versus hermaphrodites is fundamental to examining the spectrum between subdioecy and dioecy; however, this is challenging, as multiple interacting factors, such as pollen limitation and resource availability, affect plant reproduction. We compared the female reproductive success of females and hermaphrodites via a field experiment in which we hand-pollinated individuals of the subdioecious shrub Eurya japonica of similar size growing under similar light conditions. Effects of pollen limitation and seed quality were also evaluated through comparing the results of hand- and natural-pollination treatments and performing additional laboratory and greenhouse experiments. Overall, females had higher fruit set and produced heavier fruit and more seeds than hermaphrodites, and these results were more pronounced for hand-pollinated than for natural-pollinated plants of both sexes. We also found that seeds naturally produced by females had a higher mean germination rate. These results indicate that females had a pronounced advantage in female reproductive success under conditions of no pollen limitation. The sexual difference in the degree of pollen limitation suggests a pollinator-mediated interaction, whereas the higher female reproductive success of females even under natural conditions implies that E. japonica is a good model species for elucidating the later stages of the gynodioecy-dioecy pathway.

Estimation of outcrossing rates at hierarchical levels of fruits, individuals, populations and species in Magnolia stellata.

In plant species with mixed mating systems, differences in diverse factors-including their pollination system, flowering phenology, life form and susceptibility to inbreeding depression-cause variation in outcrossing rates among fruits within individuals, among individuals within populations and among populations within species. To quantify this hierarchical variation, we examined outcrossing rates at the seed stage in five populations of Magnolia stellata, a self-compatible, insect-pollinated and protogynous tree species. For this purpose, we sampled 1498 seeds within 204 fruits obtained from 56 individuals of the five populations, determined genotypes of the sampled seeds and maternal individuals at six polymorphic microsatellite loci, then estimated outcrossing rates and their variance components at four hierarchical levels (fruits, individuals, populations and species) using a nested analysis of variance-type linear model with a Bayesian approach. The species-level outcrossing rate was 0.730 (95% credible interval, 0.595-0.842), indicating that this species has a mixed mating system. Outcrossing rates were not significantly different among populations, but were significantly different among individuals within populations. Variance components at the levels of individual and fruit were statistically supported and were highest for the former. Thus, factors influencing outcrossing rates at the individual level, such as differences in flowering phenology and early-stage inbreeding depression, appear to have important effects within these M. stellata populations, but not among them. The method of hierarchically estimating outcrossing rates using a Bayesian approach, as applied in this study, is compared with conventional methods for estimating outcrossing rates, and the statistical properties of the Bayesian approach are discussed.

Effects of seed- and pollen-mediated gene dispersal on genetic structure among Quercus salicina saplings.

We evaluated the effects of seed- and pollen-mediated gene dispersal on genetic structure among Quercus salicina saplings. Parentage analysis using 10 microsatellite markers indicated that the 111 adult trees located within a 11.56 ha plot in the Tatera Forest Reserve, Japan, included only one parent of 44.2% and both parents of 40.7% of the 226 saplings located in a 1-ha core plot at its center. Coancestry (F(ij)) estimates indicated that there was strong genetic structure among the saplings. The numbers of pairs of full- and half-siblings were high among neighboring saplings, suggesting that there was strong maternal half-sibling family structure among the saplings around their seed parents, probably generated by the spatially limited seed dispersal and the small extent of overlapping seed shadows owing to the low density of adults. The frequencies also suggest that the maternal half-sibling families are interspersed with full-siblings, produced by correlated mating, probably because pollination frequency depends on the distance between parents. The frequencies of pairs of half-siblings decreased as the distance between saplings increased, but did not fall to zero even at distances up to the 90-95 m class, suggesting that paternal half-siblings originating from correlated paternity were widely distributed owing to extensive pollen flow. We separately examined the genetic structure for maternal and paternal alleles in the saplings. Unsurprisingly, very strong genetic structure was detected for maternal alleles. However, weak (but significant) genetic structure was also detected for paternal alleles. Therefore, pollen dispersal may affect the extent of genetic structure as well as seed dispersal.

Genetic variation and differentiation in populations of a threatened tree, Magnolia stellata: factors influencing the level of within-population genetic variation.

Genetic variation and differentiation in Magnolia stellata were studied in 20 populations distributed across most of the species' range using 10 microsatellite markers, and the factors influencing their levels of within-population genetic variation were examined. Generally, populations distributed intermittently from southern Gifu to central Aichi Prefectures showed substantially higher levels of genetic variation (exceptions included populations located at unusually high altitude sites or western and southern edges of the range) than more isolated populations on the Atsumi Peninsula of southern Aichi Prefecture and in northern Mie Prefecture. Significant isolation-by-distance patterns were detected in genetic differentiation among the studied populations, and a neighbor-joining tree based on D(A) distances among the populations reflected well the geographical positions of the populations. The level of within-population genetic variation was significantly influenced not only by the size of the populations (represented by the number of reproductive individuals) but also by their degree of isolation (represented by the number of populations within a radius of 0.5 km around them). Therefore, areas within radii of 0.5 km may encompass M. stellata metapopulations, in which gene flow may usually occur. We suggest that this area may be a suitable standard for constructing conservation units for the species.

Effects of kin-structured seed dispersal on the genetic structure of the clonal dioecious shrub ilex leucoclada.

Nonrandom patterns of gene dispersal have been identified as possible causes of genetic structuring within populations. Attempts to model these patterns have generally focused solely on the effects of isolation by distance, but the processes involved are more complex than such modeling suggests. Here, we extend considerations of gene dispersal processes beyond simple isolation by distance effects by directly evaluating the effects of kin-structured gene dispersal mediated by the group dispersal of related seeds within fruits (i.e., kin-structured seed dispersal) by birds on genetic structure in Ilex leucoclada, a clonal dioecious shrub. To examine the genetic structure patterns, we established two 30x30 m plots (one with immature soils in old-growth forest and one in secondary forest, designated IM and SC, respectively) with different I. leucoclada stem densities. In these two plots 145 and 510 stems were found, representing 78 and 85 genets, respectively, identified by analyzing their genotypes at eight microsatellite loci. The clonal structure was stronger in the SC plot than in the IM plot. Correlograms of coancestry for genets in both plots exhibited significant, positive, high values in the shortest distance class, indicating the presence of strong genetic structure. However, Sp statistics revealed that the pattern of the genetic structure differed between the plots. In addition, to estimate the family structure within fruits, we sampled forty fruits, in total, from 15 randomly selected plants in the area around the IM and SC plots, and found that 80% of the fruits were multiseeded and 42-100% of the multiseeded fruits contained at least one pair of full sibs. Simulations based on these estimates demonstrated that the group dispersal of related seeds produced through correlated mating both within and across fruits, but not unstructured half-sib dispersal, could generate the observed magnitude and trends of genetic structure found in the IM plot. Furthermore, in addition to kin-structured seed dispersal, isolation by distance processes is also likely to promote genetic substructuring in the SC plot. After discussing possible ecological factors that may have contributed to the observed genetic structure, we contrast our results with those predicted by general isolation by distance models, and propose that kin-structured seed dispersal should promote some evolutionary phenomena, and thus should be incorporated, where appropriate, in models of gene dispersal in natural plant populations.

Effects of canopy gaps on the genetic structure of Camellia japonica saplings in a Japanese old-growth evergreen forest.

The genetic structure of Camellia japonica saplings was investigated in relation to canopy conditions in an old-growth evergreen forest in Tsushima, Japan. To elucidate effects of canopy gaps on genetic structure, a 1 ha study site was divided into 20 x 20 m quadrats, which were classified into a gap quadrats (GAP), closed canopy quadrats (CLS) and mixed quadrats. Five GAP quadrats and six CLS quadrats were analyzed separately. Isolation-by-distance was tested by examining the correlation between genetic distance and geographic distance. A significant positive correlation was detected for GAP quadrats, whilst that for CLS quadrats was significantly smaller and not significantly different from zero. On the other hand, an analysis using Moran's I spatial autocorrelation coefficients indicates that the genetic structure is weaker in GAP quadrats than in CLS quadrats in short distance classes. The values were significantly positive for both types of quadrat. These results, along with our field observations on flowering, suggest that canopy gaps affect the genetic structure of C. japonica saplings in two distinct ways. First, canopy gaps may promote flowering and mating in an isolation-by-distance manner within canopy gaps. Second, canopy gaps may promote seed production and resulting overlap in seed shadows may weaken fine-scale genetic structures.

Interannual genetic heterogeneity of pollen pools accepted by Quercus salicina individuals.

Since flowering often varies among years in wind-pollinated woody species, the genetic composition of pollen pools accepted by seed parents can differ between years. The interannual heterogeneity of pollen flow may be important for maintaining genetic diversity within populations because it can increase genetic variation within populations and the effective sizes of the populations. In this study we examined heterogeneity, using paternity analysis and analysis of molecular variance, in the genetic composition of pollen pools among different reproductive years for six Quercus salicina seed parents in an 11.56-ha plot in a temperate old-growth evergreen broadleaved forest. The genotypes at seven microsatellite loci were determined for 111 adult trees and 777 offspring of the six seed parents in 2-5 reproductive years. Genetic differentiation of pollen pools among different reproductive years for each seed parent was significant over all seed parents and for each of four seed parents that were analysed for more than 2 years, but not for either of the other two seed parents (analysed for 2 years). For both the pollen pools originating from inside the plot and those originating from outside it, genetic differentiation among different reproductive years for each seed parent was significant over all seed parents. However, among-year genetic differentiation in the pollen pools originating from within the plot was detected for all four of the seed parents that were analysed for more than 2 years, but for only one of the four in the pools originating from outside the plot. Genetic diversity (estimated as allelic richness and gene diversity) was higher for pollen pools over all reproductive years than for pollen pools in single years. These results indicate that the year-to-year genetic variation of pollen pools increases genetic diversity in offspring and is strongly affected by the variation in pollen parents within the plot because of their high pollination contributions. The high year-to-year variation in pollen parents within the plot and overall supports the hypothesis that the offspring produced across years represent a larger genetic neighbourhood.

Half-sib family structure of Fagus crenata saplings in an old-growth beech-dwarf bamboo forest.

Half-sib family structure of Fagus crenata saplings was examined in an old-growth beech-dwarf bamboo (Sasa spp.) forest using microsatellite genotypes in parentage analysis to identify the half-sib families in two 50 x 50 m plots: one with 36 adults, 641 saplings and no Sasa cover, the other with 21 adults, 61 saplings and Sasa cover. For large proportions of the saplings (44.6% and 75.4%, respectively) both of their parents were found within the same plot, indicating that pollination events frequently involved short-range pollen dispersal, probably because of the high density of adults in the study population. Although almost all of the adults had half-sib families, the number of offspring in the families was highly variable. In the plot with no Sasa cover, the variation in the number of offspring was significantly explained by the size of parents, i.e. the reproductive success is higher for large adults than for small adults. The half-sibs were aggregately distributed around their parents and the distribution overlapped among different half-sib families, which may be due to the limited seed dispersal and overlapping seed shadows of this species. As expected, there was weak genetic structure in the plot. By contrast, in the plot with Sasa cover, the half-sibs were distributed sparsely with a low density, and the degree of genetic structure was very weak. The difference in the half-sib family structure and genetic structure among saplings presumably reflects the difference in density that should be affected by regeneration dynamics associated with environmental conditions.

Size distribution and genetic structure in relation to clonal growth within a population of Magnolia tomentosa Thunb. (Magnoliaceae).

To establish a baseline for conservation of a threatened clonal tree, Magnolia tomentosa, we investigated size distribution and genetic structure within a population, using six microsatellite markers. Within the study site, 1044 living ramets (stems) were distinguished into 175 genets (individuals). The mean number of ramets per genet was 5.97, and 76% of all genets had multiple ramets. Genets, which apparently produced new ramets through sprouting and layering, were generally composed of several large ramets and many small ramets. Spatial autocorrelation analysis of microsatellite alleles revealed positive autocorrelation over short distances for both ramets and genets. The Moran's I-value of ramets in the shortest distance class was 3.8 times larger than that of genets, reflecting the effect of clonal growth. To analyse the size-class differences in genetic structure, the 175 genets were separated into two size classes, small and large. The correlogram for the small genets exhibited positive spatial autocorrelation in the shortest distance class, but this was not the case for the correlogram for the large genets, indicating that genetic structure is weakened or lost through self-thinning as the genets grow. The FIS value over all loci for the small genets was positive and deviated significantly from zero, while the corresponding value for the large genets was close to zero. The excess homozygotes in the small genets may be the result of genetic substructuring and/or inbreeding, and the reduction in homozygote frequency from the small to large genets may be because of loss of genetic structure and/or inbreeding depression.

Heterogeneous genetic structure in a Fagus crenata population in an old-growth beech forest revealed by microsatellite markers.

The within-population genetic structure of Fagus crenata in a 4-ha plot (200 x 200 m) of an old-growth beech forest was analysed using microsatellite markers. To assess the genetic structure, Moran's I spatial autocorrelation coefficient was calculated. Correlograms of Moran's I showed significant positive values less than 0.100 for short-distance classes, indicating weak genetic structure. The genetic structure within the population is created by limited seed dispersal, and is probably weakened by overlapping seed shadow, secondary seed dispersal, extensive pollen flow and the thinning process. Genetic structure was detected in a western subplot of 50 x 200 m with immature soils and almost no dwarf bamboos (Sasa spp.), where small and intermediate-sized individuals were distributed in aggregations with high density because of successful regeneration. By contrast, genetic structure was not found in an eastern subplot of the same size with mature soils and Sasa cover, where successful regeneration was prevented, and the density of the small and intermediate-sized individuals was low. Moreover, genetic structure of individuals in a small-size class (diameter at breast height < 12 cm) was more obvious than in a large-size class (diameter at breast height >/= 12 cm). The apparent genetic structure detected in the 4-ha plot was therefore probably the result of the structure in the western portion of the plot and in small and intermediate-sized individuals that successfully regenerated under the favourable environment. The heterogeneity in genetic structure presumably reflects variation in the density that should be affected by differences in regeneration dynamics associated with heterogeneity in environmental conditions.

Comparison of the fine-scale genetic structure of three dipterocarp species.

We investigated the fine-scale genetic structure of three tropical-rainforest trees, Hopea dryobalanoides, Shorea parvifolia and S. acuminata (Dipterocarpaceae), in Peninsular Malaysia, all of which cooccurred within a 6-ha plot in Pasoh Forest Reserve. A significant genetic structure was found in H. dryobalanoides, weaker (but still significant) genetic structure in S. parvifolia and nonsignificant structure in S. acuminata. Seeds of all three species are wind dispersed, and their flowers are thought to be insect pollinated. The most obvious difference among these species is their height: S. parvifolia and S. acuminata are canopy species, whereas H. dryobalanoides is a subcanopy species. Clear differences were also found among these species in their range of seed dispersal, which depends on the height of the release point; so taller trees disperse their seed more extensively. The estimates of seed dispersal area were consistent with the degree of genetic structure found in the three species. Therefore, tree height probably had a strong influence on the fine-scale genetic structure of the three species.

Genetic diversity of nuclear and mitochondrial genomes in Pinus parviflora Sieb. & Zucc. (Pinaceae) populations.

Genetic diversities of the nuclear and mitochondrial genomes in Pinus parviflora were studied in 16 populations, which were distributed across most of the species' range in Japan. Six mitochondrial DNA haplotypes were identified among the 16 populations. The intrapopulation diversity of allozymes was similar to that of other endemic woody species (H(S)=0.259). Although P. parviflora is distributed in discrete populations, differentiation between these was very low (G(ST)=0.044). In addition, the extent of genetic differentiation between two varieties (var. pentaphylla and var. parviflora) was extremely low (G(VT)=0.001). Intrapopulation diversity of mitochondrial DNA was also very low (H(S)=0.098), but population differentiation was high (G(ST)=0.863). Moreover, the distribution of haplotypes reflected the taxonomic differences between P. parviflora var. pentaphylla and var. parviflora. The populations of var. pentaphylla and var. parviflora contained different haplotypes. Differing modes of inheritance may account for the differences in nuclear and mitochondrial genetic diversity.

Clonal diversity and genetic differentiation in Ilex leucoclada M. patches in an old-growth beech forest.

We investigated clonal diversity within patches of Ilex leucoclada and genetic variation within and among patches using random amplified polymorphic DNA (RAPD) markers in a 1-ha plot within an old-growth beech forest. We found 38 patches that exhibited a clumped distribution in the middle of the plot. We identified a total of 166 RAPD phenotypes among the 215 stems sampled from 27 patches that were completely within the plot. The population showed high clonal diversity within patches (mean number of genets relative to number of stems = 0.79; mean Simpson's D = 0.89). Variation in RAPD phenotypes among patches was highly significant (PhiST in the molecular variance analysis = 0.316, P < 0.001), indicating genetic differentiation among patches. Pairwise genetic distances, PhiST, among patches did not correlate with geographical distances among patches. The cluster analysis based on the genetic distances showed few clear clusters of patches, indicating no spatial genetic structure among patches. High levels of clonal diversity both within patches and within the population may be explained by multiple founders, seedling recruitment during patch-formation, and somatic mutation. The significant genetic differentiation among patches may be caused by separate founding events and/or kin structuring within patches.

Size-class differences in genetic structure and individual distribution of Camellia japonica L. in a Japanese old-growth evergreen forest.

Size-class differences in genetic structure and individual spatial distribution were investigated for Camellia japonica within a 1-ha plot in a Japanese old-growth evergreen forest using microsatellite markers. Three size-classes were considered containing plants that were: 30-32.5 cm tall, 103.8 cm-200 cm tall and those that had a diameter at breast height > or =5 cm, designated JV1, JV2, and ADL, respectively. Each size-class contained 174 individuals. Morisita's index of dispersion indicated clumping of individuals was present within all size-classes, with JV2 displaying the highest level. The clumped distribution of JV1 individuals may be a result of limited seed dispersal, while that of JV2 may be attributed to heterogenieties of favourable microsites, such as canopy gaps. There were no significant differences in allele frequencies among size-classes. There were, however, some differences in spatial genetic structure among them. Moran's I spatial autocorrelation analysis revealed clear spatial genetic structure in class JV1 probably due to limited seed dispersal. In class JV2, genetic structure was not observed. Overlapping seed shadows, probably in canopy gaps, may lead to blurred genetic structure in JV2.

Genetic structure of Camellia japonica L. in an old-growth evergreen forest, Tsushima, Japan.

The spatial genetic structure of Camellia japonica was investigated, using microsatellite markers, in a 4-ha permanent plot within an old-growth forest. Spatial distribution of individuals was also assessed to obtain an insight into spatial relationships between individuals and alleles. Morisita's index of dispersion showed that 518 C. japonica individuals in the plot were clumped, and Moran's I spatial autocorrelation coefficient revealed weak genetic structure, indicating a low level of allele clustering. Average I correlograms showed that there was stronger genetic structure over short-distance classes. The clumped distribution of individuals and the positive autocorrelation over short-distance classes may result from the limited seed dispersal and microsite heterogeneity of the stand, while the genetic structure may be weakened by overlapping seed shadow and extensive pollen flow, mediated by animal vectors, and the high outcrossing rate found in C. japonica.

Genetic diversity and differentiation of Taxodium in the south-eastern United States using cleaved amplified polymorphic sequences.

Two taxa of Taxodium, bald cypress and pond cypress, occur in the south-eastern United States. The ranges of these taxa overlap in the south-eastern Coastal Plain, with the range of the latter being more restricted. Although these taxa co-occur throughout a portion of the more expansive range of bald cypress (Taxodium distichum (L.) L. C. Rich), the habitats of the two taxa appear to differ. Consequently, considerable debate has occurred regarding the taxonomic status of pond cypress. Some authors recognize pond cypress as a distinct species (Taxodium ascendens Brongn.), whereas others recognize it as a variety/ecotype (Taxodium distichum var. imbricarium (Nutt.) Croom). In this study, the genetic diversity of these two taxa was investigated using 10 DNA markers based on sequences from cDNA clones of Cryptomeria japonica. Cryptomeria is a monospecific genus native to Japan, and is a close relative of Taxodium. These markers were codominant in Cryptomeria and were presumed to be codominant in Taxodium. DNA was extracted from leaf tissue collected from six populations of bald cypress and seven populations of pond cypress throughout Florida and Georgia. The average heterozygosities of bald cypress and pond cypress were 0.386 (SE 0.040) and 0.380 (SE 0.040), respectively. Most of the genetic variation (91.9%) was found within populations, 4.9% was found between populations and 3.2% between taxa. Results of DNA analysis using cleaved amplified polymorphic sequences (CAPS) in this study did not suggest that pond cypress was a species distinct from bald cypress. Our conclusion is that the two taxa of Taxodium should be given varietal status.

Intraspecific variation and phylogeographic patterns of Fagus crenata (Fagaceae) mitochondrial DNA.

Mitochondrial (mt) DNA variation in Japanese beech, Fagus crenata (Fagaceae), was studied in 17 populations distributed throughout the species' range. Total genomic DNA of samples from single trees representing each of 12 populations were digested with 18 restriction enzymes and hybridized with three probes containing coxI, coxIII, and atpA gene sequences. Thirty-four of the 54 enzyme/probe combinations showed polymorphisms and all the individuals were subsequently analyzed with six combinations of three probes and two enzymes. Restriction fragment length polymorphisms were evident around all three genes, allowing the identification of eight distinct haplotypes. Haplotype diversity within the populations was found to be very low (HS = 0.031), but population differentiation to be much higher (GST = 0.963). The mtDNA variation was strikingly different from allozyme variation (HS = 0.209; GST = 0.039). Gene flow for maternally inherited mtDNA should be restricted to seed dispersal while nuclear gene flow occurs by both seed and pollen dispersal. Therefore, the difference in the variation between mtDNA and allozymes may be largely a result of the much higher rate of gene flow associated with pollen dispersal than with seed dispersal. The mtDNA variation displayed strong geographic structure, which may reflect the species' distribution in the last glacial maximum and subsequent colonization, and probably also reflects intraspecific phylogeography of the species.

Molecular phytogeny of conifers using RFLP analysis of PCR-amplified specific chloroplast genes.

We investigated the molecular phylogeny of conifers using restriction endonuclease fragment length polymorphism of six polymerase chain reaction-amplified chloroplast genes - frxC, rbcL, psbA, psbD, trnK, and 16S. We detected 227 total site changes among species, representing 23, 26, 38, 48, 67, and 25 site changes in frxC, psbA, psbD, rbcL, trnK and 16S, respectively. The mean nucleotide substitution was 10.75% (SD 0.573) among species in five families. Forty maximally parsimonious trees were obtained using the Wagner parsimony method, and a 50% majority-rule consensus tree was obtained from them. Data analysis produced similar basic patterns when both the Wagner parsimony and the neighbor-joining methods were applied, and the main lineages were clearly separated. Taxaceae and Cephalotaxaceae species were used as the out-groups when applying Wagner parsimony methods. With the Wagner method, the consistency index was 0.510, the retention index was 0.879, and tree length was 435 steps. Our results indicated that Cupressaceae and Taxodiaceae are closely related families and that Sciadopitys verticillata is the basal lineage of Cupressaceae and Taxodiaceae. The neighbor-joining tree is similar to the 50% majority-rule consensus of the 40 Wagner parsimony trees except for the position of Keteleeria daversifolia, the Picea and Cedrus group, and the divergence within Cupressaceae.

A linkage map for sugi (Cryptomeria japonica) based on RFLP, RAPD, and isozyme loci.

A linkage map for sugi was constructed on the basis of restriction fragment length polymorphism (RFLP), random amplified polymorphic DNA (RAPD), and isozyme loci using a three-generation pedigree prepared for genetic analysis of heartwood color. A total of 128 RFLP (123 cDNA and 5 genomic probes), 33 RAPD, 2 isozyme, and 1 morphological (dwarf) loci segregated in 73 progeny. Of the 164 segregating loci, 145 loci were distributed in 20 linkage groups. Of these loci, 91 with confirmed map positions were assigned to 13 linkage groups, covering a total of 887.3 cM. A clustering of markers with distorted segregation was observed in 6 linkage groups. In the four clusters, distortions with a reduction in the number of homozygotes from one parent only were found.