Search

We epigenotyped 211 individuals from 17 Zingiber kawagoii populations using methylation-sensitive amplification polymorphism (MSAP) and investigated the associations of methylated (mMSAP) and unmethylated (uMSAP) loci with 16 environmental variables. Data regarding genetic variation based on amplified fragment length polymorphism (AFLP) were obtained from an earlier study. We found a significant positive correlation between genetic and epigenetic variation. Significantly higher mean mMSAP and uMSAP uHE (unbiased expected heterozygosity: 0.223 and 0.131, respectively, p < 0.001) per locus than that estimated based on AFLP (uHE = 0.104) were found. Genome scans detected 10 mMSAP and 9 uMSAP FST outliers associated with various environmental variables. A significant linear fit for 11 and 12 environmental variables with outlier mMSAP and uMSAP ordination, respectively, generated using full model redundancy analysis (RDA) was found. When conditioned on geography, partial RDA revealed that five and six environmental variables, respectively, were the most important variables influencing outlier mMSAP and uMSAP variation. We found higher genetic (average FST = 0.298) than epigenetic (mMSAP and uMSAP average FST = 0.044 and 0.106, respectively) differentiation and higher genetic isolation-by-distance (IBD) than epigenetic IBD. Strong epigenetic isolation-by-environment (IBE) was found, particularly based on the outlier data, controlling either for geography (mMSAP and uMSAP βE = 0.128 and 0.132, respectively, p = 0.001) or for genetic structure (mMSAP and uMSAP βE = 0.105 and 0.136, respectively, p = 0.001). Our results suggest that epigenetic variants can be substrates for natural selection linked to environmental variables and complement genetic changes in the adaptive evolution of Z. kawagoii populations.

Ecological and evolutionary processes linking adaptation to environment are related to species’ range shifts. In this study, we employed amplified-fragment-length-polymorphism-based genome scan methods to identify candidate loci among Zingiber kawagoii populations inhabiting varying environments distributed at low to middle elevations (143–1488 m) in a narrow latitudinal range (between 21.90 and 25.30° N). Here, we show evidence of selection driving the divergence of Z. kawagoii. Twenty-six FST outliers were detected, which were significantly correlated with various environmental variables. The allele frequencies of nine FST outliers were either positively or negatively correlated with the population mean FST. Using several independent approaches, we found environmental variables act in a combinatorial fashion, best explaining outlier genetic variation. Nonetheless, we found that adaptive divergence was affected mostly by annual temperature range, and it is significantly positively correlated with latitude and significantly negatively correlated with the population mean FST. This study addresses a latitudinal pattern of changes in annual temperature range (which ranged from 13.8 °C in the Lanyu population to 18.5 °C in the Wulai population) and emphasizes the pattern of latitudinal population divergence closely linked to the allele frequencies of adaptive loci, acting in a narrow latitudinal range. Our results also indicate environmentally dependent local adaptation for both leading- and trailing-edge populations.

The testing association of environmental variables with genetic and epigenetic variation could be crucial to deciphering the effects of environmental factors playing roles as selective drivers in ecological speciation. Although ecological speciation may occur in closely related species, species boundaries may not be established over a short evolutionary timescale. Here, we investigated the genetic and epigenetic variations using amplified fragment length polymorphism (AFLP) and methylation-sensitive amplification polymorphism (MSAP), respectively, and tested their associations with environmental variables in populations of four closely related species in the R. pseudochrysanthum complex. No distinctive species relationships were found using genetic clustering analyses, neighbor-joining tree, and neighbor-net tree based on the total AFLP variation, which is suggestive of the incomplete lineage sorting of ancestral variation. Nonetheless, strong isolation-by-environment and adaptive divergence were revealed, despite the significant isolation-by-distance. Annual mean temperature, elevation, normalized difference vegetation index, and annual total potential evapotranspiration were found to be the most important environmental variables explaining outlier genetic and epigenetic variations. Our results suggest that the four closely related species of the R. pseudochrysanthum complex share the polymorphism of their ancestor, but reproductive isolation due to ecological speciation can occur if local environmental divergence persists over time.

Population diversification can be shaped by a combination of environmental factors as well as geographic isolation interacting with gene flow. We surveyed genetic variation of 243 samples from 12 populations of Calocedrus formosana using amplified fragment length polymorphism (AFLP) and scored a total of 437 AFLP fragments using 11 selective amplification primer pairs. The AFLP variation was used to assess the role of gene flow on the pattern of genetic diversity and to test environments in driving population adaptive evolution. This study found the relatively lower level of genetic diversity and the higher level of population differentiation in C. formosana compared with those estimated in previous studies of conifers including Cunninghamia konishii, Keteleeria davidiana var. formosana, and Taiwania cryptomerioides occurring in Taiwan. BAYESCAN detected 26 FST outlier loci that were found to be associated strongly with various environmental variables using multiple univariate logistic regression, latent factor mixed model, and Bayesian logistic regression. We found several environmentally dependent adaptive loci with high frequencies in low- or high-elevation populations, suggesting their involvement in local adaptation. Ecological factors, including relative humidity and sunshine hours, that are generally not altitude related could have been the most important selective drivers for population divergent evolution in C. formosana. The present study provides fundamental information in relation to adaptive evolution and can be useful for assisted migration program of C. formosana in the future conservation of this species.

Abstract Genetic variation evolves during postglacial range expansion of a species and is important for adapting to varied environmental conditions. It is crucial for the future survival of a species. We investigate the nuclear DNA sequence variation to provide evidence of postglacial range expansion of Musa basjoo var. formosana, a wild banana species, and test for adaptive evolution of amplified fragment length polymorphic (AFLP) loci underlying local adaptation in association with environmental variables. Postglacial range expansion was suggested by phylogeographical analyses based on sequence variation of the second intron of copper zinc superoxide dismutase 2 gene. Two glacial refugia were inferred by the average F ST parameter (mean F ST of a population against the remaining populations). Using variation partitioning by redundancy analysis, we found a significant amount of explained AFLP variation attributed to environmental and spatially-structured environmental effects. By combining genome scan methods and multiple univariate logistic regression, four AFLP loci were found to be strongly associated with environmental variables, including temperature, precipitation, soil moisture, wet days, and surface coverage activity representing vegetation greenness. These environmental variables may have played various roles as ecological drivers for adaptive evolution of M. basjoo var. formosana during range expansion after the last glacial maximum.

Populations can be genetically isolated by differences in their ecology or environment that hampered efficient migration, or they may be isolated solely by geographic distance. Moreover, mountain ranges across a species’ distribution area might have acted as barriers to gene flow. Genetic variation was quantified using amplified fragment length polymorphism (AFLP) and 13 selective amplification primer combinations used generated a total of 482 fragments. Here, we tested the barrier effects of mountains on gene flow and environmentally dependent local adaptation of Cunninghamia konishii occur in Taiwan. A pattern of genetic isolation by distance was not found and variation partitioning revealed that environment explained a relatively larger proportion of genetic variation than geography. The effect of mountains as barriers to genetic exchange, despite low population differentiation indicating a high rate of gene flow, was found within the distribution range of C. konishii. Twelve AFLP loci were identified as potential selective outliers using genome-scan methods (BAYESCAN and DFDIST) and strongly associated with environmental variables using regression approaches (LFMM, Samβada, and rstanarm) demonstrating adaptive divergence underlying local adaptation. Annual mean temperature, annual precipitation, and slope could be the most important environmental factors causally associated with adaptive genetic variation in C. konishii. The study revealed the existence of physical barriers to current gene flow and environmentally dependent adaptive divergence, and a significant proportion of the rate of gene flow may represent a reflection of demographic history.

The question of what determines divergence both between and within species has been the central topic in evolutionary biology. Neutral drift and environmentally dependent divergence are predicted to play roles in driving population and lineage divergence. However, neutral drift may preclude adaptation if the rate of gene flow between populations is high. Here, we sampled populations of three Taiwania (Taiwania cryptomerioides) lineages occurring in Taiwan, the mainland of China (Yunnan-Myanmar border), and northern Vietnam, and tested the relative strength of neutral drift and divergent selection in shaping divergence of those populations and lineages. We quantified genetic and epigenetic variation, respectively, using amplified fragment length polymorphism (AFLP) and methylation-sensitive amplification polymorphism (MSAP). Analysis of 1413 AFLP and 462 MSAP loci using frequency-based genome scan methods and generalized linear models (GLMs) found no potential selective outliers when only Taiwanese populations were examined, suggesting that neutral drift was the predominant evolutionary process driving differentiation between those populations. However, environmentally associated divergence was found when lineages were compared. Thirty-two potential selective outliers were identified based on genome scans and their associations with environmental variables were tested with GLMs, generalized linear mixed effect models (GLMMs), and model selection with a model averaging approach. Ten loci (six AFLP and four MSAP) were found to be strongly associated with environmental variables, particularly monthly temperature variation and normalized difference vegetation index (NDVI) using model selection and a model averaging approach. Because only a small portion of genetic and epigenetic loci were found to be potential selective outliers, neutral evolutionary process might also have played crucial roles in driving lineage divergence, particularly between geographically and genetically isolated island and mainland Asia lineages. Nevertheless, the vast amount of neutral drift causing genetic and epigenetic variations might have the potential for adaptation to future climate changes. These could be important for the survival of Taiwania in different geographic areas.

Double digest restriction site-associated DNA sequencing (ddRADseq) is a tool for delivering genome-wide single nucleotide polymorphism (SNP) markers for non-model organisms useful in resolving fine-scale population structure and detecting signatures of selection. This study performs population genetic analysis, based on ddRADseq data, of a coniferous species, Keteleeria davidiana var. formosana, disjunctly distributed in northern and southern Taiwan, for investigation of population adaptive divergence in response to environmental heterogeneity. A total of 13,914 SNPs were detected and used to assess genetic diversity, FST outlier detection, population genetic structure, and individual assignments of five populations (62 individuals) of K. davidiana var. formosana. Principal component analysis (PCA), individual assignments, and the neighbor-joining tree were successful in differentiating individuals between northern and southern populations of K. davidiana var. formosana, but apparent gene flow between the southern DW30 population and northern populations was also revealed. Fifteen of 23 highly differentiated SNPs identified were found to be strongly associated with environmental variables, suggesting isolation-by-environment (IBE). However, multiple matrix regression with randomization analysis revealed strong IBE as well as significant isolation-by-distance. Environmental impacts on divergence were found between populations of the North and South regions and also between the two southern neighboring populations. BLASTN annotation of the sequences flanking outlier SNPs gave significant hits for three of 23 markers that might have biological relevance to mitochondrial homeostasis involved in the survival of locally adapted lineages. Species delimitation between K. davidiana var. formosana and its ancestor, K. davidiana, was also examined (72 individuals). This study has produced highly informative population genomic data for the understanding of population attributes, such as diversity, connectivity, and adaptive divergence associated with large- and small-scale environmental heterogeneity in K. davidiana var. formosana.

The present study investigated the genetic diversity, population structure, F ST outliers, and extent and pattern of linkage disequilibrium in five populations of Keteleeria davidiana var. formosana, which is listed as a critically endangered species by the Council of Agriculture, Taiwan. Twelve amplified fragment length polymorphism primer pairs generated a total of 465 markers, of which 83.74% on average were polymorphic across populations, with a mean Nei's genetic diversity of 0.233 and a low level of genetic differentiation (approximately 6%) based on the total dataset. Linkage disequilibrium and HICKORY analyses suggested recent population bottlenecks and inbreeding in K. davidiana var. formosana. Both STRUCTURE and BAPS observed extensive admixture of individual genotypes among populations based on the total dataset in various clustering scenarios, which probably resulted from incomplete lineage sorting of ancestral variation rather than a high rate of recent gene flow. Our results based on outlier analysis revealed generally high levels of genetic differentiation and suggest that divergent selection arising from environmental variation has been driven by differences in temperature, precipitation, and humidity. Identification of ecologically associated outliers among environmentally disparate populations further support divergent selection and potential local adaptation.

Expressed sequence tags-derived simple sequence repeats (EST-SSRs) are related to functional genes and can be useful in crop breeding. Parental lines selected in sweet potato (Ipomoea batatas (L.) Lam.) breeding including controlled hybridization and polycross may be originated from a limited number of genetic sources. We surveyed genetic variation of three breeding populations (hybrid, landrace, and polycross) of 89 cultivated sweet potato varieties derived from controlled hybridization (41 hybrids) and polycross (38 varieties) breeding population and landrace (10 varieties) maintained by various geographic regions and main agronomic traits of local farmers, using 23 EST-SSRs. Six wild sweet potato (I. trifida) samples were also included in the study. High genetic diversity in the three types of breeding populations was detected. High proportions of shared alleles leading to low levels of differentiation between breeding populations were discovered. The neighbor-joining tree and individual assignments showed intermingling of varieties among the three breeding populations. Additionally, genetic clustering showed little power of grouping varieties of different breeding populations into separate clades. These results indicate limits of sweet potato breeding. Nevertheless, ten alleles of nine loci showed variance contributed significantly to the delimitation of varieties in different breeding populations. [ABSTRACT FROM AUTHOR]

Inferring the evolutionary history of a group of closely related species can be challenging. Genetic diversity, structure, and multilocus phylogeny were studied to shed light on the evolutionary processes that shapes diversity in three closely related Lecanorchis species (L. cerina, L. suginoana, and L. thalassica). In this study, we surveyed nucleotide and genetic variation, respectively, using nuclear internal transcribed spacer 1 (nrITS1) and amplified fragment length polymorphism (AFLP) in individuals of Lecanorchis species. Both nrITS1 parsimonious and AFLP neighbor-joining trees showed paraphyly indicated that lineage sorting for these species may be incomplete. ITS1 haplotype sharing was found within L. cerina, within L. suginoana, and between L. suginoana and L. thalassica. Higher AFLP diversity was observed in L. cerina than that in L. suginoana and in L. thalassica, indicating larger population size in L. cerina. Two contrasting patterns of AFLP genetic differentiation between geographic neighbors relative to geographic distant populations were observed. Two major AFLP genetic components were found in individuals of Lecanorchis species examined. Individuals of L. cerina and of L. thalassica were classified into either component with high probabilities. Most individuals of L. suginoana were classified into one of these components, but two individuals showed apparent admixture of these components. The evidence of paraphyly, haplotype sharing, higher L. cerina diversity, sharing and admixture of genetic components, and the contrasting patterns of genetic differentiation suggest that both incomplete lineage sorting and secondary gene flow could have occurred in these closely related species, though at different temporal scale. [ABSTRACT FROM AUTHOR]

An attempt is made to study genetic differentiation between populations of Querus glauca in Taiwan using nuclear microsatellite markers to infer the potential refugium in the last glaciation stage. It was found that the southern part of Taiwan showed the highest population genetic divergence compared to all the others which suggested the existence of a potential major refugium.

The genus Chamaecyparis (Cupressaceae) is comprised of six taxa native to East Asia and North America. Two independent lineages from eastern North America to East Asia and from western North America to East Asia have been identified from phylogenetic analyses based on variation in plastid DNA (ptDNA). However, the trees inferred from a nuclear ribosomal internal transcribed spacer (nrITS), the NEEDLY intron 2, and combined data from genetic and morphological characters of other studies are incongruent with the tree inferred from ptDNA. In this study, we sequenced the plastid matK gene as well as other sequences acquired from GenBank including the matK, rbcL, and nrITS sequences to reconstruct the phylogenetic trees of Chamaecyparis. The phylogenetic topology inferred from matK was congruent with the previous inference obtained using other ptDNA markers while the rbcL and ITS trees were congruent with previously inferred trees using the combined data and ITS sequences, respectively. Relatively less- (rbcL) and more-informative sites (ITS) may lead to different lineage sorting and incongruent phylogenetic topologies, which were evidenced by rejection of evolutionary homogeneity between ptDNA and ITS sequences in the partition homogeneity test. The phylogenetic tree reconstructed using matK and other ptDNA strongly supports the geographically disjunct distribution of Chamaecyparis in North America and East Asia. A dispersal-vicariance analysis, and geologic and fossil evidence indicated that at least two independent dispersal events occurred from North America to East Asia, which support the previous biogeographic inference by ptDNA. The use of biparental inherited markers for biogeographic inferences should be done with caution. [ABSTRACT FROM AUTHOR]

Cinnamomum kanehirae Hayata (Lauraceae), the most valuable subtropical and temperate broadleaf timber tree in Taiwan, is rapidly disappearing from the wild. Taking advantage of a scion garden established by the Taiwan Forestry Research Institute, we examined patterns of chloroplast DNA (cpDNA) variations in 19 populations including 94 individuals. By sequencing two cpDNA fragments using universal primers (the trnL- trnF and petG- trnP intergenic spacers), we found eight polymorphic sites, six haplotypes, and extremely low nucleotide diversity ( π = 0.00016) from 792 bp aligned sequences. The ancestral haplotype is widely distributed. Among the populations studied, three separated populations, at Yungfeng, Fuli, and Tahu have high nucleotide diversity. No phylogeographical structures of haplotypes were revealed because the tests of NST− GST for populations did not differ from zero in any situations; a ‘star-like’ genealogy is characteristic when all haplotypes rapidly coalesce and is a general outcome of population expansion. The neutrality test also suggested demographic expansion. The genetic divergence and diversity analyses suggested that two potential refugia existed during the last glaciation with a major one located in southeastern Taiwan and a minor one located in Tahu in north-central Taiwan in the Hsuehshan Range, west of the Central Mountain Range. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM PUBLISHER]

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, A VR-pita. The genotype variation and resistant/susceptible phenotype at the Pi-ta locus of wild rice (Oryza rufipogon), the ancestor of cultivated rice (O. saliva), 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 πnon/π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. [ABSTRACT FROM AUTHOR]

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 (GST = 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. [ABSTRACT FROM AUTHOR]

Aim The purpose of this paper was to study the patterns of genetic variation, demographic history, haplotype relationships and potential location of diversity centres of two closely related species, Machilus thunbergii and Machilus kusanoi. Location The phylogeography of M. thunbergii and M. kusanoi was examined by sampling 110 and 106 individuals from 25 and 16 sampling sites, respectively, across their distributional range in Taiwan. Machilus thunbergii is distributed on the Asian mainland, South Korea, southern Japan, the Ryukyus, Taiwan and the Philippines, whereas M. kusanoi is endemic to Taiwan. These two species are closely related, and both are widely distributed in Taiwan but occupy different altitudinal zones and habitats. Methods The range-wide variation of M. thunbergii and M. kusanoi in Taiwan was studied using chloroplast DNA (cpDNA) variations. A haplotype network was constructed with the computer programtcs. Nested clade analysis was conducted with the computer programceodis, and various parameters of genetic diversity were calculated and neutrality tested by the computer program DnaSP. Population differentiation was estimated using the programsarlequin andhapstep. The contribution of the populations to gene diversity and to allelic richness was calculated using the softwarecontrib. The level of divergence for each population from the remaining populations was calculated as the mean value of pairwise FST for each population against the rest of the populations. Results Extremely low levels of genetic differentiation were found for both species. This result suggested that these two species probably survived in multiple relict refugia with different population sizes throughout the island during low-temperature periods of the Pleistocene. In addition, nested clade analysis (NCA) of cpDNA haplotypes indicated that restricted gene flow with isolation-by-distance characterized the recolonization after the Pleistocene by Tashueshan and Shiouhluan populations of M. thunbergii in the north-central area west of the Central Mountain Range (CMR). In contrast, NCA analysis indicated that a major diversity centre on the southern tip of the island (Kending population) and contiguous range expansion characterized the recolonization by M. kusanoi of northern areas along the east side of the CMR. The major diversity centres found for the two species examined were further supported by the results of the mean FST for individual populations in comparison with other populations, and of the contribution of the divergence component to the total diversity. Main conclusions This research supports the multiple relict refugia hypothesis for both species investigated. Populations of M. thunbergii at Shiouhluan and Tashueshan in the north-central area west of the CMR represent a diversity centre currently expanding its size. A diversity centre at the southern-edge population of M. kusanoi, and a contiguous range expansion from Kending, were found. These results indicate that the M. thunbergii populations at Tashueshan and Shiouhluan and the M. kusanoi population at Kending, and even Soukar, are evolutionarily significant units for conservation programmes. [ABSTRACT FROM AUTHOR]

In this study, we examined spatial patterns of chloroplast DNA (cpDNA) variation in a total of 30 populations ofCastanopsis carlesiiHayata (Fagaceae), a subtropical and temperate tree species, including 201 individuals sampled throughout Taiwan. By sequencing two cpDNA fragments using universal primers (thetrnL intron and thetrnV-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 ofC. carlesii. Among the populations studied, three separated populations, at Lienhuachih, Fushan and Lichia, have high nucleotide diversity. Estimates ofNST– 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 forQuercus 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. [ABSTRACT FROM AUTHOR]

This paper described current phylogeographical patterns of chloroplastic DNA variation of Trochodendron aralioides, a temperate tree species, and inferred its possible refugium in Taiwan. This information was compared with the known phylogeographical pattern of subtropical tree species. A total of 24 populations were sampled including 20 from Taiwan, two each from the Ryukyus and Japan. A haplotype network was constructed by computer program TCS, various parameters of genetic diversity were calculated and neutrality was tested by computer program DnaSP. To examine the similarity of genetic structure among populations, a maximum parsimony tree was reconstructed by computer program PAUP*. The results of isozyme of T. aralioides from a previous publication were incorporated into this study to infer the phylogeographical history. Nine haplotypes according to six substitutions, two indels and one inversion of the two cpDNA intergenic spacer fragments ( petG- trnP and petA- psbJ) of T. aralioides were recognized. Genetic structure of the population of Japan is totally different from those of Taiwan and the Ryukyus. In Taiwan, the genetic structure was differentiated among populations revealed by Gst = 0.700 and Nst = 0.542, and the population genetics was clearly spatially structured. Two population groups were recognized. The first group was distributed islandwide and extended to the Ryukyus. The second group contained five of the seven known haplotypes, and was restricted to the area between latitude 24°46′ and 24°06′ N. In Taiwan, north-central area between latitude 24°46′ and 24°06′ N is potentially a refugium during the last glaciations. This finding is contradicted to subtropical species as Cyclobalanopsis glauca. [ABSTRACT FROM AUTHOR]

Abstract Phylogeographical and mismatch analysis of chloroplast DNA (cpDNA) variation were used to infer the temporal dynamics of distributional and demographic history of Taiwan fir (Cunninghamia konishii ). We examined 64 and 52 trees from 17 populations of C. konishii and 14 provenances of C. lanceolata , respectively, by sequencing three intergenic spacers and one intron using cpDNA universal primers. Of the aligned 1888 base pairs (bp) sequence, 30 varied among 28 haplotypes, which consisted of three transitions, 14 transversions and 13 indels. One ancestral haplotype was found in 86 individuals across the surveyed range of both species, C. konishii and C. lanceolata , which was distributed in all populations and provenances. The 28 haplotypes also included 15 C. konishii specific and 12 C. lanceolata -specific haplotypes. Ancestral haplotype was found fixed in five populations of C. konishii and five provenances of C. lanceolata . Other haplotypes occurred mainly as singletons. The levels of population differentiation studied are relatively low in both Cunninghamia species. The nucleotide diversity (θ ) of chloroplast DNA sequences within C. konishii was slightly higher than that of C. lanceolata . Excess in singletons as well as star-like phylogeny of haplotypes suggested no clearcut migration patterns of C. konishii after glacial maximum. One probable demographic history of C. konishii is the postglacial population growth of C. konishii after a glacial bottleneck event. This inference is supported by the combined results of fossil pollen record, low nucleotide diversity, significant Tajima's d -value, phylogeographical analysis and unimodal mismatch distribution. Similarities and discrepancies between our results and those of Lu et al. (2001) are discussed. [ABSTRACT FROM AUTHOR]

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 α-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. [ABSTRACT FROM AUTHOR]