UBC18 mediates ERF1 degradation under light-dark cycles.

Ethylene Response Factor 1 (ERF1) plays a crucial role in biotic and abiotic stress responses. Previous studies have shown that ERF1 regulates stress-responsive gene expression by binding to different cis-acting elements in response to various stress signals. ERF1 was also reported to be unstable in the dark, and it regulates hypocotyl elongation. Here, we elucidated the mechanism underlying degradation of ERF1. Yeast two-hybrid screening showed that UBIQUITIN-CONJUGATING ENZYME 18 (UBC18) interacted with ERF1. The interaction between ERF1 and UBC18 was verified using pull-down assays and coimmunoprecipitation analyses. We then compared the ERF1 protein abundance in the UBC18 mutant and overexpression plants. Based on the results of protein degradation and in vivo ubiquitination assays, we proposed that UBC18 mediates ERF1 ubiquitination and degradation. ERF1 was more stable in UBC18 mutants and less stable in UBC18 overexpression lines compared with that in wild-type plants. ERF1 was degraded by the 26S proteasome system via regulation of UBC18 and promotes dark-repression of downstream genes and proline accumulation. UBC18 negatively regulated drought and salt stress responses by altering the abundance of ERF1 and the expression of genes downstream of ERF1.

Increased glutathione contributes to stress tolerance and global translational changes in Arabidopsis.

Although glutathione is well known for its reactive oxygen species (ROS) scavenging function and plays a protective role in biotic stress, its regulatory function in abiotic stress still remains to be elucidated. Our previous study showed that exogenously applied reduced glutathione (GSH) could improve abiotic stress tolerance in Arabidopsis. Here, we report that endogenously increased GSH also conferred tolerance to drought and salt stress in Arabidopsis. Moreover, both exogenous and endogenous GSH delayed senescence and flowering time. Polysomal profiling results showed that global translation was enhanced after GSH treatment and by the induced increase of GSH level by salt stress. By performing transcriptomic analyses of steady-state and polysome-bound mRNAs in GSH-treated plants, we reveal that GSH has a substantial impact on translation. Translational changes induced by GSH treatment target numerous hormones and stress signaling molecules, which might contribute to the enhanced stress tolerance in GSH-treated plants. Our translatome analysis also revealed that abscisic acid (ABA), auxin and jasmonic acid (JA) biosynthesis, as well as signaling genes, were activated during GSH treatment, which has not been reported in previously published transcriptomic data. Together, our data suggest that the increased glutathione level results in stress tolerance and global translational changes.

ORA47 (octadecanoid-responsive AP2/ERF-domain transcription factor 47) regulates jasmonic acid and abscisic acid biosynthesis and signaling through binding to a novel cis-element.

ORA47 (octadecanoid-responsive AP2/ERF-domain transcription factor 47) of Arabidopsis thaliana is an AP2/ERF domain transcription factor that regulates jasmonate (JA) biosynthesis and is induced by methyl JA treatment. The regulatory mechanism of ORA47 remains unclear. ORA47 is shown to bind to the cis-element (NC/GT)CGNCCA, which is referred to as the O-box, in the promoter of ABI2. We proposed that ORA47 acts as a connection between ABA INSENSITIVE1 (ABI1) and ABI2 and mediates an ABI1-ORA47-ABI2 positive feedback loop. PORA47:ORA47-GFP transgenic plants were used in a chromatin immunoprecipitation (ChIP) assay to show that ORA47 participates in the biosynthesis and/or signaling pathways of nine phytohormones. Specifically, many abscisic acid (ABA) and JA biosynthesis and signaling genes were direct targets of ORA47 under stress conditions. The JA content of the P35S:ORA47-GR lines was highly induced under wounding and moderately induced under water stress relative to that of the wild-type plants. The wounding treatment moderately increased ABA accumulation in the transgenic lines, whereas the water stress treatment repressed the ABA content. ORA47 is proposed to play a role in the biosynthesis of JA and ABA and in regulating the biosynthesis and/or signaling of a suite of phytohormone genes when plants are subjected to wounding and water stress.

The Arabidopsis ETHYLENE RESPONSE FACTOR1 regulates abiotic stress-responsive gene expression by binding to different cis-acting elements in response to different stress signals.

ETHYLENE RESPONSE FACTOR1 (ERF1) is an upstream component in both jasmonate (JA) and ethylene (ET) signaling and is involved in pathogen resistance. Accumulating evidence suggests that ERF1 might be related to the salt stress response through ethylene signaling. However, the specific role of ERF1 in abiotic stress and the molecular mechanism underlying the signaling cross talk still need to be elucidated. Here, we report that ERF1 was highly induced by high salinity and drought stress in Arabidopsis (Arabidopsis thaliana). The salt stress induction required both JA and ET signaling but was inhibited by abscisic acid. ERF1-overexpressing lines (35S:ERF1) were more tolerant to drought and salt stress. They also displayed constitutively smaller stomatal aperture and less transpirational water loss. Surprisingly, 35S:ERF1 also showed enhanced heat tolerance and up-regulation of heat tolerance genes compared with the wild type. Several suites of genes activated by JA, drought, salt, and heat were found in microarray analysis of 35S:ERF1. Chromatin immunoprecipitation assays found that ERF1 up-regulates specific suites of genes in response to different abiotic stresses by stress-specific binding to GCC or DRE/CRT. In response to biotic stress, ERF1 bound to GCC boxes but not DRE elements; conversely, under abiotic stress, we observed specific binding of ERF1 to DRE elements. Furthermore, ERF1 bound preferentially to only one among several GCC box or DRE/CRT elements in the promoter region of its target genes. ERF1 plays a positive role in salt, drought, and heat stress tolerance by stress-specific gene regulation, which integrates JA, ET, and abscisic acid signals.

Functional characterization of an abiotic stress-inducible transcription factor AtERF53 in Arabidopsis thaliana.

AP2/ERF proteins play crucial roles in plant growth and development and in responses to biotic and abiotic stresses. ETHYLENE RESPONSE FACTOR 53 (AtERF53) belongs to group 1 in the ERF family and is induced in the early hours of dehydration and salt treatment. The functional study of AtERF53 is hampered because its protein expression in Arabidopsis is vulnerable to degradation in overexpressed transgenic lines. Taking advantage of the RING domain ligase1/RING domain ligase2 (rglg1rglg2) double mutant in which the AtERF53 can express stably, we investigate the physiological function of AtERF53. In this study, we demonstrate that expression of AtERF53 in wild-type Arabidopsis was responsive to heat and abscisic acid (ABA) treatment. From results of the cotransfection experiment, we concluded that AtERF53 has positive transactivation activity. Overexpression of AtERF53 in the rglg1rglg2 double mutant conferred better heat-stress tolerance and had resulted in higher endogenous ABA and proline levels compared to rglg1rglg2 double mutants. AtERF53 also has a function to regulate guard-cell movement because the stomatal aperture of AtERF53 overexpressed in rglg1rglg2 double mutant was smaller than that in the rglg1rglg2 double mutant under ABA treatment. In a global gene expression study, we found higher expressions of many stress-related genes, such as DREB1A, COR15A, COR15B, PLC, P5CS1, cpHSC70 s and proline and ABA metabolic-related genes. Furthermore, we identified several downstream target genes of AtERF53 by chromatin immunoprecipitation assay. In conclusion, the genetic, molecular and biochemical result might explain how AtERF53 serving as a transcription factor contributes to abiotic stress tolerance in Arabidopsis.

Arabidopsis RGLG2, functioning as a RING E3 ligase, interacts with AtERF53 and negatively regulates the plant drought stress response.

Transcriptional activities of plants play important roles in responses to environmental stresses. ETHYLENE RESPONSE FACTOR53 (AtERF53) is a drought-induced transcription factor that belongs to the AP2/ERF superfamily and has a highly conserved AP2 domain. It can regulate drought-responsive gene expression by binding to the GCC box and/or the dehydration-responsive element in the promoter of downstream genes. Overexpression of AtERF53 driven by the cauliflower mosaic virus 35S promoter resulted in an unstable drought-tolerant phenotype in T2 transgenic Arabidopsis (Arabidopsis thaliana) plants. Using a yeast two-hybrid screen, we identified a RING domain ubiquitin E3 ligase, RGLG2, which interacts with AtERF53 in the nucleus. The copine domain of RGLG2 exhibited the strongest interacting activity. We also demonstrated that RGLG2 could move from the plasma membrane to the nucleus under stress treatment. Using an in vitro ubiquitination assay, RGLG2 and its closest sequelog, RGLG1, were shown to have E3 ligase activity and mediated AtERF53 ubiquitination for proteasome degradation. The rglg1rglg2 double mutant but not the rglg2 or rglg1 single mutant exhibited a drought-tolerant phenotype when compared with wild-type plants. AtERF53-green fluorescent proteins expressed in the rglg1rglg2 double mutants were stable. The 35S:AtERF53-green fluorescent protein/rglg1rglg2 showed enhanced AtERF53-regulated gene expression and had greater tolerance to drought stress than the rglg1rglg2 double mutant. In conclusion, RGLG2 negatively regulates the drought stress response by mediating AtERF53 transcriptional activity in Arabidopsis.

Drought and salt stress tolerance of an Arabidopsis glutathione S-transferase U17 knockout mutant are attributed to the combined effect of glutathione and abscisic acid.

Although glutathione S-transferases (GSTs) are thought to play major roles in oxidative stress metabolism, little is known about the regulatory functions of GSTs. We have reported that Arabidopsis (Arabidopsis thaliana) GLUTATHIONE S-TRANSFERASE U17 (AtGSTU17; At1g10370) participates in light signaling and might modulate various aspects of development by affecting glutathione (GSH) pools via a coordinated regulation with phytochrome A. Here, we provide further evidence to support a negative role of AtGSTU17 in drought and salt stress tolerance. When AtGSTU17 was mutated, plants were more tolerant to drought and salt stresses compared with wild-type plants. In addition, atgstu17 accumulated higher levels of GSH and abscisic acid (ABA) and exhibited hyposensitivity to ABA during seed germination, smaller stomatal apertures, a lower transpiration rate, better development of primary and lateral root systems, and longer vegetative growth. To explore how atgstu17 accumulated higher ABA content, we grew wild-type plants in the solution containing GSH and found that they accumulated ABA to a higher extent than plants grown in the absence of GSH, and they also exhibited the atgstu17 phenotypes. Wild-type plants treated with GSH also demonstrated more tolerance to drought and salt stresses. Furthermore, the effect of GSH on root patterning and drought tolerance was confirmed by growing the atgstu17 in solution containing l-buthionine-(S,R)-sulfoximine, a specific inhibitor of GSH biosynthesis. In conclusion, the atgstu17 phenotype can be explained by the combined effect of GSH and ABA. We propose a role of AtGSTU17 in adaptive responses to drought and salt stresses by functioning as a negative component of stress-mediated signal transduction pathways.

Structural and functional assays of AtTLP18.3 identify its novel acid phosphatase activity in thylakoid lumen.

The membrane protein AtTLP18.3 of Arabidopsis (Arabidopsis thaliana) contains a domain of unknown function, DUF477; it forms a polysome with photosynthetic apparatuses in the thylakoid lumen. To explore the molecular function of AtTLP18.3, we resolved its crystal structures with residues 83 to 260, the DUF477 only, and performed a series of biochemical analyses to discover its function. The gene expression of AtTLP18.3 followed a circadian rhythm. X-ray crystallography revealed the folding of AtTLP18.3 as a three-layer sandwich with three α-helices in the upper layer, four ß-sheets in the middle layer, and two α-helices in the lower layer, which resembles a Rossmann fold. Structural comparison suggested that AtTLP18.3 might be a phosphatase. The enzymatic activity of AtTLP18.3 was further confirmed by phosphatase assay with various substrates (e.g. p-nitrophenyl phosphate, 6,8-difluoro-4-methylumbelliferyl phosphate, O-phospho-L-serine, and several synthetic phosphopeptides). Furthermore, we obtained the structure of AtTLP18.3 in complex with O-phospho-L-serine to identify the binding site of AtTLP18.3. Our structural and biochemical studies revealed that AtTLP18.3 has the molecular function of a novel acid phosphatase in the thylakoid lumen. DUF477 is accordingly renamed the thylakoid acid phosphatase domain.

Historical spatial range expansion and a very recent bottleneck of Cinnamomum kanehirae Hay. (Lauraceae) in Taiwan inferred from nuclear genes.

BACKGROUND: Species in the varied geographic topology of Taiwan underwent obvious demographic changes during glacial periods. Cinnamomum kanehirae has been exploited for timber and to obtain medicinal fungi for the past 100 years. Understanding anthropogenic factors influencing the demography of this species after the last glacial maximum (LGM) is critically important for the conservation of this species. RESULTS: Populations of C. kanehirae were classified into four geographic regions: northwestern (NW), west-central (WC), southwestern (SW), and southeastern (SE). In total, 113 individuals from 19 localities were sampled, and variations in the chalcone synthase gene (Chs) intron and leafy (Lfy) intron-2 sequences of nuclear DNA were examined in order to assess phylogeographic patterns, the timescales of demographic and evolutionary events, and recent anthropogenic effects. In total, 210 Chs and 170 Lfy sequences, which respectively constituted 36 and 35 haplotypes, were used for the analyses. Estimates of the migration rate (M) through time revealed a pattern of frequent gene flow during previous and the present interglacials. The isolation-by-distance test showed that there generally was no significant correlation between genetic and geographic distances. The level of among-region genetic differentiation was significant when comparing eastern to western populations. However, no significant among-region genetic differentiation was found in comparisons among the four geographic regions. Moreover, essentially no genetic structuring was found for the three regions west of the CMR. A fit of spatial range expansion was found for pooled and regional samples according to the non-significant values of the sum of squared deviations. Using the Bayesian skyline plot (BSP) method, a recent bottleneck after the LGM expansion was detected in both regional and pooled samples. CONCLUSIONS: Common haplotype distributions among geographic regions and the relatively shallow genetic structuring displayed are the result of historical gene flows. Southward dispersals in an earlier time frame from the NW region and in a later time frame from the SE region were inferred. The BSP analysis suggested a postglacial expansion event. Recent trends, however, refer to a bottleneck due to human interventions observed for both pooled and regional C. kanehirae samples.

Duplication of the class I cytosolic small heat shock protein gene and potential functional divergence revealed by sequence variations flanking the {alpha}-crystallin domain in the genus Rhododendron (Ericaceae).

BACKGROUND AND AIMS: Positive selection in the -crystallin domain (ACD) of the chloroplast small heat shock protein (CPsHSP) gene was found in a previous study and was suggested to be related to the ecological adaptation of Rhododendron species in the subgenus Hymenanthes. Consequently, it was of interest to examine whether gene duplication and subsequent divergence have occurred in other sHSP genes, for example class I cytosolic sHSP genes (CT1sHSPs) in Rhododendron in Taiwan, where many endemic species have evolved as a result of habitat differentiation. METHODS: A phylogeny of CT1sHSP amino acid sequences was built from Rhododendron, Arabidopsis thaliana, Oryza sativa, Populus trichocarpa, Vitis vinifera and other species for elucidation of the phylogenetic relationships among CT1sHSPs. Phylogenies of Rhododendron CT1sHSP nucleotide and amino acid sequences were generated for positive selection and functional divergence analysis, respectively. Positively selected sites and amino acid differences between types of Rhododendron CT1sHSPs were mapped onto the wheat sHSP16.9 protein structure. Average genetic distance (Dxy) and dN/dS ratios between types of Rhododendron CT1sHSP genes were analysed using sliding window analysis. Gene conversion was also assessed between types of Rhododendron CT1sHSPs. KEY RESULTS: Two types of Rhododendron CT1sHSP were identified. A high level of genetic similarity and diversity within and flanking the ACD, respectively, between types of Rhododendron CT1sHSP were found. Main differences between the two types of Rhododendron CT1sHSPs were: (1) increased hydrophobicity by two positively selected amino acid sites and a seven-amino-acid insertion in the N-terminal arm; and (2) increased structural flexibility and solubility by a seven-amino-acid insertion in the N-terminal arm and one positively selected amino acid site in the C-terminal extension. CONCLUSIONS: Functional conservation of the ACD of Rhododendron CT1sHSP genes was inferred b

Molecular evolution of the Pi-ta gene resistant to rice blast in wild rice (Oryza rufipogon).

Rice blast disease resistance to the fungal pathogen Magnaporthe grisea is triggered by a physical interaction between the protein products of the host R (resistance) gene, Pi-ta, and the pathogen Avr (avirulence) gene, AVR-pita. The genotype variation and resistant/susceptible phenotype at the Pi-ta locus of wild rice (Oryza rufipogon), the ancestor of cultivated rice (O. sativa), was surveyed in 36 locations worldwide to study the molecular evolution and functional adaptation of the Pi-ta gene. The low nucleotide polymorphism of the Pi-ta gene of O. rufipogon was similar to that of O. sativa, but greatly differed from what has been reported for other O. rufipogon genes. The haplotypes can be subdivided into two divergent haplogroups named H1 and H2. H1 is derived from H2, with nearly no variation and at a low frequency. H2 is common and is the ancestral form. The leucine-rich repeat (LRR) domain has a high pi(non)/pi(syn) ratio, and the low polymorphism of the Pi-ta gene might have primarily been caused by recurrent selective sweep and constraint by other putative physiological functions. Meanwhile, we provide data to show that the amino acid Ala-918 of H1 in the LRR domain has a close relationship with the resistant phenotype. H1 might have recently arisen during rice domestication and may be associated with the scenario of a blast pathogen-host shift from Italian millet to rice.

Constitutive components and induced gene expression are involved in the desiccation tolerance of Selaginella tamariscina.

Selaginella tamariscina, one of the most primitive vascular plants, can remain alive in a desiccated state and resurrect when water becomes available. To evaluate the nature of desiccation tolerance in this plant, we compared the composition of soluble sugars and saturation ratios of phospholipids (PLs) between hydrated and desiccated tissues of S. tamariscina using gas chromatography. In this study, differences in gene expression and ABA contents were also analyzed during dehydration. The results revealed that trehalose (at >130 mg g(-1) DW) was the major soluble sugar, and low saturated fatty acid content in PLs (0.31) was maintained in both hydrated and desiccated tissues. In addition, the ABA content of S. tamariscina increased 3-fold, and genes involved in ABA signaling and cellular protection were up-regulated while photosystem-related genes were down-regulated during dehydration. The biochemical and molecular findings suggest that both constitutive and inducible protective molecules contribute to desiccation tolerance of S. tamariscina.

Molecular population genetics and gene expression analysis of duplicated CBF genes of Arabidopsis thaliana.

BACKGROUND: CBF/DREB duplicate genes are widely distributed in higher plants and encode transcriptional factors, or CBFs, which bind a DNA regulatory element and impart responsiveness to low temperatures and dehydration. RESULTS: We explored patterns of genetic variations of CBF1, -2, and -3 from 34 accessions of Arabidopsis thaliana. Molecular population genetic analyses of these genes indicated that CBF2 has much reduced nucleotide diversity in the transcriptional unit and promoter, suggesting that CBF2 has been subjected to a recent adaptive sweep, which agrees with reports of a regulatory protein of CBF2. Investigating the ratios of Ka/Ks between all paired CBF paralogus genes, high conservation of the AP2 domain was observed, and the major divergence of proteins was the result of relaxation in two regions within the transcriptional activation domain which was under positive selection after CBF duplication. With respect to the level of CBF gene expression, several mutated nucleotides in the promoters of CBF3 and -1 of specific ecotypes might be responsible for its consistently low expression. CONCLUSION: We concluded from our data that important evolutionary changes in CBF1, -2, and -3 may have primarily occurred at the level of gene regulation as well as in protein function.

Uniform genetic diversity, low differentiation, and neutral evolution characterize contemporary refuge populations of Taiwan fir (Abies kawakamii, Pinaceae).

Based on fossil pollen, the distribution range of Taiwan fir [Abies kawakamii (Hay.) Ito] (Pinaceae) is smaller than it was 50 000 years ago. To characterize the present refuge populations of A. kawakamii, which survive only in subalpine forests in Taiwan, we surveyed nuclear genes and chloroplast intergenic spacers to assess the genetic diversity of Taiwan fir. Populations maintain high genetic diversity and contain similar numbers of haplotypes for the GapC (cytosolic glyceraldehyde 3-phosphate dehydrogenase) fragment. Haplotypes for GapC are generally widespread, and population-specific haplotypes accounted for 2.5% of the total. Differentiation among populations is very low (G(ST) = 0.01). Only three haplotypes were detected for the cpDNA marker, and every population had one or two haplotypes. In a neutrality test, the variation in nucleotides did not deviate from that expected with neutral evolution for either marker. A retreat route to higher elevations was not evident from either the GapC or cpDNA markers. Hsuehshan was the site of the most divergent population in Taiwan. We concluded that uniform genetic diversity, low differentiation, low numbers of population-specific haplotypes, and neutral evolution characterize contemporary refuge populations of Taiwan fir.

Divergent evolution of the chloroplast small heat shock protein gene in the genera Rhododendron (Ericaceae) and Machilus (Lauraceae).

BACKGROUND AND AIMS: Evolutionary and ecological roles of the chloroplast small heat shock protein (CPsHSP) have been emphasized based on variations in protein contents; however, DNA sequence variations related to the evolutionary and ecological roles of this gene have not been investigated. In the present study, a basal angiosperm, Machilus, together with the eudicot Rhododendron were used to illustrate the evolutionary dynamics of gene divergence in CPsHSPs. METHODS: Degenerate primers were used to amplify CPsHSP-related sequences from 16 Rhododendron and eight Machilus species that occur in Taiwan. Manual DNA sequence alignment was carried out according to the deduced amino acid sequence alignment performed by CLUSTAL X. A neighbour-joining tree was generated in MEGA using conceptual translated amino acid sequences from consensus sequences of cloned CPsHSP genes from eight Machilus and 16 Rhododendron species as well as amino acid sequences of CPsHSPs from five monocots and seven other eudicots acquired from GenBank. CPsHSP amino acid sequences of Funaria hygrometrica were used as the outgroups. The aligned DNA and amino acid sequences were used to estimate several parameters of sequence divergence using the MEGA program. Separate Bayesian inference of DNA sequences of Rhododendron and Machilus species was analysed and the resulting gene trees were used for detection of putative positively selected amino acid sites by the Codeml program implemented in the PAML package. Mean hydrophobicity profile analysis was performed with representative amino acid sequences for both Rhododendron and Machilus species by the Bioedit program. The computer program SplitTester was used to examine whether CPsHSPs of Rhododendron lineages and duplicate copies of the Machilus CPsHSPs have evolved functional divergence based on the hydrophobicity distance matrix. KEY RESULTS: Only one copy of the CPsHSP was found in Rhododendron. However, a higher evolutionary rate of amino acid substitutions in the Hymenanthes lineage of Rhododendron was inferred. Two positively selected amino acid sites may have resulted in higher hydrophobicity in the region of the alpha-crystallin domain (ACD) of the CPsHSP. By contrast, the basal angiosperm, Machilus, possessed duplicate copies of the CPsHSP, which also differed in their evolutionary rates of amino acid substitutions. However, no apparent relationship of ecological relevance toward the positively selected amino acid sites was found in Machilus. CONCLUSIONS: Divergent evolution was found for both Rhododendron lineages and the paralogues of CPsHSP in Machilus that were directed to the shift in hydrophobicity in the ACD and/or methionine-rich region, which might have played important roles in molecular chaperone activity.

Phylogeographic study reveals the origin and evolutionary history of a Rhododendron species complex in Taiwan.

This study infers a single origin and a once-widespread distribution of the Rhododendron pseudochrysanthum species complex in Taiwan based on chloroplast DNA (cpDNA) variation. In total, 124 individuals from five endemic Rhododendron species were used for amplifications of two chloroplast intergenic spacers: trnL-trnF and atpB-rbcL. The haplotype and nucleotide diversities were much lower for the R. pseudochrysanthum complex, comprised of the species R. pseudochrysanthum, R. morii, Rhododendron rubropunctatum, and Rhododendron hyperythrum, than for Rhododendron formosanum. Two measures of pairwise population differentiation, N(ST) and F(ST), consistently revealed mostly non-significant levels of genetic divergence between populations of the R. pseudochrysanthum complex. No genetic difference was found among the four species of the R. pseudochrysanthum complex by analysis of molecular variance (AMOVA), which is concordant with the parsimonious topology of cpDNA haplotypes for the complex. Nested clade analysis (NCA) of the cpDNA haplotypes indicated that restricted gene flow with isolation-by-distance characterized the recolonization pattern of the R. pseudochrysanthum complex. In contrast, the NCA analysis indicated a contiguous range expansion for cpDNA haplotypes of R. formosanum. This research suggests a once-widespread distribution of the R. pseudochrysanthum complex probably via north-to-south colonization of mid-elevations during low-temperature periods of the Pleistocene. Population fragmentation followed the warmer climate which began in the Holocene and resulted in the present-day range contraction into high elevations.

Allozyme variation of populations of Castanopsis carlesii (Fagaceae) revealing the diversity centres and areas of the greatest divergence in Taiwan.

BACKGROUND AND AIMS: The genetic variation and divergence estimated by allozyme analysis were used to reveal the evolutionary history of Castanopsis carlesii in Taiwan. Two major questions were discussed concerning evolutionary issues: where are the diversity centres, and where are the most genetically divergent sites in Taiwan? METHODS: Twenty-two populations of C. carlesii were sampled throughout Taiwan. Starch gel electrophoresis was used to assay allozyme variation. Genetic parameters and mean FST values of each population were analysed using the BIOSYS-2 program. Mean F(ST) values of each population against the remaining populations, considered as genetic divergence, were estimated using the FSTAT program. KEY RESULTS: Average values of genetic parameters describing the within-population variation, the average number of alleles per locus (A=2.5), the effective number of alleles per locus (Ae=1.38), the allelic richness (Ar=2.38), the percentage of polymorphic loci (P=69%), and the expected heterozygosity (He=0.270) were estimated. High levels of genetic diversity were found for C. carlesii compared with other local plant species. Genetic differentiation between populations was generally low. CONCLUSIONS: From the data of expected heterozygosity, one major diversity centre was situated in central Taiwan corroborating previous reports for other plant species. According to the mean FST value of each population, the most divergent populations were situated in two places. One includes populations located in north central Taiwan between 24.80 degrees N and 24.20 degrees N. The other is located in south-eastern Taiwan between 22.40 degrees N and 23.10 degrees N. These two regions are approximately convergent with the most divergent locations determined for several other plant species using chloroplast DNA markers published previously. An important finding obtained from this study is that unordered markers like allozymes can be used to infer past population histories as well as chloroplast DNA markers do.

Source populations of Quercus glauca in the last glacial age in Taiwan revealed by nuclear microsatellite markers.

In this work, we attempted to study genetic differentiation between populations of Quercus glauca in Taiwan using nuclear microsatellite markers to infer the potential refugium in the last glaciation stage. Four microsatellite loci for 20 individuals each in 10 populations of Taiwan were analyzed. We found that Q. glauca has relatively high within-population diversity (H(E) = 0.741) and low population differentiation (F(ST) = 0.042) but shows isolation by distance. The most divergent populations, according to the average F(ST) for individual populations in comparison with every other population, were found in populations Cy, Sa, and Hy in southern Taiwan and Pa in north-central Taiwan. Moreover, populations Cy, Sa, and Pa were recognized as being the source populations for gene recolonization after the last glaciation stage. In addition, the three sites of Wu, Ym, and Cy exhibited the highest gene diversities that coincided with populations with the highest chloroplast DNA variations. This may have resulted from an admixture of colonization routes. In conclusion, observations of the most divergent populations and source populations suggest that southern and probably north-central Taiwan may have potentially been refugia for Q. glauca in the last glaciation. This agrees with the possible refugium in southern Taiwan revealed by a previous study using chloroplast DNA markers.

Potential refugia in Taiwan revealed by the phylogeographical study of Castanopsis carlesii Hayata (Fagaceae).

In this study, we examined spatial patterns of chloroplast DNA (cpDNA) variation in a total of 30 populations of Castanopsis carlesii Hayata (Fagaceae), a subtropical and temperate tree species, including 201 individuals sampled throughout Taiwan. By sequencing two cpDNA fragments using universal primers (the trnL intron and the trnV-trnM intergenic spacer), we found a total of 1663 bp and 21 polymorphic sites. These gave rise to a total of 28 cpDNA haplotypes. The level of differentiation among the populations studied was relatively high (GST = 0.723). Two ancestral haplotypes are widely distributed. The Central Mountain Ridge (CMR) of Taiwan represents an insurmountable barrier to the east-west gene flow of C. carlesii. Among the populations studied, three separated populations, at Lienhuachih, Fushan and Lichia, have high nucleotide diversity. Estimates of NST-GST for populations on both sides of the CMR indicate that no phylogeographical structure exists. According to the genealogical tree, number of rare haplotype and population genetic divergence, this study suggests that two potential refugia existed during the last glaciation: the first refugium was located in a region to the north of Hsuehshan Range (HR) and west of the CMR; the second refugium was located in south, especially southeastern Taiwan. In fact, the second refugium happens to be the same as that reported for Quercus glauca. A 'star-like' genealogy is characteristic when all haplotypes rapidly coalesce and is a general outcome of population expansion. The neutrality test and mismatch distribution also suggest demographic expansion recovering from a bottleneck.

Spatial pattern of chloroplast DNA variation of Cyclobalanopsis glauca in Taiwan and East Asia.

This study examined the spatial pattern of chloroplast DNA (cpDNA) variation in Cyclobalanopsis glauca (Thunb. ex Murray) Oerst. (Fagaceae) in 140 trees from Taiwan (25 populations), Japan (three), Ryukyus (two), Hong Kong (one) and Mainland China (one). By sequencing three cpDNA intergenic spacer fragments using universal primers (trnT-trnL, trnV-trnM, including the trnV intron, and petG-trnP), we found a total of 1,980 bp and 15 polymorphic sites. Among them, 12 sites were caused by point mutation, and three resulted from insertion. This gives rise to a total of 13 cpDNA haplotypes. The level of differentiation among the populations studied is relatively high (GST = 0.612). Two ancestral haplotypes (A and B) are distributed widely in East Asia. Interestingly, all the derived cpDNA variations are found only in Taiwan but not in other areas. The Central Mountain Ridge (CMR) of Taiwan creates an unsurpassed barrier to the east-west gene flow of C. glauca. Among the populations on the west of CMR, only three separated populations, Yangmingshan, Wushe and Chinshuiying, have high haplotype diversity, each consisting of sister haplotypes all mutated from the same ancestral haplotype. Thus, they have probably originated from de novo mutation after the last glaciation. This inference agrees with the observation that no spatial autocorrelation existed on the west side. Two unrelated dominant lineages on the east of the CMR (haplotypes D and F) showed significant spatial genetic structure. Estimate of NST - GST was -0.090 and differed significantly from zero. Thus at the local scale, the phylogeographical component of the genetic structure is significant on the east of the CMR. Accompanied by published palynological records of the last glaciation, this study suggests the possibility that these two types were colonized northward from the southeastern part of Taiwan. 'Star-like' genealogy is characterized, with all the haplotypes coalescing rapidly and as a general outcome of population expansion (Page & Holmes 1998). A neutrality test also suggested a demographic expansion recovered from a bottleneck. We therefore inferred that the southeastern part of Taiwan might be a potential refugium for C. glauca.