BACKGROUND
the analysis of divergence patterns between species has been a major focus in systematics and molecular ecology  <cit> .
species descending from a common ancestor are expected to differentiate from each other and eventually attain reciprocal monophyly.
coalescence theory predicts that sister species are likely to be polyphyletic/paraphyletic in their genetic composition during the early stages after they split from their common ancestor .
as sufficient time goes by, they become reciprocally monophyletic with the loss of complementary haplotypes.
in contrast, lack of monophyly can be a result of several factors, such as incomplete lineage sorting, gene duplication, or hybridization  <cit> .
of them, hybridization leading to speciation has been known as one of the major forces for the diversification in plants  <cit> .
as many hybridizing species remain morphologically discrete  <cit> , hybrids usually share traits with parental species, often resulting in taxonomic difficulties.
interspecific gene flow usually leads to a mixed genetic composition in hybrids, but may have low impacts on parental species' genomes .
although molecular tools provide sufficient power in detecting genetic variation, empirical data on species delimitation are often difficult to interpret, since both gene introgression after speciation and shared ancestral polymorphisms can result in species paraphyly  <cit> .
a complex of incipient species with different degrees of morphological or ecological differentiation provides an ideal model for studying species divergence.
taiwan, an island of the island-arc system along the western edge of the pacific ocean, first emerged from the waters via collision between eurasian and pacific plates about  <dig> million years before present , but attained its modern shape only 5- <dig> million years ago   <cit> .
its rugged topography is characterized by hundreds of steep mountains, which in turn provide diverse habitats along distinct vegetation zones, including tropical, coastal evergreen forests, subalpine shrubs and alpine tundra along the central mountain range  <cit> .
taiwan was connected to mainland se asia by a land bridge during lower sea levels of the last glacial maximum  <cit> ; some species may have reached taiwan by this route and subsequently become endemic  <cit> .
phylogeographical patterns and genetic characteristics of populations/species are shaped by the interaction between historical vicariance and recurrent genetic exchanges.
these evolutionary events would leave evolutionary footprints in the genetic polymorphisms within and between populations across the distributional range.
of the geological events, regular historical glacial cycles in the eurasia continent had prevalent influences on survival and recolonization of populations/species  <cit> .
postglacial expansion of alpine species may follow the phalanx model, which describes the effects of slower expansions from refugia due to habitat constraints  <cit> .
that is, during interglacial periods, alpine species moved upward to peaks and became fragmented.
in contrast, down-slope movement of alpine species during glacial periods could connect populations that were previously isolated and allow for genetic exchange among populations  <cit> .
on taiwan island, only the peaks of high mountains over  <dig> m in elevation were covered by ice sheets during the glacial maxima.
rhododendron, one of the largest and most widespread woody plant genera, is distributed from the northern temperate zone, throughout tropical southeast asia, to northeastern australia and consists of over  <dig> species  <cit> .
hymenanthes, one of eight subgenera comprising rhododendron, consists of  <dig> species, most of them distributed in the himalaya-southwest china region  <cit> .
in taiwan,  <dig> taxa occur within rhododendron  <cit> .
the r. pseudochrysanthum complex consists of r. pseudochrysanthum hayata sensu stricto , r. morii hayata, and r. rubropunctatum hayata, all in subgenus hymenanthes subsection maculifera, and all are phylogenetically related based on molecular analyses  <cit> .
the remaining taxa of subsection maculifera occur in west and central china, disjunctly from those in taiwan  <cit> .
r. pseudochrysanthum s.s. is an alpine species , while r. morii has a wider range from 1650- <dig> m a.s.l., with some overlap in distribution.
in contrast, r. rubropunctatum is restricted to elevations of 600- <dig> m a.s.l.
in northern taiwan.
the three species of r. pseudochrysanthum complex, all distributed along the central mountain range  and other mountain ranges west of the cmr, share morphological polymorphisms in leaf shape, petiole and flower shape, thus making taxonomy difficult, despite the differentiation of ecological distributions  <cit> .
possibly the closest relative of this species complex is r. hyperythrum hayata, which is also endemic to taiwan, occurring 1400- <dig> m a.s.l.
<cit> .
this species was previously placed in subsection pontica  <cit> , but phylogenetic work indicates that it is closely related to r. pseudochrysanthum  <cit> .
these two might be sister species, but at present a phylogeny of subsection maculifera is lacking, so this hypothesis cannot yet be tested.
rhododendron hyperythrum is also known to hybridise with members of the r. pseudochrysanthum complex  <cit> .
closely related to r. hyperythrum, the three species of the r. pseudochrysanthum complex with low morphological differentiations may be recognized as a single species   <cit> .
the power of molecular markers in elucidating plant phylogeography has been detected in many studies  <cit> .
comparative phylogenies between biparental  and uniparental  markers are useful for resolving phylogeographic patterns, origin of hybridization, and assessing the migratory routes of species  <cit> .
because one or both organelle genomes are often maternally inherited in plants  <cit> , they are particularly suitable for investigating processes associated with seed dispersal, such as range expansions  <cit> , contribution of seed movement to total gene flow  <cit> , refugial isolation  <cit> , and origin of maternal or paternal lineage  <cit> .
previous studies have investigated the phylogeography of the r. pseudochrysanthum in taiwan, revealing close relationships among these taxa  <cit>  and a lack of geographical subdivision.
however, little attention has been paid to the demographic dynamics of r. pseudochrysanthum across the mountain ranges.
in contrast to previous studies, sampling in this study covered vaster areas, with the aid of the bureau of the yushan national park.
the present study analyzed genetic variation of the r. pseudochrysanthum s.l.
at atpb-rbcl intergenic spacer of cpdna and the internal transcribed spacer regions  of nrdna.
a phylogeographical approach was used to infer past evolutionary history and processes and to examine species delimitation and population structure of the species complex.
the objectives of this study were to investigate the following:
<dig>  do the lineages of cpdna correspond with those of nrdna?
<dig>  given geographical barriers, does significant genetic differentiation exist among populations of different mountain ranges?
<dig>  have the populations of different mountain ranges of r. pseudochrysanthum s.l.
experienced different demographic dynamics?
RESULTS
genetic variability of cpdna and nrdna and genetic differentiation in the rhododendron pseudochrysanthum s.l
the cpdna atpb-rbcl intergenic spacer of r. formosanum, r. hyperythrum, and r. pseudochrysanthum s.l.
ranged from  <dig> to  <dig> bp.
the sequence of alignment of r. pseudochrysanthum s.l.
had  <dig> bp, including  <dig> polymorphic sites, of which  <dig> were parsimony-informative.
most of these polymorphic sites were single-base and scattered along a sequence.
the nrdna its sequences of r. formosanum, r. hyperythrum, and r. pseudochrysanthum s.l.
were aligned with a consensus length of  <dig> bp.
no large indels or rearrangements were detected.
the sequence alignment of r. pseudochrysanthum s.l.
had  <dig> bp, including  <dig> polymorphic sites, of which  <dig> were parsimony informative.
with regard to cpdna,  <dig> haplotypes  were identified among  <dig> samples of the r. pseudochrysanthum s.l.
.
nucleotide diversity  detected in r. pseudochrysanthum s.l.
was  <dig>  with high level of haplotype diversity  .
high levels of nucleotide diversity  and haplotype diversity  were detected in the populations of the yushan mountain range, while the single population sampled from the alishan mountain range was fixed at a single haplotype .
ph <dig> was the most common haplotype shared among the alishan, central, and sheishan mountain ranges.
in addition, no haplotypes in the yushan mountain range  were shared by other mountain ranges.
samples size, nucleotide diversity, haplotype number, haplotype diversity, and tajima's d are indicated.
significance values are given in asterisks .
for nrdna,  <dig> haplotypes  were identified among  <dig> samples of the r. pseudochrysanthum s.l.
.
nucleotide diversity detected in r. pseudochrysanthum s.l.
was  <dig>  with a high level of haplotype diversity .
a high level of nucleotide diversity  with a high level of haplotype diversity  were detected in the populations of the yushan mountain range, while the lowest level of nucleotide diversity  was detected in a single population of the alishan mountains .
to investigate the genetic structure of r. pseudochrysanthum s.l., a samova  was applied to define groups and to identify any locations of genetic uniqueness among the  <dig> populations.
for cpdna, an assumption with two groups  displayed the greatest value of fct and a maximal variance .
the groups at k =  <dig> exactly matched those found in the phylogenetic and network analyses, with those from the yushan mountain range forming one group  and those from the other three mountain ranges forming the other.
genetic structure was significant at three levels based on the spatial apportionment of genetic variation, i.e., among geographical groups , among populations within groups , and within populations .
alishan mountain range: aa; central mountain range: cb, ck, cy, cs, cn, ch; sheishan mountain range: ss, sl, sb, sc; yushan mountain range: yt, yy, yp.
the samova of nrdna sequences also revealed the greatest value of fct and a statistically maximal variance between two groups  .
likewise, a large proportion of the genetic variability was partitioned among groups, although nonsignificantly .
in contrast, when the three geographical groups of sl vs. sc vs. all others were compared, structuring with significant variabilities distributed among groups , among populations within groups , and within populations  was identified.
significant genetic differentiation of cpdna was detected between r. pseudochrysanthum s.l.- r. formosanum  and r. pseudochrysanthum s.l.
- r. hyperythrum .
within r. pseudochrysanthum s.l., significant genetic differentiation existed between the populations of the yushan mountain range and the populations of the other ranges  , whereas there was only little differentiation among the populations of the alishan, central and sheishan mountain ranges .
in contrast, lower levels of genetic differentiation of nrdna were detected between r. pseudochrysanthum s.l.
- r. formosanum  and r. pseudochrysanthum s.l.
- r. hyperythrum , and the populations of different mountain ranges of r. pseudochrysanthum s.l.
.
gene genealogies and phylogeographical analyses of the rhododendron pseudochrysanthum s. l
the ml tree of the cpdna atpb-rbcl intergenic spacer haplotypes , rooted at r. formosanum, identified two major clusters i and ii.
a minimum-spanning network was also reconstructed based on the mutational steps between haplotypes of cpdna .
rhododendron pseudochrysanthum s.l.
was identified with eight mutational steps from r. formosanum.
two major clusters i and ii, corresponding to the ml tree, were identified with two mutational steps apart.
a ml tree and a minimum-spanning network were reconstructed of the nrdna haplotypes, respectively .
in total, clusters a-c were identified, none of which corresponded to any member of species.
we also utilized the *beast program to analyze the combined data of cp- and nrdna sequences.
a population tree of r. pseudochrysanthum s.l., rooted at r. formosanum  and r. hyperythrum , was reconstructed .
in contrast to the nuclear dna tree, none of the outgroup populations were nested within r. pseudochrysanthum s.l.
mismatch distribution analyses of cpdna were perfomed to infer the long-term demographic history of populations.
for r. pseudochrysanthum s.l., populations of the central and sheishan mountain ranges had a unimodal distribution .
because the data used to produce mismatch distributions are not independent  <cit> , tajima's d was then used to detect departure from population equilibrium.
departures from the null model can be caused by many factors, such as changes in population size , which can lead to an excess of low frequency variants and negative d values.
in contrast, processes such as population subdivision, or recent population bottlenecks can cause an excess of intermediate frequency variants leading to positive d values  <cit> .
most values of tajima's d in r. pseudochrysanthum s.l.
were negative, except in yt  and yy  of the yushan mountain range .
using the program im, we were able to refine six model parameters estimated in mcmc simulations and test for the isolation with migration model in r. pseudochrysanthum s.l.
all of our parameter estimates had unimodal posterior distributions.
the posterior distribution of θother mtns  was larger than that of θyushan , whereas the posterior distribution for θa  peaked at  <dig>  .
based on these θ values, the populations of the alishan, central and sheishan mountain ranges were likely to grow substantially following divergence, whereas the populations of the yushan mountain range remained stable or slightly grew relative to the ancestral population.
when the mles for the migration parameters were transformed into population migration rates, 2nem = m* was estimated among the populations of different mountain ranges.
results from the im revealed that the number of migrants between the yushan mountain range and the other mountain ranges was effectively low.
the posterior distribution of m  between the yushan mountain range and the other mountain ranges peaked at  <dig> , i.e., both myushan to other mtns =  <dig>   and mother mtns to yushan =  <dig>  , which were much lower than the m values among alishan, central, and sheishan mountain ranges .
all of the non-zero m values, nevertheless, implied the occurrence of historical gene flow between populations .
to elucidate the contribution of hybridization/introgression to the lack of species monophyly of r. formosanum, r. hyperythrum and r. pseudochrysanthum s.l., we estimated the m values among these species.
the posterior distribution of m among these species peaked at non-zero values , suggesting possible interspecific hybridization/introgression.
molecular dating of rhododendron pseudochrysanthum s.l
bayesian estimates of the divergence time of clades of the sequences of cpdna and nrdna were obtained using beast v.  <dig> .
for cpdna, the sequences including r. formosanum, r. hyperythrum, and r. pseudochrysanthum s.l.
could be traced back to a common ancestor ca.
<dig> - <dig>  mya.
all sequences of r. hyperythrum and r. pseudochrysanthum s.l.
coalesced at ca.
<dig> - <dig>  mya.
all cpdna sequences of r. pseudochrysanthum s.l.
coalesced at ca.
<dig> - <dig>  mya, while those of r. hyperythrum coalesced at ca.
<dig> - <dig>  mya.
likewise, the nrdna sequences including r. formosanum, r. hyperythrum, and r. pseudochrysanthum s.l.
could be traced back to a common ancestor about  <dig> - <dig>  mya.
all nrdna sequences of r. pseudochrysanthum s.l.
and r. hyperythrum coalesced to a common ancestor ca.
<dig> - <dig>  mya.
all nrdna sequences of r. pseudochrysanthum s.l.
coalesced to a common ancestor ca.
<dig> - <dig>  mya, while those of r. hyperythrum coalesced at ca.
<dig> - <dig>  mya.
to elucidate the possible evolutionary forces for the lack of species monophyly of r. formosanum, r. hyperythrum and r. pseudochrysanthum s.l., we calculated the expected time to coalescence at cpdna  and at nrdna, where θ is a population mutation parameter obtained from im, and u is the geometric mean of mutation rates, and compared them with the divergence times obtained from the beast analyses.
according to the im estimates of θother mtns =  <dig>  and θyushan =  <dig> , the expected time to coalescence of r. pseudochrysanthum s.l.
ranged from  <dig>  to  <dig>  mya based on cpdna, and from  <dig>  to  <dig>  mya based on nrdna.
discussion
systematic inconsistency between nuclear and chloroplast dnas in rhododendron pseudochrysanthum s.l
inconsistency in the cpdna and nrdna phylogenies  was detected in r. formosanum.
the ml tree for the cpdna agreed with the results of chamberlain et al.
<cit>  and goetsch et al.
<cit>  that r. formosanum is a monophyletic group relative to r. pseudochrysanthum s.l.
in contrast, the ml tree of nrdna revealed that one accession of r. formosanum fell within cluster b with the ingroup species.
lack of species monophyly may result from incomplete lineage sorting, which occurs when populations or species recently diverged, thus retaining ancestral genetic polymorphisms, or from genetic exchanges via hybridization and introgression that introduce foreign alleles into a species or population, resulting in a genetic admixture  <cit> .
for given species, the time for a locus to acquire reciprocal monophyly is about  <dig>  × θ/4u at cpdna or  <dig>  × θ/2u at nrdna.
willyard et al.
<cit>  suggest that ancestral polymorphisms would likely be shared between species if the coalescence time of the alleles exceeds the divergence time at one locus.
in this study, the expected time to coalescence of r. pseudochrysanthum s.l.
overlapped the estimates of divergence time  based on cpdna, in which monophyly was attained.
in contrast, at the nrdna, the expected time for coalescence  predated the divergence time .
accordingly, nuclear loci would less likely be monophyletic in r. pseudochrysanthum s.l.
incomplete lineage sorting alone, apparently, cannot completely explain the paraphyly of r. pseudochrysanthum s.l.
interspecific hybridization that disturbs the reciprocal monophyly needs to be aroused based on the non-zero m parameters between r. formosanum-r. pseudochrysanthum s.l.
and r. hyperythrum-r. pseudochrysanthum s.l.
.
furthermore, milne et al.
<cit>  revealed that natural hybridization/introgression is not unusual in rhododendron.
thus, interspecific hybridization/introgression would cause the lack of monophyly and prolong the duration of paraphyly, as frequently occurs in most plants, e.g., ixeris  <cit> , begonia  <cit> , and ilex  <cit> .
genetic diversity and population structure of r. pseudochrysanthum s.l
in this study, we investigated genetic diversity at the cpdna and nrdna loci of r. pseudochrysanthum s.l.
high levels of nucleotide diversity were detected at both cpdna and nrdna loci, agreeing with a pattern found in many east asian plants, that is considered to have been shaped by a common glaciation history .
during the quaternary glaciations, taiwan, as a refuge, provided shelters for species from a range of biomes and thus harbored much plant biodiversity.
examples include cunninghamia konishii  <cit> , pinus luchuensis ssp.
taiwanensis  <cit> , cycas taitungensis  <cit> , and cycas revoluta  <cit> .
a potential refugium for r. pseudochrysanthum s.l.
in taiwan was inferred to harbor the genetic variations.
furthermore, highly diverse environments may also have created a wide spectrum of habitats to accommodate mutations  <cit> .
the populations of the yushan mountain range had the highest levels of genetic diversity  in the cpdna.
a similar pattern was found in taiwan fir , a dominant species of subalpine forests at 2500- <dig> m a.s.l.
<cit> .
this common result is likely to be attributable to a stable population size in the yushan mountain range.
the habitat ranges of the yushan mountain range, with the highest peak at  <dig> m, are wider than those of the other mountain ranges , providing diverse habitats for r. pseudochrysanthum s.l.
given a stable, effective population size under diverse habitats, the populations of the yushan mountain range fostered a great number of genetic polymorphisms.
in contrast, genetic diversity can be lost by chance in fragmented populations, as occurred in the population of the alishan mountain range.
ph <dig>  the interior dominant haplotype, is shared among the alishan, central, and sheishan mountain ranges.
a single population of the alishan mountain range was fixed for the dominant haplotype ph <dig>  indicating a possible genetic bottleneck due to a founder effect and/or genetic drift during a period when population sizes were small  <cit> .
genetic structure of r. pseudochrysanthum s.l
r. pseudochrysanthum s.l.
inhabits mountaintops isolated by surrounding low-elevation habitats, forming small and fragmented patches.
genetic differentiation would be detected between these isolated populations of different mountain ranges.
nevertheless, genetic analysis of the cpdna variation revealed very weak genetic differentiation among the populations within a mountain range.
the short time of isolation since the last glacial retreat is seemingly not long enough for differentiating these populations within a mountain range and thus resulted in sharing ancestral genetic polymorphisms.
specifically, spatial distribution of cpdna revealed a noteworthy pattern of high differentiation between populations of the yushan mountain range and the other mountain ranges  .
the divergence time for r. pseudochrysanthum s.l.
is about  <dig> - <dig>  mya, and given that within r. pseudochrysanthum s.l., the populations from different mountain ranges are more differentiated than the species, then it follows that this divergence time likely reflects an estimated divergence time between mountain populations.
molecular dating revealed that populations among different mountain ranges split some  <dig> - <dig>  mya, a time range approximating or predating the formation of the central mountain range of taiwan  <cit> .
independent colonizations via individuals with two different subclusters, likely from adjacent mainland over pleistocene, may result in high genetic differentiation within r. pseudochrysanthum s.l.
unexpected genetic structuring was detected with populations of the yushan mountain range highly differentiated from other mountain ranges, despite this range lying between the alishan and the central mountain ranges.
the single population of the alishan mountain range was fixed at the dominant haplotype ph <dig>  which was shared by the populations of the central and sheishan mountain ranges.
some other haplotypes may have been randomly lost from the alishan population due to possible bottleneck/founder's events.
despite the fact that the alishan mountain range is geographically distant to central and sheishan mountain ranges, sharing the dominant ph <dig> made these populations clustered together and genetically undifferentiated.
furthermore, subsequent glacial retreat and mountain uplifting triggered the geographical isolations and may have thereby enhanced the divergence of populations among different mountain ranges.
today, over an interglacial period, rhododendron populations of different mountain ranges represent some isolated patches.
given such geographical barriers, seeds are less likely to disperse across mountain ranges.
low capabilities of seed dispersal in current environments, plus possible independent colonization events followed by limited historical gene flow , may have together contributed to the high levels of genetic differentiation between populations of the yushan mountain range and the other mountain ranges .
in contrast, low levels of genetic differentiation at cpdna detected among the populations of alishan, central and sheishan mountain ranges  implied that there may exist historical gene flow during the glacial maxima, as the alpine populations migrated downwards and became connected to each other .
compared to the high genetic differentiation among populations between yushan mountain range and the other mountain ranges based on cpdna, nrdna revealed much less differentiation among populations of different mountain ranges  .
here, biparentally inherited nrdna was used to estimate the extent of pollen flow.
in contrast to the limited seed dispersal, pollen of rhododendron species seemingly can move long distances  <cit> .
long-distance pollen dispersal of r. pseudochrysanthum s.l.
across the discontinuity of mountain ranges reduced the genetic differentiation.
ancestral population size, gene flow, and population dynamics of r. pseudochrysanthum s.l
given their different inheritances, nrdna and cpdna markers provide extensive aspects to the understanding of the evolutionary forces shaping fates of r. pseudochrysanthum s.l.
combining cpdna and nrdna data, im analysis was used to evaluate the ancestral population size and gene flow.
populations of the alishan, central, and sheishan mountain ranges were characterized by a recent population expansion as θother mtns was approximately four times larger than ancestral θa  <cit> .
star-like genealogies of cpdna and negative tajima's d of cpdna and nrdna  also supported a scenario of recent demographic expansion in populations of the central and sheishan mountain ranges, a common pattern also occurring in castanopsis carlesii  <cit> , cunninghamia konishii  <cit> , and trema dielsiana  <cit> .
in contrast, the value of θyushan approximated ancestral θa, indicating that the effective population size of the yushan mountain range may have remained stable relative to the ancestral population.
given low seed dispersal capabilities of rhododendron in today's environments, non-zero migration rates  detected among populations  likely represented historical gene flow across mountain ranges.
here, compared to the gene flow among populations of alishan, central, and sheishan mountain ranges, with m =  <dig> - <dig> , much lower levels of migration rates were detected between populations of the yushan mountain range and other mountain ranges .
even so, higher levels of migration rates detected from populations of the yushan mountain range to those of the alishan  and sheishan  mountain range further indicated asymmetric gene flow among them.
in contrast, historical dispersal may have occurred freely among other populations.
that is, during glacial maxima, there may exist corridors among the mountain ranges.
these corridors allowed frequent gene flow among populations of alishan, central, and sheishan mountain ranges, and occasional gene flow between populations of the yushan mountain range and the other mountain ranges.
altogether, populations of these mountain ranges may have been mediated with different levels of historical gene flow during the glacial maxima.
CONCLUSIONS
in this study, inconsistent gene genealogies of cpdna and nrdna were detected.
in the cpdna tree, agreeing with chamberlain et al.
<cit> , r. pseudochrysanthum s.l.
is phylogenetically more closely related to r. hyperythrum than to r. formosanum.
in contrast, nrdna failed to identify reciprocal monophyly of ingroup vs. outgroup, indicating possible incomplete lineage sorting and interspecific hybridization/introgression.
for r. pseudochrysanthum s.l., the spatial distribution of cpdna revealed a noteworthy pattern of high differentiation between the populations of the yushan mountain range and the populations of the other mountain ranges.
molecular dating revealed that the split of the populations at different mountain ranges approximated or predated the formation of the central mountain range of taiwan island, likely suggesting independent colonizations via individuals from previously differentiated populations on the mainland.
at the population level, the populations of central, and sheishan mountain ranges may have undergone postglacial demographic expansion, while populations of the yushan mountain range are likely to have remained stable ever since the colonization.
in contrast, the single population of the alishan mountain range with a fixed cpdna haplotype may have experienced bottleneck/founder's events.

 the analysis of divergence patterns between species has been a major focus in systematics and molecular ecology. species descending from a common ancestor are expected to differentiate from each other and eventually attain reciprocal monophyly. coalescence theory predicts that sister species are likely to be polyphyletic/paraphyletic in their genetic composition during the early stages after they split from their common ancestor. as sufficient time goes by, they become reciprocally monophyletic with the loss of complementary haplotypes. in contrast, lack of monophyly can be a result of several factors, such as incomplete lineage sorting, gene duplication, or hybridization. of them, hybridization leading to speciation has been known as one of the major forces for the diversification in plants. as many hybridizing species remain morphologically discrete, hybrids usually share traits with parental species, often resulting in taxonomic difficulties. interspecific gene flow usually leads to a mixed genetic composition in hybrids, but may have low impacts on parental species' genomes. although molecular tools provide sufficient power in detecting genetic variation, empirical data on species delimitation are often difficult to interpret, since both gene introgression after speciation and shared ancestral polymorphisms can result in species paraphyly. a complex of incipient species with different degrees of morphological or ecological differentiation provides an ideal model for studying species divergence. taiwan, an island of the island-arc system along the western edge of the pacific ocean, first emerged from the waters via collision between eurasian and pacific plates about million years before present, but attained its modern shape only 5- million years ago. its rugged topography is characterized by hundreds of steep mountains, which in turn provide diverse habitats along distinct vegetation zones, including tropical, coastal evergreen forests, subalpine shrubs and alpine tundra along the central mountain range. taiwan was connected to mainland se asia by a land bridge during lower sea levels of the last glacial maximum; some species may have reached taiwan by this route and subsequently become endemic. phylogeographical patterns and genetic characteristics of populations/species are shaped by the interaction between historical vicariance and recurrent genetic exchanges. these evolutionary events would leave evolutionary footprints in the genetic polymorphisms within and between populations across the distributional range. of the geological events, regular historical glacial cycles in the eurasia continent had prevalent influences on survival and recolonization of populations/species. postglacial expansion of alpine species may follow the phalanx model, which describes the effects of slower expansions from refugia due to habitat constraints. that is, during interglacial periods, alpine species moved upward to peaks and became fragmented. in contrast, down-slope movement of alpine species during glacial periods could connect populations that were previously isolated and allow for genetic exchange among populations. on taiwan island, only the peaks of high mountains over m in elevation were covered by ice sheets during the glacial maxima. rhododendron, one of the largest and most widespread woody plant genera, is distributed from the northern temperate zone, throughout tropical southeast asia, to northeastern australia and consists of over species. hymenanthes, one of eight subgenera comprising rhododendron, consists of species, most of them distributed in the himalaya-southwest china region. in taiwan, taxa occur within rhododendron. the r. pseudochrysanthum complex consists of r. pseudochrysanthum hayata sensu stricto, r. morii hayata, and r. rubropunctatum hayata, all in subgenus hymenanthes subsection maculifera, and all are phylogenetically related based on molecular analyses. the remaining taxa of subsection maculifera occur in west and central china, disjunctly from those in taiwan. r. pseudochrysanthum s.s. is an alpine species, while r. morii has a wider range from 1650- m a.s.l., with some overlap in distribution. in contrast, r. rubropunctatum is restricted to elevations of 600- m a.s.l. in northern taiwan. the three species of r. pseudochrysanthum complex, all distributed along the central mountain range and other mountain ranges west of the cmr, share morphological polymorphisms in leaf shape, petiole and flower shape, thus making taxonomy difficult, despite the differentiation of ecological distributions. possibly the closest relative of this species complex is r. hyperythrum hayata, which is also endemic to taiwan, occurring 1400- m a.s.l.. this species was previously placed in subsection pontica, but phylogenetic work indicates that it is closely related to r. pseudochrysanthum. these two might be sister species, but at present a phylogeny of subsection maculifera is lacking, so this hypothesis cannot yet be tested. rhododendron hyperythrum is also known to hybridise with members of the r. pseudochrysanthum complex. closely related to r. hyperythrum, the three species of the r. pseudochrysanthum complex with low morphological differentiations may be recognized as a single species. the power of molecular markers in elucidating plant phylogeography has been detected in many studies. comparative phylogenies between biparental and uniparental markers are useful for resolving phylogeographic patterns, origin of hybridization, and assessing the migratory routes of species. because one or both organelle genomes are often maternally inherited in plants, they are particularly suitable for investigating processes associated with seed dispersal, such as range expansions, contribution of seed movement to total gene flow, refugial isolation, and origin of maternal or paternal lineage. previous studies have investigated the phylogeography of the r. pseudochrysanthum in taiwan, revealing close relationships among these taxa and a lack of geographical subdivision. however, little attention has been paid to the demographic dynamics of r. pseudochrysanthum across the mountain ranges. in contrast to previous studies, sampling in this study covered vaster areas, with the aid of the bureau of the yushan national park. the present study analyzed genetic variation of the r. pseudochrysanthum s.l. at atpb-rbcl intergenic spacer of cpdna and the internal transcribed spacer regions of nrdna. a phylogeographical approach was used to infer past evolutionary history and processes and to examine species delimitation and population structure of the species complex. the objectives of this study were to investigate the following: do the lineages of cpdna correspond with those of nrdna? given geographical barriers, does significant genetic differentiation exist among populations of different mountain ranges? have the populations of different mountain ranges of r. pseudochrysanthum s.l. experienced different demographic dynamics?  genetic variability of cpdna and nrdna and genetic differentiation in the rhododendron pseudochrysanthum s.l the cpdna atpb-rbcl intergenic spacer of r. formosanum, r. hyperythrum, and r. pseudochrysanthum s.l. ranged from to bp. the sequence of alignment of r. pseudochrysanthum s.l. had bp, including polymorphic sites, of which were parsimony-informative. most of these polymorphic sites were single-base and scattered along a sequence. the nrdna its sequences of r. formosanum, r. hyperythrum, and r. pseudochrysanthum s.l. were aligned with a consensus length of bp. no large indels or rearrangements were detected. the sequence alignment of r. pseudochrysanthum s.l. had bp, including polymorphic sites, of which were parsimony informative. with regard to cpdna, haplotypes were identified among samples of the r. pseudochrysanthum s.l.. nucleotide diversity detected in r. pseudochrysanthum s.l. was with high level of haplotype diversity. high levels of nucleotide diversity and haplotype diversity were detected in the populations of the yushan mountain range, while the single population sampled from the alishan mountain range was fixed at a single haplotype. ph was the most common haplotype shared among the alishan, central, and sheishan mountain ranges. in addition, no haplotypes in the yushan mountain range were shared by other mountain ranges. samples size, nucleotide diversity, haplotype number, haplotype diversity, and tajima's d are indicated. significance values are given in asterisks. for nrdna, haplotypes were identified among samples of the r. pseudochrysanthum s.l.. nucleotide diversity detected in r. pseudochrysanthum s.l. was with a high level of haplotype diversity. a high level of nucleotide diversity with a high level of haplotype diversity were detected in the populations of the yushan mountain range, while the lowest level of nucleotide diversity was detected in a single population of the alishan mountains. to investigate the genetic structure of r. pseudochrysanthum s.l., a samova was applied to define groups and to identify any locations of genetic uniqueness among the populations. for cpdna, an assumption with two groups displayed the greatest value of fct and a maximal variance. the groups at k = exactly matched those found in the phylogenetic and network analyses, with those from the yushan mountain range forming one group and those from the other three mountain ranges forming the other. genetic structure was significant at three levels based on the spatial apportionment of genetic variation, i.e., among geographical groups, among populations within groups, and within populations. alishan mountain range: aa; central mountain range: cb, ck, cy, cs, cn, ch; sheishan mountain range: ss, sl, sb, sc; yushan mountain range: yt, yy, yp. the samova of nrdna sequences also revealed the greatest value of fct and a statistically maximal variance between two groups. likewise, a large proportion of the genetic variability was partitioned among groups, although nonsignificantly. in contrast, when the three geographical groups of sl vs. sc vs. all others were compared, structuring with significant variabilities distributed among groups, among populations within groups, and within populations was identified. significant genetic differentiation of cpdna was detected between r. pseudochrysanthum s.l.- r. formosanum and r. pseudochrysanthum s.l. - r. hyperythrum. within r. pseudochrysanthum s.l., significant genetic differentiation existed between the populations of the yushan mountain range and the populations of the other ranges, whereas there was only little differentiation among the populations of the alishan, central and sheishan mountain ranges. in contrast, lower levels of genetic differentiation of nrdna were detected between r. pseudochrysanthum s.l. - r. formosanum and r. pseudochrysanthum s.l. - r. hyperythrum, and the populations of different mountain ranges of r. pseudochrysanthum s.l.. gene genealogies and phylogeographical analyses of the rhododendron pseudochrysanthum s. l the ml tree of the cpdna atpb-rbcl intergenic spacer haplotypes, rooted at r. formosanum, identified two major clusters i and ii. a minimum-spanning network was also reconstructed based on the mutational steps between haplotypes of cpdna. rhododendron pseudochrysanthum s.l. was identified with eight mutational steps from r. formosanum. two major clusters i and ii, corresponding to the ml tree, were identified with two mutational steps apart. a ml tree and a minimum-spanning network were reconstructed of the nrdna haplotypes, respectively. in total, clusters a-c were identified, none of which corresponded to any member of species. we also utilized the *beast program to analyze the combined data of cp- and nrdna sequences. a population tree of r. pseudochrysanthum s.l., rooted at r. formosanum and r. hyperythrum, was reconstructed. in contrast to the nuclear dna tree, none of the outgroup populations were nested within r. pseudochrysanthum s.l. mismatch distribution analyses of cpdna were perfomed to infer the long-term demographic history of populations. for r. pseudochrysanthum s.l., populations of the central and sheishan mountain ranges had a unimodal distribution. because the data used to produce mismatch distributions are not independent, tajima's d was then used to detect departure from population equilibrium. departures from the null model can be caused by many factors, such as changes in population size, which can lead to an excess of low frequency variants and negative d values. in contrast, processes such as population subdivision, or recent population bottlenecks can cause an excess of intermediate frequency variants leading to positive d values. most values of tajima's d in r. pseudochrysanthum s.l. were negative, except in yt and yy of the yushan mountain range. using the program im, we were able to refine six model parameters estimated in mcmc simulations and test for the isolation with migration model in r. pseudochrysanthum s.l. all of our parameter estimates had unimodal posterior distributions. the posterior distribution of θother mtns was larger than that of θyushan, whereas the posterior distribution for θa peaked at. based on these θ values, the populations of the alishan, central and sheishan mountain ranges were likely to grow substantially following divergence, whereas the populations of the yushan mountain range remained stable or slightly grew relative to the ancestral population. when the mles for the migration parameters were transformed into population migration rates, 2nem = m* was estimated among the populations of different mountain ranges. results from the im revealed that the number of migrants between the yushan mountain range and the other mountain ranges was effectively low. the posterior distribution of m between the yushan mountain range and the other mountain ranges peaked at, i.e., both myushan to other mtns = and mother mtns to yushan =, which were much lower than the m values among alishan, central, and sheishan mountain ranges. all of the non-zero m values, nevertheless, implied the occurrence of historical gene flow between populations. to elucidate the contribution of hybridization/introgression to the lack of species monophyly of r. formosanum, r. hyperythrum and r. pseudochrysanthum s.l., we estimated the m values among these species. the posterior distribution of m among these species peaked at non-zero values, suggesting possible interspecific hybridization/introgression. molecular dating of rhododendron pseudochrysanthum s.l bayesian estimates of the divergence time of clades of the sequences of cpdna and nrdna were obtained using beast v.. for cpdna, the sequences including r. formosanum, r. hyperythrum, and r. pseudochrysanthum s.l. could be traced back to a common ancestor ca. - mya. all sequences of r. hyperythrum and r. pseudochrysanthum s.l. coalesced at ca. - mya. all cpdna sequences of r. pseudochrysanthum s.l. coalesced at ca. - mya, while those of r. hyperythrum coalesced at ca. - mya. likewise, the nrdna sequences including r. formosanum, r. hyperythrum, and r. pseudochrysanthum s.l. could be traced back to a common ancestor about - mya. all nrdna sequences of r. pseudochrysanthum s.l. and r. hyperythrum coalesced to a common ancestor ca. - mya. all nrdna sequences of r. pseudochrysanthum s.l. coalesced to a common ancestor ca. - mya, while those of r. hyperythrum coalesced at ca. - mya. to elucidate the possible evolutionary forces for the lack of species monophyly of r. formosanum, r. hyperythrum and r. pseudochrysanthum s.l., we calculated the expected time to coalescence at cpdna and at nrdna, where θ is a population mutation parameter obtained from im, and u is the geometric mean of mutation rates, and compared them with the divergence times obtained from the beast analyses. according to the im estimates of θother mtns = and θyushan =, the expected time to coalescence of r. pseudochrysanthum s.l. ranged from to mya based on cpdna, and from to mya based on nrdna. discussion systematic inconsistency between nuclear and chloroplast dnas in rhododendron pseudochrysanthum s.l inconsistency in the cpdna and nrdna phylogenies was detected in r. formosanum. the ml tree for the cpdna agreed with the results of chamberlain et al. and goetsch et al. that r. formosanum is a monophyletic group relative to r. pseudochrysanthum s.l. in contrast, the ml tree of nrdna revealed that one accession of r. formosanum fell within cluster b with the ingroup species. lack of species monophyly may result from incomplete lineage sorting, which occurs when populations or species recently diverged, thus retaining ancestral genetic polymorphisms, or from genetic exchanges via hybridization and introgression that introduce foreign alleles into a species or population, resulting in a genetic admixture. for given species, the time for a locus to acquire reciprocal monophyly is about × θ/4u at cpdna or × θ/2u at nrdna. willyard et al. suggest that ancestral polymorphisms would likely be shared between species if the coalescence time of the alleles exceeds the divergence time at one locus. in this study, the expected time to coalescence of r. pseudochrysanthum s.l. overlapped the estimates of divergence time based on cpdna, in which monophyly was attained. in contrast, at the nrdna, the expected time for coalescence predated the divergence time. accordingly, nuclear loci would less likely be monophyletic in r. pseudochrysanthum s.l. incomplete lineage sorting alone, apparently, cannot completely explain the paraphyly of r. pseudochrysanthum s.l. interspecific hybridization that disturbs the reciprocal monophyly needs to be aroused based on the non-zero m parameters between r. formosanum-r. pseudochrysanthum s.l. and r. hyperythrum-r. pseudochrysanthum s.l.. furthermore, milne et al. revealed that natural hybridization/introgression is not unusual in rhododendron. thus, interspecific hybridization/introgression would cause the lack of monophyly and prolong the duration of paraphyly, as frequently occurs in most plants, e.g., ixeris, begonia, and ilex. genetic diversity and population structure of r. pseudochrysanthum s.l in this study, we investigated genetic diversity at the cpdna and nrdna loci of r. pseudochrysanthum s.l. high levels of nucleotide diversity were detected at both cpdna and nrdna loci, agreeing with a pattern found in many east asian plants, that is considered to have been shaped by a common glaciation history. during the quaternary glaciations, taiwan, as a refuge, provided shelters for species from a range of biomes and thus harbored much plant biodiversity. examples include cunninghamia konishii, pinus luchuensis ssp. taiwanensis, cycas taitungensis, and cycas revoluta. a potential refugium for r. pseudochrysanthum s.l. in taiwan was inferred to harbor the genetic variations. furthermore, highly diverse environments may also have created a wide spectrum of habitats to accommodate mutations. the populations of the yushan mountain range had the highest levels of genetic diversity in the cpdna. a similar pattern was found in taiwan fir, a dominant species of subalpine forests at 2500- m a.s.l.. this common result is likely to be attributable to a stable population size in the yushan mountain range. the habitat ranges of the yushan mountain range, with the highest peak at m, are wider than those of the other mountain ranges, providing diverse habitats for r. pseudochrysanthum s.l. given a stable, effective population size under diverse habitats, the populations of the yushan mountain range fostered a great number of genetic polymorphisms. in contrast, genetic diversity can be lost by chance in fragmented populations, as occurred in the population of the alishan mountain range. ph the interior dominant haplotype, is shared among the alishan, central, and sheishan mountain ranges. a single population of the alishan mountain range was fixed for the dominant haplotype ph indicating a possible genetic bottleneck due to a founder effect and/or genetic drift during a period when population sizes were small. genetic structure of r. pseudochrysanthum s.l r. pseudochrysanthum s.l. inhabits mountaintops isolated by surrounding low-elevation habitats, forming small and fragmented patches. genetic differentiation would be detected between these isolated populations of different mountain ranges. nevertheless, genetic analysis of the cpdna variation revealed very weak genetic differentiation among the populations within a mountain range. the short time of isolation since the last glacial retreat is seemingly not long enough for differentiating these populations within a mountain range and thus resulted in sharing ancestral genetic polymorphisms. specifically, spatial distribution of cpdna revealed a noteworthy pattern of high differentiation between populations of the yushan mountain range and the other mountain ranges. the divergence time for r. pseudochrysanthum s.l. is about - mya, and given that within r. pseudochrysanthum s.l., the populations from different mountain ranges are more differentiated than the species, then it follows that this divergence time likely reflects an estimated divergence time between mountain populations. molecular dating revealed that populations among different mountain ranges split some - mya, a time range approximating or predating the formation of the central mountain range of taiwan. independent colonizations via individuals with two different subclusters, likely from adjacent mainland over pleistocene, may result in high genetic differentiation within r. pseudochrysanthum s.l. unexpected genetic structuring was detected with populations of the yushan mountain range highly differentiated from other mountain ranges, despite this range lying between the alishan and the central mountain ranges. the single population of the alishan mountain range was fixed at the dominant haplotype ph which was shared by the populations of the central and sheishan mountain ranges. some other haplotypes may have been randomly lost from the alishan population due to possible bottleneck/founder's events. despite the fact that the alishan mountain range is geographically distant to central and sheishan mountain ranges, sharing the dominant ph made these populations clustered together and genetically undifferentiated. furthermore, subsequent glacial retreat and mountain uplifting triggered the geographical isolations and may have thereby enhanced the divergence of populations among different mountain ranges. today, over an interglacial period, rhododendron populations of different mountain ranges represent some isolated patches. given such geographical barriers, seeds are less likely to disperse across mountain ranges. low capabilities of seed dispersal in current environments, plus possible independent colonization events followed by limited historical gene flow, may have together contributed to the high levels of genetic differentiation between populations of the yushan mountain range and the other mountain ranges. in contrast, low levels of genetic differentiation at cpdna detected among the populations of alishan, central and sheishan mountain ranges implied that there may exist historical gene flow during the glacial maxima, as the alpine populations migrated downwards and became connected to each other. compared to the high genetic differentiation among populations between yushan mountain range and the other mountain ranges based on cpdna, nrdna revealed much less differentiation among populations of different mountain ranges. here, biparentally inherited nrdna was used to estimate the extent of pollen flow. in contrast to the limited seed dispersal, pollen of rhododendron species seemingly can move long distances. long-distance pollen dispersal of r. pseudochrysanthum s.l. across the discontinuity of mountain ranges reduced the genetic differentiation. ancestral population size, gene flow, and population dynamics of r. pseudochrysanthum s.l given their different inheritances, nrdna and cpdna markers provide extensive aspects to the understanding of the evolutionary forces shaping fates of r. pseudochrysanthum s.l. combining cpdna and nrdna data, im analysis was used to evaluate the ancestral population size and gene flow. populations of the alishan, central, and sheishan mountain ranges were characterized by a recent population expansion as θother mtns was approximately four times larger than ancestral θa. star-like genealogies of cpdna and negative tajima's d of cpdna and nrdna also supported a scenario of recent demographic expansion in populations of the central and sheishan mountain ranges, a common pattern also occurring in castanopsis carlesii, cunninghamia konishii, and trema dielsiana. in contrast, the value of θyushan approximated ancestral θa, indicating that the effective population size of the yushan mountain range may have remained stable relative to the ancestral population. given low seed dispersal capabilities of rhododendron in today's environments, non-zero migration rates detected among populations likely represented historical gene flow across mountain ranges. here, compared to the gene flow among populations of alishan, central, and sheishan mountain ranges, with m = -, much lower levels of migration rates were detected between populations of the yushan mountain range and other mountain ranges. even so, higher levels of migration rates detected from populations of the yushan mountain range to those of the alishan and sheishan mountain range further indicated asymmetric gene flow among them. in contrast, historical dispersal may have occurred freely among other populations. that is, during glacial maxima, there may exist corridors among the mountain ranges. these corridors allowed frequent gene flow among populations of alishan, central, and sheishan mountain ranges, and occasional gene flow between populations of the yushan mountain range and the other mountain ranges. altogether, populations of these mountain ranges may have been mediated with different levels of historical gene flow during the glacial maxima.  in this study, inconsistent gene genealogies of cpdna and nrdna were detected. in the cpdna tree, agreeing with chamberlain et al., r. pseudochrysanthum s.l. is phylogenetically more closely related to r. hyperythrum than to r. formosanum. in contrast, nrdna failed to identify reciprocal monophyly of ingroup vs. outgroup, indicating possible incomplete lineage sorting and interspecific hybridization/introgression. for r. pseudochrysanthum s.l., the spatial distribution of cpdna revealed a noteworthy pattern of high differentiation between the populations of the yushan mountain range and the populations of the other mountain ranges. molecular dating revealed that the split of the populations at different mountain ranges approximated or predated the formation of the central mountain range of taiwan island, likely suggesting independent colonizations via individuals from previously differentiated populations on the mainland. at the population level, the populations of central, and sheishan mountain ranges may have undergone postglacial demographic expansion, while populations of the yushan mountain range are likely to have remained stable ever since the colonization. in contrast, the single population of the alishan mountain range with a fixed cpdna haplotype may have experienced bottleneck/founder's events. 
for r. pseudochrysanthum s.l., the spatial distribution of the cpdna had a noteworthy pattern showing high genetic differentiation between populations in the yushan mountain range and populations of the other mountain ranges. both incomplete lineage sorting and interspecific hybridization/introgression may have contributed to the lack of monophyly among r. hyperythrum, r. formosana and r. pseudochrysanthum s.l. a complex of incipient species with different degrees of morphological or ecological differentiation provides an ideal model for studying species divergence. rooted at r. hyperythrum and r. formosana, both trees lacked reciprocal monophyly for all members of the complex. independent colonizations, plus low capabilities of seed dispersal in current environments, may have resulted in the genetic differentiation between populations of different mountain ranges. at the population level, the populations of central, and sheishan mountains may have undergone postglacial demographic expansion, while populations of the yushan mountain range are likely to have remained stable ever since the colonization. we examined the phylogeography and the evolutionary history of the rhododendron pseudochrysanthum s. l. systematic inconsistency was detected between gene genealogies of the cpdna and nrdna. in contrast, the single population of the alishan mountain range with a fixed cpdna haplotype may have experienced bottleneck/founder's events.
 a complex of incipient species with different degrees of morphological or ecological differentiation provides an ideal model for studying species divergence. we examined the phylogeography and the evolutionary history of the rhododendron pseudochrysanthum s. l.  systematic inconsistency was detected between gene genealogies of the cpdna and nrdna. rooted at r. hyperythrum and r. formosana, both trees lacked reciprocal monophyly for all members of the complex. for r. pseudochrysanthum s.l., the spatial distribution of the cpdna had a noteworthy pattern showing high genetic differentiation between populations in the yushan mountain range and populations of the other mountain ranges.  both incomplete lineage sorting and interspecific hybridization/introgression may have contributed to the lack of monophyly among r. hyperythrum, r. formosana and r. pseudochrysanthum s.l. independent colonizations, plus low capabilities of seed dispersal in current environments, may have resulted in the genetic differentiation between populations of different mountain ranges. at the population level, the populations of central, and sheishan mountains may have undergone postglacial demographic expansion, while populations of the yushan mountain range are likely to have remained stable ever since the colonization. in contrast, the single population of the alishan mountain range with a fixed cpdna haplotype may have experienced bottleneck/founder's events. 